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Preface 


The  purpose  of  this  study  was  to  examine  the  some  of 
the  problems  encountered  in  high  angle-of  attack  flight, 
namely  wing  rock.  This  was  done  both  analytically, 
through  bifurcation  analysis  and  conputer  simulation,  and 
experimentally,  through  flight  test.  In  all  phases  of 
the  research,  I  am  indebted  to  those  who  have  lent 
invaluable  assistance. 
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without  whose  assistance  this  study  would  not  be 
possible.  He  kept  this  research  on  track  during  its 
beginning  phases.  His  knowledge  of  aircraft  dynamics  and 
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indispensable,  as  were  the  simulation  and  driver  progreuns 
he  lent  me.  I  would  also  like  to  thank  Capt  Jeff  Beck, 
who  developed  both  the  12  and  35  state  driver  models  for 
the  bifurcation  analysis,  for  lending  assistance  when  the 
research  was  snagged  by  a  problem  with  the  driver 
programs  or  with  the  AUTO  software.  I  would  also  like  to 
thank  Capt  Dan  Baumann,  for  his  insights  on  the 
aerodynamic  model  he  developed. 

In  the  flight  test  portion,  I  owe  a  great  deal  of 
gratitude  to  Capt  Paul  Lockhart  and  Capt  Bob  Mlnarik 
(USMC),  who  suffered  through  my  attenpts  at  managing  a 
test,  and  flew  the  majority  of  the  test  missions.  I  know 
both  of  them  found  flying  the  F-15D  difficult.  I  would 
also  like  to  thank  the  3246th  Test  Wing  at  Eglin  AFB  FL, 


and  most  notably  Maj  Mark  "Roman*  Polansky  and  Maj  Randy 
Neville,  who  graciously  loaned  us  the  F-15D,  and  worked 
hard  to  make  this  research  a  success.  Finally,  I  would 
like  to  thank  Maj  Dan  Gleason,  who  guided  me  throughout 
the  two  years  both  AFIT  and  TPS  and  without  whose  support 
this  research  would  never  have  been  conpleted. 
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wing  rock  in  the  F-15  was  examined  both  analytically  and 
experimentally.  Using  a  previously  developed  model  for  the 
F-15,  bifurcation  analysis  and  continuation  techniques  were 
used  to  map  out  the  periodic  wing  rock  solutions  and  the 
equilibrium  solutions  leading  up  to  wing  rock.  This  was 
done  for  four  maneuvers;  a  1  g  stall,  rudder  sweeps, 
constant  bank  turns  and  a  symmetric  pullup.  To  supplement 
this  research,  time  history  simulations  were  used  to 
examine  the  time  history  of  wing  rock.  A  study  of 
stability  derivatives  was  also  done,  to  determine  the 
critical  parameters  in  wing  rock.  Bifurcation  was  also 
used  to  study  candidate  feedback  architectures  used  to 
suppress  wing  rock.  It  was  found  that  feeding  back  roll 
rate  was  effective  in  delaying  wing  rock  onset  and 
suppressing  the  subsequent  motion,  but  this  made  the 
aircraft  more  departure  prone. 

The  results  from  the  1  g  stall,  constant  bank  turns, 
and  the  symmetric  pullup  were  experimentally  tested  through 
flight  test.  Wing  rock  onset  differed  4  degrees  AOA  from 
predicted  in  all  maneuvers  but  the  symmetric  pullup,  where 
the  flight  mach  number  correlated  with  the  computer  model's 
flight  condition.  Wing  rock  was  found  to  be  highly  random 
and  non  periodic,  directly  contradicting  computer 
predictions  and  prior  research. 
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I.  Wing  Rock  Analysis  of  the  F-15 


Throughout  the  history  of  aviation,  man  has  been 
extending  the  performance  capabilities  of  fighter  aircraft. 
Traditionally  this  was  done  by  increasing  the  aircraft's 
speed,  altitude,  range,  and  armament,  all  while  keeping  the 
aircraft's  turning  radius  as  small  as  possible.  While 
flight  envelopes  have  grown,  certain  regions  in  the  flight 
envelope  have  been  off-limits,  namely  the  region  on  the 
lift  curve  approaching  stall.  As  aviation  evolved  -.nto  the 
high  performance  jet  fighters  of  today,  this  region  held 
particular  dangers  for  the  pilots:  aircraft  departed 
violently,  sometimes  rolling  or  yawing  viciously  as  they 
stalled,  leaving  the  pilots  in  unrecoverable  flat  spins. 
Thus,  the  region  near  .stall  became  a  very  dangerous  region 
that  was  best  avoided. 

Despite  the  tremendous  increase  in  speed  brought  about 
by  even  more  advanced  turbojet  engines,  combat  experience 
showed  that  dogfighting  between  two  opposing  aircraft  was 
still  prevalent.  This  was  graphically  borne  out  in  the 
Vietnam  conflict,  where  heavier  and  sluggish  F-4's  and  F- 
105 's  were  forced  to  engage  the  slower,  yet  more  nimble 
MiG-17 's  and  MiG-21 's. 

Dogfighting,  or  Air  Combat  Maneuvering  (ACM)  involves 
often  tight  turns  that  push  the  fighters  operational 
capability.  One  of  the  boundaries  reached  by  fighter 
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aircraft  in  the  ACM  environment  is  the  limit  of  useful 
angle  of  attack  vAOA)  .  Often,  this  is  exceeded  in  ACM, 
with  disastrous  results.  The  aircraft  with  a  wider  AOA 
envelope  gives  the  pilot  an  advantage  over  his  adversary. 

Before  stall  and  departure  are  reached  in  combat 
aircraft  however,  a  variety  of  nonlinear  phenomena  are 
encountered,  leading  to  a  degradation  of  the  aircraft's 
flying  qualities.  The  first  of  the  nonlinear  phenomena 
encountered  is  wing  rock.  Wing  rock  is  defined  as  a 
lightly  dcimped  rolling  motion,  and  is  exhibited  in  many 
combat  aircraft  as  angle  of  attack  is  increased  (1:1). 
Wing  rock  imposes  tracking  limitations  on  aircraft  at  high 
AOA,  and  can  cause  other  safety  and  flying  qualities 
problems.  For  example,  wing  rock  can  be  easily  excited  by 
lateral  control  inputs,  and  severe  pilot  induced 
oscillations  (PIO)  can  be  experienced,  resulting  in 
controllability  problems  ranging  from  degradation  of 
tactical  effectiveness  during  combat  maneuvers  to  flight 
safety  during  landing  (2:3). 

Wing  rock  may  be  experienced  throughout  the  entire 
velocity  envelope,  but  the  root  causes  of  the  phenomenon 
vary  over  the  flight  envelope.  In  slow,  subsonic,  high  AOA 
flight,  wing  rock  can  be  triggered  by  the  interaction  of 
the  aircraft  forebody  flow  with  the  downstream  components; 
i.e.  non  symmetrical  vortex  shedding  off  the  forebody  and 
impinging  on  the  vertical  stabilizer.  As  Mach  number  is 
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increased  however,  the  onset  AOA  decreases  and  the  flow 


asyrtmetry  changes.  For  example,  in  the  F-5  aircraft  at 
transonic  speeds,  wing  rock  is  driven  by  shock  induced 
separation,  and  is  a  con^letely  wing  dominated  phenomenon. 
This  characteristic  is  not  unique  to  the  F-5,  as  many 
current  combat  aircraft  exhibit  wing  rock  over  a  wide  range 
of  flight  conditions  in  which  the  character  of  the  motions 
and  the  aerodynamics  that  cause  them  vary  significantly 
(2:5)  . 


Overview  of  Investigation 

The  purpose  of  this  investigation  is  to  analytically 
study  the  wing  rock  phenomena  of  the  F-15,  using 
bifurcation  and  continuation  techniques.  The  validity  of 
the  analysis  and  aerodynamic  model  was  then  verified 
through  a  series  of  flight  tests  performed  at  the  Air  Force 
Test  Pilot  School  at  Edwards  AFB  CA.  To  support  the 
objectives,  the  high  AOA  aerodynamic  model  developed  by 
Bc\umann  (10)  is  combined  with  the  control  system  models 
developed  by  Beck  (13)  for  a  realistic  rr.odel  of  the 
aircraft  dyneunics. 

The  F-15  was  picked  for  this  investigation  primarily 
due  to  the  work  previously  done  by  Eeck  and  Baumann  on  the 
control  models  and  aerodynamics  of  the  F-15.  These  models 
lend  themselves  to  a  thorough  analysis  of  the  wing  rock 
phenomenon  because  they  need  only  further  refinement,  and 
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not  a  full  development,  which  can  be  a  research  project  in 
and  of  itself. 

There  are  four  distinct  areas  that  this  investigation 
covers : 

I.  Investigation  of  wing  rock  throughout  the  F-15 
maneuvering  envelope.  This  involves  three 
different  maneuvers:  1-g  approach  to  stall, 
steady  state  turns,  and  steady  pullups.  The  1-g 
approach  to  stall  flight  condition  was 
investigated  in  reference  13,  so  it  provides  a 
good  starting  point  as  well  as  a  check  to  the 
current  investigation.  The  other  two  flight 
conditions  were  chosen  as  they  represent 
maneuvers  commonly  encountered  in  flight,  and 
especially  in  ACM. 

A  specific  flight  test  program  then 
investigate  wing  rock  during  these  maneuvers, 
with  flight  parameters  set  as  close  as  possible 
to  model  parameters,  in  an  effort  to  validate  and 
expose  any  weaknesses  in  the  analytical  work  done 
on  the  phenomena . 

II.  Investigate  the  effects  of  the  F-15  Control 
Au^mentacion  System  (CAS)  on  wing  rock.  There  is 
■evidence  in  reference  13  that  a  slightly  modified 
CAS  system  has  a  profound  effect  on  suppressing 
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wing  rock.  This  is  further  investigated  by 
modeling  the  F-15  with  and  without  its  CAS 
system. 

III.  Investigate  the  changes  in  aerodynamic 
co<*f ficients  as  AOA  is  increased.  It  is  known 
that  roll  danping  is  an  important  parameter  in 
wing  rock  and  the  roll  danping  derivative  is  the 
primary  determinant  of  dutch  roll  damping  at  high 
AOA.  By  observing  what  changes  occur  to  the 
critical  stability  derivatives  as  a  computi’r 
simulator  of  the  aircraft  approaches  and  enters 
wing  rock,  the  dominant  parameters  can  be 
identified,  in  the  hope  that  the  results  will 
yield  additional  insight  into  the  previous 
research  in  this  area,  and  aid  in  selecting  an 
adequate  flight  control  scheme  to  solve  the 
problem. 

IV.  Investigate  control  system  solutions  to  wing 
rock.  This  phase  consists  of  researching  the 
literature  for  existing  high  AOA  flight  control 
data  and  applying  this  data  to  the  F-15  model,  to 
find  candidate  feedback  architectures  that  will 
suppress  the  wing  rock  problem.  Simple  model 
feedback  architectures  will  be  developed  in 
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conjunction  with  the  results  of  section  three  and 
the  literature  search,  to  modify  the  coefficients 
that  are  most  greatly  affected  by  the  limit 
cycles . 

Overview  of  Thesis 

The  following  chapter  reviews  some  of  the  basic 
principles  of  bifurcation  theory,  to  orient  the  reader  with 
the  terminology  and  the  basic  principles  of  this  new 
analysis  method.  Chapter  3  reviews  the  research  done  on 
wing  rock,  covering  some  of  the  basic  aerodynamics  driving 
the  phenomena,  reviewing  past  research,  and  covering  some 
of  the  benefits  of  suppressing  wing  rock.  Chapter  4  goes 
:nto  greater  depth  on  the  structure  of  wing  rock,  using 
corrputer  simulation  tc  analyze  the  dynamics  of  the  F-15B  as 
it  enters  wing  rock.  The  bifurcation  model  is  also 
introduced  as  a  way  to  detect  the  onset  of  wing  rock,  as 
well  as  the  length  and  stability  of  the  periodic  branch 
that  maps  wing  rock.  Furthermore,  flight  test  results  are 
compared  to  both  the  conputer  simulation  of  wing  rock 
dynamics  and  the  prediction  of  wing  rock  onset  point  as 
modeled  by  the  bifurcation  method.  Covered  at  the  end  of 
the  chapter  is  an  analysis  of  some  of  the  more  important 
stability  derivatives  that  previous  researchers  have 
identified  as  important,  to  see  how  they  contribute  to  the 
phenomena  in  the  F-15B.  This  part  of  the  analysis  is  aided 


6 


by  conparing  the  F-15B  to  a  variety  of  other  aircraft  that 
exhibit  the  same  phenomena,  to  see  if  there  is  one 
stability  derivative  dominant  in  wing  rock. 

Chapter  5  examines  wing  rock  in  maneuvering  flight. 
This  expands  on  the  1-g  stall  that  served  as  the  basis  for 
previous  studies  of  wing  rock  by  examining  the  limit  cycle 
solutions  that  occur  while  the  aircraft  is  performing  basic 
maneuvers,  such  as  steady  turns  and  pullups.  Computer 
models  are  coirpared  to  flight  test  results  to  examine  how 
close  both  correlate.  Also  examined  is  when  wing  rock  is 
encountered  as  only  the  rudder  is  varied. 

Chapter  6  covers  various  single  feedback  architectures 
that  were  tried  to  examine  their  effect  on  wing  rock.  Both 
roll  rate  feedback  to  the  aileron  and  differential  tail  are 
examined  as  well  as  angle  of  attack  and  roll  rate  feedback 
to  the  rudder.  Also  examined  is  the  35  state  reduced  order 
model  of  the  F-15  control  augmentation  system  (CAS)  to 
examine  the  effects  of  a  modified  (but  still  realistic) 
F-15  flight  control  system  on  high  angle  of  attack 
aerodynamics . 

Finally  Chapter  7  covers  the  results  of  the  research, 
and  makes  recommendations  as  to  what  directions  future 


research  could  take. 


To  understand  hew  an  aircraft  transitions  from 
equilibrium  flight  to  the  steady  periodic  motion  known  as 
wing  rock,  this  study  has  used  bifurcation  theor^^. 
Bifurcation  theory  is  the  classical  mathematical  discipline 
that  treats  non-linear  phenomena  (20:xi).  This  chapter  is 
to  orient  the  reader  on  some  of  \.he  basic  principles  and 
terminology  of  this  discipline  that  was  applied  to  this 
research.  Most  of  this  chapter  is  derived  from  Seydel 
(20),  which  the  reader  is  referred  to  for  an  excellent 
discussion  on  this  topic.  The  driver  program  used  in  this 
research,  AUTO,  will  also  be  outlined. 


Equilibrium  and  Stability 

The  motion  of  non-time  dependent  systems  can  be 
mathematically  written  as 


where  u  is  the  state  vector.  This  system  would  said  to  be 
in  equilibrium  if 

0=f(u)  (2) 
Physically,  equilibrium  represents  a  system  at  rest  or  in 
uniform  motion,  such  as  an  aircraft  with  no  translational 
or  angular  accelerations,  and  with  its  roll  and  pitch 
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angles  constant.  The  points  that  satisfy  equation  (2)  are 
known  as  •quillbriua  points. 

The  stability  of  a  system  is  determined  by  its 
behavior  near  an  equilibrium  point.  The  system  is 
considered  stable  if  the  response  to  a  small  perturbation 
is  small  as  time  goes  to  infinity.  If  the  response  goes  to 
zero  as  time  goes  to  infinity  then  the  system  is  said  to  be 
asysptctically  stable.  A  system  that  gets  large  over  time 
is  considered  unstable.  A  system  that  does  not  grow  or  go 
to  zero  as  time  goes  to  infinity  is  said  to  be  neutral. 

Stability  of  a  system  can  be  found  by  linearizing  the 
system  around  the  equilibrium  point.  The  stability  of  the 
system  can  then  be  found  by  looking  at  the  eigenvalues  of 
the  resulting  Jacobian  matrix  of  the  linear  system.  The 
system  is  considered  stable  if  the  eigenvalues  reside  in 
the  left-half  plane;  i.e.  if  they  are  negative  or  zero. 
The  nature  of  the  eigenvalues  determines  the  type  of 
equilibrium  point  and  the  behavior  of  trajectories  near  the 
equilibrium  point. 
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Parameter  Dependence 


The  dependence  of  a  system  on  some  parameter  can 
be  found  by  varying  the  parameter  and  finding  any  new 
equilibrium  points.  Equation  (1)  can  be  rewritten  as 


d  ■  f(u,c) 


(3) 


where  c  is  the  control  parameter,  and  is  called  the 
bifurcation  paraaatar.  In  an  aircraft  model,  these 
parameters  would  be  such  things  as  thrust,  weight,  or 
control  surface  deflections.  By  plotting  the  bifurcation 
parameter  versus  a  state  variable  of  interest,  it  is 
possible  to  obtain  a  qualitative  idea  of  the  system's 
dependence  on  this  parameter.  This  plot  is  known  as  a 
bifurcation  diagraa.  An  exan^le  of  this  would  be  an 
aircraft  model  with  elevator  deflection  as  the  bifurcation 
parameter.  As  the  elevator  is  deflected  from  one  stop  to 
another,  the  bifurcation  diagram  provides  a  global  view  of 
the  aircraft's  behavior.  Unfortunately,  bifurcation 
diagrams  only  provide  stability  information,  and  do  not 
give  information  of  the  aircraft's  dynamic  response  over 
time.  For  time  response  information,  the  researcher  still 
has  to  resort  to  more  traditional  methods,  such  as 
numerical  integration  schemes  to  obtain  this  information. 
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Limit  Points 

One  of  the  branch  points  that  occur  when  an  eigenvalue 
of  the  Jacobian  is  zero  is  a  limit  point  (or  turning 


Figure  1.  Bifurcation  Diagram  with  Limit  Point 

point) .  Figure  1  represents  a  bifurcation  diagram  of  a 
limit  point.  Limit  points  are  characterized  by  having  two 
solutions  on  one  side  of  Cq  and  none  on  the  other  side  of 
Co#  where  Co  is  the  value  of  the  parameter  at  the  limit 
point . 

Because  the  limit  point  occurs  as  a  real  eigenvalue 
moves  across  the  imaginary  axis,  it  represents  a  change  in 
stability.  At  least  one  portion  of  the  branch  departing 


the  limit  point  is  unstable.  Typically,  in  bifurcation 
diagrams  stable  portions  are  represented  with  solid  lines 
and  unstable  portions  are  shown  with  dashed  lines. 


Sometimes,  a  branch  loses  stability  at  one  limit 
point,  only  to  become  stable  at  another  limit  point.  This 
is  known  as  hystarasis,  and  an  example  of  this  is  shown  in 
figure  2,  The 


Figure  2.  Lialt  Points  Showing  Byttaraaia 


arrows  indicate  the  path  the  equilibrium  solution  takes  as 
c  is  increased  or  decreased.  Each  limit  point  has  two 
equilibrium  solutions  that  it  can  exist  at,  and  the 
transition  between  these  two  solutions  is  called  a  juaqp. 
Other  types  of  equilibrium  points  can  exist  when  the 
Jacobian  Matrix  has  a  zero  eigenvalue,  but  these  are  beyond 
the  scope  of  this  research.  The  interested  reader  is 
referred  to  (20)  for  further  information. 
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Limit  Cycles  and  Hoof  Bifurcation  Points 

Sometimes,  systems  tend  towards  a  periodic  or  cyclic 
motion  instead  of  an  equilibrium  point.  A  system  described 
by  a  limit  cycle  is  one  that  remains  in  a  cyclic  motion 
with  period  T,  such  that  u(t)  =  u(t+T)  where  u(t)  is  a 
solution  to  equation  (1)  .  Some  typical  limit  cycles  are 
those  describing  "nerve  impulses,  currents  in  electrical 
circuits,  vibrations  in  violin  strings,  the  flutter  of 
panels  and  laser  light.  '  (20:22)  Figure  3  shows  the 

trajectory  of  a  system  described  by  the  Van  der  Pol 
equation. 

A  -  o(i  -  a^)A  *  a  •  0 

(4) 

with 

o  ■  1 

Uj  «  a 

u,  -  di 

The  system  is  leaving  an  unstable  focus  and  approaching  a 
limit  cycle.  Note  that  for  stable  limit  cycles, 
trajectories  approach  from  the  inside  and  the  outside. 
Limit  cycles  are  of  particular  interest  in  this  study  as 
wing  rock  is  described  by  limit  cycles. 
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Figure  3.  Trajectory  of  the  Van  der  Pol  Bguation 

Limit  cycles  are  detected  on  bifurcation  diagrams 
through  Hopf  points.  A  Hopf  point  is  the  point  where  the 
periodic  motion  emanates  from  a  branch  of  stable 
equilibrium.  Hopf  bifurcations  occur  when  a  coitplex 
conjugate  pair  of  eigenvalues  crosses  the  imaginary  axis. 

Figure  4  shows  a  typical  Hopf  bifurcation.  The 
branch  is  represented  by  plotting  the  maximum  value  of  the 
state  attains  during  the  limit  cycle.  Closed  circles 
represent  stable' values  and  open  circles  represent  unstable 
portions.  Other  bifurcations  that  occur  from  limit  cycle 
branches  are  beyond  the  scope  of  this  review. 
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Bifurcation 


AUTO  Software 

The  software  package  used  in  this  research  to  trace 
out  the  equilibrium  branches  and  find  the  limit  cycles  is 
the  progreim  AUTO  written  by  E.  Dodel  (11)  .  From  a  known 
equilibrium  condition  which  satisfies  equation  (3),  Dodel 
uses  pseudo  arc  length  continuation  to  complete  the 
equilibrium  branches  for  each  new  value  of  the  bifurcation 
parameter.  The  pseudo  arc  length  technique  varies  the 
stepsize  along  the  branch  and  using  the  direction  vector 
(u,C)  a  predictor-corrector  method  algorithm  finds  the  next 
solution.  The  predictor-corrector  method  algorithm  used  is 
the  Newton  method.  The  pseudo  arc  length  technique  allows 
the  algorithm  to  be  scaled  so  it  can  conpute  near  and  past 
limit  points  where  the  direction  vector  is  infinite.  AUTO 
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also  incorporates  adaptive  stepsize,  to  ease  computation 
time  when  the  solution  is  converging  rapidly.  If  the 
solution  is  not  converging,  the  stepsize  is  halved  until  a 
minimum  stepsize  is  reached.  If  convergence  still  does  not 
occur,  the  program  will  also  signal  non-convergence. 

Bifurcation  and  limit  points  are  identified  in  AUTO  by 
monitoring  the  Jacobian  matrix  at  each  solution  and 
identifying  sign  changes  in  the  eigenvalues.  Using 
bifurcation  analysis,  these  changes  are  identified  as  limit 
points,  bifurcation  points,  or  Hopf  points.  AUTO  has  the 
capability  to  calculate  branches  emanating  from  bifurcation 
points,  and  compute  limit  cycles  that  emanate  from  Hopf 
bifurcations.  More  information  on  the  capabilities  of  AUTO 
can  be  found  in  the  AUTO  user  manual  (11) . 
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Ill .  Overview  of  Wing  Rock 


As  modern  fighter  aircraft  engage  in  combat,  they  are 
pushed  to  their  operational  limits.  One  limit  commonly 
encountered  during  the  violent  turns  that  occur  in  the  air 
combat  arena  is  the  AOA  limit,  beyond  which  the  aircraft 
begins  to  exhibit  undesirable  nonlinear  behavior, 
potentially  leading  to  departure  or  spin.  To  gain  a 
maneuverability  advantage,  many  modern  aircraft  designers 
are  trying  to  eliminate  this  undesirable  behavior  at  high 
AOA.  The  first  undesirable  phenomenon  encountered  as  AOA 
is  increased  is  "wing  rock."  It  is  a  common  phenomenon 
found  in  many  low  aspect  ratio  aircraft,  including  such 
aircraft  as  the  F-5,  T-38,  F-14,  F-15  and  the  AV-8B 
Harrier . 

Wing  rock  i*^  defined  to  be  an  uncommanded  rolling 
motion  that  occurs  as  angle  of  attack  is  increased.  It 
builds  up  to  a  limit  cycle,  in  which  the  aircraft  will  roll 
in  one  direction,  stop,  begin  rolling  back  to  the  opposite 
direction  until  a  maximum  roll  angle  is  reached  and  then 
reverse  direction  and  roll  back.  In  studies  done  on  the  F- 
14,  this  limit  cycle  developed  with  amplitudes  on  the  order 
of  ±35®/sec  roll  rate  and  +.10®  of  ride  slip.  This  limit 
cycle  had  a  period  of  "4.5  sec,  indicating  that  it  was  seen 
and  felt  by  the  pilot  (1:150). 
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(1:164) 

IR  indicated  by  figure  5,  wing  rock  onset  occurs  well 
before  reached,  effectively  limiting  the  useful 
AOA  of  the  aircraft.  Wing  rock  occurs  at  a  lower  AOA  than 
the  other  non-linear  phenomena,  and  the  AOA  of  wing  rock 
onset  is  often  the  indication  that  maximum  useful  lift  has 
been  reached  (2:3). 

One  problem  brought  on  by  wing  rock  is  the  loss  of 
tracking  performance,  with  possible  PIO  encountered  as  the 
pilot  tries  to  track  the  target  in  the  region  where  wing 
rock  exists.  An  example  of  this  is  outlined  in  reference 
1,  which  describes  an  F-14  in  ACM.  As  the  pilot  pulled  up 
into  the  wing  rock  region,  he  encountered  PIO  as  he 
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attempted  to  maintain  the  target  aircraft  in  his  gunsight. 
Although  the  anplitudes  were  not  large,  they  were  enough  to 
measurably  degrade  tracking  performance  (1:158). 

Wing  rock  can  also  cause  flight  safety  problems  in 
other  regions  of  flight.  If  it  is  encountered  in  landing, 
the  potentially  high  roll  rates  (as  high  as  S0°/sec)  and 
large  bank  angles  (±90°)  could  lead  to  severe  PIO  close  to 
the  ground  (3:3).  Furthermore,  the  existence  of  wing  rock 
can  cause  problems  in  flight  control  systems  that  are 
designed  for  enhancement  of  flying  qualities  at  high  AOA. 
For  example,  use  of  a  high  gain  aileron-to-rudder 
interconnect  for  elimination  of  adverse  yaw  at  high  AOA  can 
aggravate  and  amplify  the  wing  rock  motions  to  the  extent 
that  the  airplane/pilot  combination  becomes  tactically 
ineffective  (3:3). 

Wing  rock  does  have  one  advantage  in  certain 
configurations;  it  serves  as  a  natural  stall  warning 
mechanism.  Wing  rock  is  the  first  phenomenon  encountered 
as  AOA  is  increased,  and  as  previously  indicated,  onset  of 
wing  rock  is  the  indication  that  maximum  useful  lift  has 
been  reached.  An  excellent  example  of  this  is  the  T-38A 
aircraft.  The  buffet  boundary  proceeding  the  high  sink 
rate  characteristic  of  stall  in  this  aircraft  is  too  wide 
to  use  as  a  stall  warning.  However,  as  the  aircraft  gets 
close  to  stall,  a  moderate  wing  roll  off  (wing  rock)  and 
lateral  oscillations  occur  just  prior  to  the  development  of 
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the  high  sink  rate  with  full  back  stick.  The  wing  rock 
lateral  oscillations  provide  a  satisfactory  stall  warning, 
in  both  cruise  and  landing  configuration  in  this  aircraft 
(4:11)  . 

Despite  this  potential  benefit,  studies  have  shown 
that  the  elimination  of  wing  rock  can  lead  to  greater 
tactical  effectiveness  in  a  high  AOA  ACM  environment. 
Reference  1  describes  one  such  stud^.  In  this  test  two  F- 
14  aircraft,  one  with  a  special  high  AOA  stability  and 
control  system,  and  one  without  were  engaged  in  one-on-one 
ACM.  Despite  various  initial  conditions  of  the  engagement, 
it  eventually  evolved  into  a  low  speed,  high  AOA  turning 
fight  in  which  each  aircraft  was  attempting  to  gain  an 
advantageous  position.  For  con^>ari8on,  these  evaluations 
were  flown  with  and  without  the  stability  and  control 
system. engaged.  The  overall  results  obtained  reconfirmed 
earlier  studies  which  have  consistently  shown  that  a 
properly  designed  high  AOA  control  system  can  provide 
substantial  benefits  during  high  AOA  ACM  in  terms  of 
reduced  pilot  workload,  lower  skill  requirements,  and 
increased  pilot  confidence.  Figure  6  (1:179)  illustrates 
this  improvement.  Shown  are  the  time  histories  from  two 
simulated  engagements  made  with  and  without  the  high  AOA 
control  system.  The  Basic  system  data  indicate  pool 
control  of  roll  and  sideslip  as  the  airplane  was  maneuvered 
above  about  20°  AOA.  Under  very  similar  conditions,  much 


20 


better  control  is  indicated  with  the  system  active  (1:158- 
159) .  Of  course,  this  system  is  not  exclusively  designed 
to  solve  the  wing  rock  problem,  but  by  eliminating  the  wing 
rock  problem  and  iriproving  control  effectiveness  at  high 
AOA,  the  aircraft's  maneuvering  envelope,  and  tactical 
effectiveness,  is  improved. 

Although  wing  rock  is  found  throughout  the  subsonic 
and  transonic  flight  regimes,  the  aerodynamics  driving  the 
phenomenon  vary  in  these  two  velocity  ranges.  Subsonic 
wing  rock,  or  the  wing  rock  that  occurs  in  the  absence  of 
appreciable  compressibility  effects  is  by  far  the  more 
researched  of  the  two  types  of  wing  rock .  The  mechanics  of 
the  phenomenon  are  described  by  Hsu  and  Lan  (5:921): 
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Figure  6.  Coirqparison  of  Airplane  Motions  of  Aircraft  During 


Simulated  ACM  With  and  Without  SCAS  (1:179) 

■Wing  rock  is  an  uncommanded  roll -yaw 
oscillation  dominated  by  roll  motion  oscillating 
with  constant  amplitude.  It  may  be  initiated 
either  with  a  sideslip  or  during  a  zero  sideslip 
flight  with  some  flow  asymmetries  over  the 
aircraft  flying  at  high  angles  of  attack.  Once 
the  asymmetric  flow  starts,  a  rolling  oscillation 
anplitude  will  keep  building  up  if  the  roll 
damping  is  negative.  The  transient  amplitude  of 
wing  rock  will  grow  gradually  over  some 
oscillation  cycles  because  of  roll  instability 
and  negligible  dihedral  effect  at  low  roll  angle. 
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Although  the  roll  danqping  is  negative  at 
small  roll  angles,  it  is  positive  at  larger  roll 
angles  for  a  sustained  wing  rock.  Both  the 
effective  dihedral  effect  and  positive  roll 
damping  via  aerodynamic  non-linearities  at  large 
roll  angles  will  gradually  reduce  roll  rate.  As 
these  restoring  moments  become  stronger,  the 
aircraft  will  reach  a  threshold  roll  angle  and 
finally  switch  the  rolling  direction." 

This  description  of  wing  rock  hits  upon  a  number  of 
in^jortant  points.  First,  it  is  driven  by  flow  asymmetries 
at  high  angle  of  attack.  These  flow  asymmetries  are 
configuration  dependent.  For  aircraft  with  long,  slender 
forebodes,  such  as  the  F-5A  or  X-29,  the  wing  rock  is 
caused  by  asymmetrical  vortex  shedding  from  the  Lorebody 
and  interacting  with  the  downstream  conponents.  For 
aircraft  with  low  to  moderate  sweep,  such  as  the  F-14,  the 
wing  rock  is  driven  by  wing  stall.  Airfoil  sections 
experiencing  leading  edge  stall  cause  a  negative  damping  in 
plunge.  This  will  cause  a  moving  wall  effect  that  promotes 
separation  and  loss  of  lift,  generating  a  rolling  moment 
that  drives  the  motion.  The  delayed  stall  due  to  the 
downstream  moving  wall  effects  on  the  opposite  wing  adds  to 
the  rolling  motion.  Finally,  when  the  roll  angle  is 
increased  to  the  point  where  on  the  down  going 
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wing  half,  the  flow  reattaches,  providing  the  lift  needed 
to  produce  the  restoring  rolling  moment  (6:49)  .  For  a  more 
complete  explanation  of  these  aerodynamic  drivers  of  wing 
roc)c,  the  reader  is  referred  to  Reference  6. 

As  previously  mentioned,  wing  rock  is  an  uncomtnanded 
roll-yaw  oscillation  dominated  by  roll.  There  has  been 
some  thought  that  significant  yaw  occurs  in  wing  rock,  but 
the  roll  observed  is  greater  than  any  yawing  motion.  An 
exanple  of  this  is  found  in  the  T-38A  aircraft.  T-38A 
flight  test  results  (appendix  D)  repeatedly  showed  that 
roll  rate  was  at  least  five  times  the  magnitude  of  yaw 
rates  in  developed  wing  rock  for  all  maneuvers  flown. 
Pilot  observations  correlate  this  data;  no  yaw  was  observed 
in  the  developed  wing  rock  limit  cycles  for  the  T-38A 
aircraft  (18:5) . 

From  this  description  it  is  observed  that  wing  rock 
manifests  itself  in  limit  cycle  oscillations.  After  onset, 
the  oscillations  begin  and  eventually  build  up  to  constant 
airplitude  oscillations.  This  characteristic  of  wing  rock 
lends  itself  to  bifurcation  analysis,  as  the  onset  is 
easily  predicted  by  the  detection  of  Hopf  points  and  the 
calculation  of  the  wing  rock  periodic  branch.  This  will 
comprise  the  majority  of  the  theoretical  research  of  this 
thesis . 
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Conclusions 

Wing  rock  is  a  highly  coitiplex,  non-linear  phenomenon 
that  is  largely  configuration  and  Mach  number  dependent. 
Despite  this,  wing  rock  appears  to  the  pilot  to  be  an 
uncoirenanded  roll,  appearing  to  be  the  same  phenomenon 
throughout  the  flight  envelope,  even  though  it  can  be 
driven  by  radically  different  aerodynamics.  In  air  combat, 
the  aircraft  with  the  lower  angle  of  attack  capability  will 
be  at  a  disadvantage,  with  wing  rock  defining  the  effective 
maximum  lift  available.  Wing  rock  destroys  tracking  and 
can  prove  to  be  a  flight  safety  hazard  in  landing.  Thus, 
it  becomes  inportant  to  mathematically  predict  this 
phenomena,  and  study  what  the  dominant  aircraft  states  are 
so  that  effective  feedback  mechanisms  can  be  developed  to 
suppress  wing  rock. 
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IV.  The  Structure  of  Wing  Rock 


To  fully  examine  the  phenomenon  of  wing  rock,  two 
different  analysis  techniques  have  been  exrqployed  in  this 
study:  bif'.rrcation  analysis  and  computer  simulations  of  the 
aircraft  dynamics.  The  dynamic  simulations  give  the 
researcher  a  means  to  visualize  the  motion  of  the  aircraft 
during  wing  rock,  and  the  bifurcation  analysis  completes 
the  picture  by  providing  a  concise  way  of  finding  the  onset 
point,  the  envelope  of  stable  solutions  and  the  amplitude 
of  the  periodic  solutions.  Bifurcation  analysis  also  more 
easily  lends  itself  to  analysis  of  feedback  architectures 
to  suppress  wing  rock. 

The  aircraft  being  analyzed  in  this  study  is  the  F- 
15B,  the  two  place,  air  superiority  fighter  currently  in 
the  Air  Force  inventory.  The  following  configuration  was 
put  into  the  aircraft  model: 

1 .  The  aircraft  was  modeled  in  cruise 
configuration  with  landing  gear,  flaps 
up  and  speedbrake  in. 

2.  No  stores,  pylons,  or  conformal 
fuel  tanks  were  modeled  on  the 
aircraft . 
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3 .  There  was  no  external  test 
equipment  modeled  on  the  aircraft, 
such  as  test  pitot  booms  or  spin 
chutes . 

4.  Aircraft  Command  Argumentation 
Systems  (CAS)  was  modeled  as  being 
off. 


5.  The''e  were  no  fuel  asymmetries 
modeled. 

A  full  description  of  the  aircraft,  complete  with 
physical  dimensions,  is  given  in  Attachment  A. 

The  aerodynamic  model  eitployed  in  this  study  was 
developed  by  Baumann,  and  represents  the  aircraft  flying  at 
M=0.6,  20,000  ft  pressure  altitude  (PA).  The  aerodynamic 
model  was  formed  by  fitting  the  values  of  the  aircraft 
stability  derivatives  to  polynomials,  so  that  the 
aerodynamics  formed  a  smooth  curve  over  the  entire  AOA 
region.  This  was  done  because  the  bifurcation  driver 
program,  AUTO,  cannot  wor)c  with  functions  or  derivatives  of 
functions  that  have  discontinuities.  A  more  conplete 
description  of  the  aerodynamic  model  can  be  found  in  (10) . 

For  the  simulation,  these  aerodynamic  derivative 
approximations  were  used  with  a  numerical  integration 
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scheme  to  obtain  the  full  aircraft  dynamic  simulator.  The 
aerodyneunics  were  fed  into  an  eight  state  dynaunic  model, 
and  integrated  in  time  to  produce  the  aircraft  dynamics. 
The  dynamics  are  changed  from  the  initial  aircraft  steady 
state  by  increasing  the  aircraft's  elevator  deflection 
steadily  by  simulating  an  infinitely  slow  pull  on  the  stick 
by  the  pilot.  At  -29  degrees,  the  elevator  hits  its  lower 
stop  and  stays  at  its  final  value.  This  decrease  in 
elevator  deflection  causes  an  increase  in  AOA,  and  when  the 
simulator  hits  the  critical  onset  value,  the  dynamics 
simulate  wing  rock. 

The  bifurcation  analysis  was  run  using  the  AUTO 
bifurcation  software  package  developed  by  E.  Dodel  (11)  . 
This  was  run  with  a  driver  program  developed  by  Beck  (13), 
which  was  a  12  state  model  representing  the  dynamics  of  the 
F-15.  This  was  generated  from  an  expansion  of  Newton's  laws 
of  motion.  By  assximing  a  rigid  aircraft,  constant  air 
density,  constant  gravity  and  a  flat  earth,  and  by 
transforming  the  force  and  moment  equations  from  the 
inertial  frame  to  the  body  axis  frame,  the  order  of  the 
model  is  reduced  to  nine.  By  assuming  the  xz  plane  is  a 
plane  of  symmetry,  the  following  equations  are  formed  (Ref 
24) 
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Equation  (13)  was  eliminated,  since  the  other  equations  of 
motion  did  not  depend  on  it.  A  12  state  model  was 
developed,  using  these  equations.  Eight  of  the  states  were 
dynamic  parameters  of  the  F-15  (««,  p, p, q,  r,  TAS,  6,  (j))  and  the 
last  four  states  represented  the  control  surface 
deflections:  aileron,  rudder,  stabilator  and  differential 
stabilator.  While  the  dynamics  of  the  12  state  model  are 
identical  to  the  eight  state  model  developed  by  Barth  (12), 
the  12  state  model  proved  to  be  easier  to  use  when  running 
simple  feedback  experiments.  These  experiments,  as  well  as 
the  more  complicated  thirty  five  state  model  developed  by 
Beck  (13)  will  be  discussed  later. 

The  main  difference  between  the  12  state  model  used  in 
this  study  and  the  12  state  model  used  in  reference  13  is 
the  incorporation  of  the  Baumann  aerodynamic  model.  This  is 
the  same  aerodynamic  model  of  the  F-15  used  in  the  dynamic 
simulator,  and  it  represents  a  significant  improvement  over 
the  aerodynamic  model  developed  by  Barth  (12) . 

Flight  Test  Program 

A  flight  test  program  was  developed  at  the  Air  Force 
Test  Pilot  School  to  experimentally  verify  the  bifurcation 


30 


program  and  coitputer  dynamic  simulation.  For  details  of 
the  test  plan  see  (21)  .  The  tests  erployed  an 
uninstrumented  F-15D,  flown  in  cruise  configuration,  with 
no  stores,  pylons,  or  conformal  fuel  tanks.  The  aircraft 
was  flown  as  close  as  possible  to  the  conditions  modeled  in 
the  bifurcation  analysis.  The  flight  test  techniques  used 
are  found  in  appendices  F,  and  the  tolerance  on  each  test 
point  is  found  in  appendix  G.  Three  basic  maneuvers  were 
flown,  a  1  g  stall,  steady  banked  turns,  and  a  symmetric 
pullup.  The  results  for  the  1  g  stall  are  below,  and  the 
results  for  the  other  maneuvers  can  be  found  with  their 
theoretical  analysis  in  the  following  chapter. 

One  note  about  the  flight  test.  The  F-15D  AOA  gauge 
is  not  calibrated  in  degrees,  but  is  calibrated  in  units. 
The  exact  calibration  of  the  gauge  in  the  test  aircraft  is 
unknown,  so  all  flight  test  results  are  reported  in  units. 
To  compare  the  analytical  results  to  the  flight  test 
results,  the  calibration  between  the  production  AOA  gauge 
and  the  pitot  boom  AOA  gauge  found  in  reference  7  (7:45) 

was  used  before  the  flight  test  progreun  was  initiated  to 
obtain  predictions  on  wing  rock  10°  difference.  In  the  body 
of  the  text,  computed  AOA  is  presented  in  degrees  with  the 
corresponding  cockpit  units  following  in  parentheses.  All 
gr.aphs  are  left  in  degrees,  due  to  the  fact  that  some  of 
the  scales  exceed  the  range  of  the  F-15  production  gauge. 
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and  it  would  be  inaccurate  to  extrapolate  the  calibration 
beyond  the  gr^ges  upper  limit  of  45  units  AOA. 

1  g  Stall  Results 

The  1  g  stall  was  the  first  maneuver  studied,  and  it 
provides  a  good  comparison  to  prior  research  (13) . 
Bifurcation  analysis  was  initially  used  to  plot  the  stable 
equilibrium  branch  of  the  aircraft  motion,  and  detect  the 
region  of  wing  rock.  At  20  degrees  (30  units)  AOA,  a  Hopf 
bifurcation  point  was  detected,  signalling  the  onset  of 
wing  rock.  The  periodic  branch,  which  is  represented  in 
figure  7  by  the  circles  shows  the  trace  of  the  wing  rock 
branch.  Each  circle  on  the  periodic  branch  represents  the 
maximxun  amplitude  in  AOA  that  the  periodic  solution  reaches 
on  that  cycle.  This  representation  also  shows  the 
stability  limits  of  the  branch.  Beyond  an  elevator 
deflection  of  -27  degrees  the  periodic  solution  is 
unstable.  Since  unstable  solutions  are  'sically 
unrealiazible,  the  solution  "jumps*  to  a  higher  t»olution, 
which  in  this  case  is  a  flat  spin.  The  flat  spin  results 
are  not  presented  here,  mcstly  because  the  bifurcation  runs 
concentrated  on  AOA  <  60  degrees.  The  spin  phenomena  is 
explored  in  greater  detail  in  references  10  and  14,  which 
deal  extensively  with  this  phenomena  in  the  F-15B. 

At  the  onset  of  wing  rock,  the  equilibrium  branch 
changes  stability,  with  no  equilibrium  solutions  exiscing 
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Figure  7.  Bifurcation  Diagram  of  F-15  12  State  Model, 
En^hasizing  the  Periodic  Solution 


at  higher  angles  of  attack  for  the  F-15B.  Thus,  the 
bifurcation  analysis  demonstrates  some  previously  well 
known  facets  of  wing  rock.  At  flight  conditions  up  to  20 
degrees  (30  units)  AOA,  the  aircraft  can  be  flown  with 
in^)unity.  Beyond  20  degrees  (30  units)  AOA  non-linear 
dynamics  make  it  inpossible  to  fly  without  degrading  the 
tracking  performance  of  the  aircraft.  And,  like  most 
modern  fighter  aircraft,  it  is  the  first  non-linear 
phenomena  encountered  as  the  angle  of  attack  is  increased. 

An  attenpt  was  made  to  validate  the  above  analysis  and 
validate  it  through  flight  test  using  an  F-15D. 
Physically,  the  F-15D  is  no  different  than  the  F-15B, 
except  it  is  capable  of  carrying  more  internal  fuel  and  has 
inproved  engines.  The  differences  between  the  tvfo  aircraft 
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are  outlined  in  appendix  A.  The  F-15D  was  flown  in  cruise 
configuration,  with  landing  gear  and  flaps  up,  and 
speedbrake  in.  The  aircraft  was  flown  at  20,000  ft  PA,  but 
not  at  M=0.6,  as  it  was  realized  that  the  aircraft  would 
slow  down  as  it  approached  stall.  However,  it  was  assumed 
that  the  aerodynamic  differences  between  M=0.6  and  the  mach 
number  where  the  aircraft  stalled  would  be  small  enough  to 
neglect  and  that  there  would  be  no  mach  effects  on  any  wing 
rock  characteristics. 

It  was  discovered  during  this  flight  test  that  the 
onset  of  wing  rock  occurred  consistently  four  units  below 
what  was  predicted  by  the  bifurcation  model.  This  was 
attributed  to  the  difference  between  the  aerodynamic 
coefficients  at  MsO.6,  where  the  computer  model  was  valid, 
and  the  mach  number  where  wing  rock  onset  occurred.  The 
model  aerodynamics  were  only  valid  at  MsO.6,  and  the 
aircraft  exhibited  wing  rock  at  lower  mach  numbers  (around 
M=0.35)  revealing  a  potential  flaw  in  the  bifurcation 
model.  While  the  eight  state  equations  were  capable  of 
accurately  modeling  the  aircraft  bleed  rate  during  a  1  g 
stall,  the  aerodynaunics  were  still  only  valid  for  M=0.6. 
As  a  result,  the  stability  nap  generated  by  bifurcation  in 
a  1  g  stall  was  not  physically  realizable.  Flight  test, 
however  did  correlate  with  some  results  of  the  theoretical 
analysis,  and  gave  some  valuable  insight  into  wing  rock. 
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Despite  the  wing  rock  being  primarily  roll,  the 
presence  of  as  much  as  8  to  10  degrees  of  yaw  in  some  wing 
rock  motion  has  led  researchers  to  classify  wing  rock  as  a 
"dutch  roll*  type  motion  (2:2;  19:366)  and  is  often 
described  by  pilots  as  such  (16:5;  7:17) .  However,  in  the 
F-15B/D,  wing  rock  is  primarily  a  rolling  motion.  This  is 
evident  in  figure  8,  which  plots  both 

roll  and  yaw  rate  for  the  F-15B  versus  time  based  on 
numerical  simulation.  In  the  plot,  roll  rate  builds  up 
much  faster  than  yaw  rate,  and  achieves  much  higher 
magnitudes.  Ihis  was  further  verified  during  flight 
testing  in  both  the  F-15D  and  the  T-38A  aircraft.  The 
motion  was  primarily  a  rolling  motion,  with  no  yaw  observed 
in  the  F-15D  and  little  yaw  recorded  in  the  T-38A  (  see 
appendix  D)  .  As  a  result,  the  motion  in  these  two  aircraft 
did  not  exhibit  the  classic  'dutch  roll*  motion. 

Also  evider*"  in  figure  8  is  the  "limit  cycle*  behavior 
of  the  wing  rock. 

This  limit  cycle  behavior  was  not  apparent  in  flight 
test,  however.  Bifurcation  did  accurately  detect  an 
unsteady  motion,  but  flight  test  showed  it  was  not 
periodic.  Flight  test  data  revealed  developed  wing  rock 
limit  cycles  to  be  erratic,  non -symmetric,  and 
inconsistent .  A  review  of  qualitative  pilot  comments  on 
the  established  wing  rock  limit  cycle  confirmed  that  the 
magnitude  of  roll  oscillations  about  an  equilibrium  AOA 
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Figure  8.  Roll  and  Yaw  Rates  for  the  F-15  at  the  Onset  of 
Wing  Rock 

were  not  constant .  The  F-15D  would  oscillate  at  some  roll 
angle  about  the  equilibrixjm  angle  of  bank  (AOB) ,  and  if 
allowed  to  persist,  would  eventually  change  to  a  new 
eguilibrum  AOB.  The  direction  of  this  magnitude  change  was 
also  random.  Pilot  comments  included  observations  such  as, 
"...aircraft  is  rocking  back  and  forth,  and  direction  it 
wishes.."  and  "...initial  symmetrical  wing  rocking  motion 
did  not  last  long.  As  AOA  reached  30  units,  the  motion  grew 
in  amplitude,  remained  low  in  frequency,  and  became 
unpredictable"  (18:7). 

This  unpredictability  of  the  motion  made  it  impossible 
to  calculate  period,  frequency,  and  in  some  cases, 
anplitude.  The  amplitudes  recorded  were  found  in  appendix 
E.  Often  the  wing  rock  motion  experienced  more  than  one 
amplitude  at  a  fixed  angle  of  attack;  in  one  case  during 
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the  1  g  stall  (Table  El)  four  different  amplitudes  were 
observed,  with  ban)c  angle  oscillations  ranging  from  ^5 
degrees  to  ^20  degrees. 

As  previously  mentioned,  this  randomness  of  the  motion 
was  not  predicted  by  either  bifurcation  or  the  dynamic 
simulator.  Nor  was  it  expected  from  prior  research. 
Numerous  sources  (1;2;5;15)  cite  that  wing  roc)c  is  a  limit 
cycle  oscillation.  Hsu  and  Lan,  quoted  previously, 
describe  wing  roclc  as  *  a  roll  motion  with  a  constant 
anplitude*  (5:921).  This  behavior  is  found  in  the  Folland 
Gnat,  whose  wing  roc)c  "  can  become  a  sustained  oscillation 
of  near  constant  anplitude,  which  can  be  held  for  several 
cycles.*  (15:5)  Also,  this  behavior  is  found  in  the  F-5 
aircraft  (2:2,10),  which  is  similar  to  the  T-38A. 

Stability  Derivative  Study 

The  F-15B  simulator  program  was  used  to  calculate 
values  of  the  aerodynamic  stability  coefficients  as  a 
function  of  angle  of  attac)c,  in  an  attenqpt  to  correlate 
these  derivatives  with  wing  rock.  Previously,  att€n?)ts 
have  been  made  by  many  researchers 

(1;2;15)  to  correlate  wing  roc)«  with  one  or  more  stability 
derivatives,  particularly  roll  damping,  Ci  .  The  three 
stability  derivatives  most  extensively  studied  are  Ci 
(dihedral  effect),  (directional  stability),  and  Ci  (roll 
damping) . 
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Of  the  three  aforementioned  stability  derivatives, 
roll  danping  has  been  most  closely  linked  with  wing  rock. 
In  a  large  number  of  aircraft,  such  as  the  F-14A  (1),  the 
Folland  Gnat  (15),  and  the  static  model  of  the  X-29A  (16) 
the  loss  of  roll  darping  is  seen  as  a  main  contributor  to 
the  wing  rock  phenomenon.  In  the  X-29  and  Gnat  aircraft, 
roll  damping  changes  from  positive  to  negative  damping  (a 
sign  change  from  negative  to  positive)  and  remains  negative 
throughout  the  wing  rock  region,  whereas  in  the  F-14A  the 
roll  daiping  remains  stable,  but  gets  very  small  in  the 
wing  rock  region.  This  has  led  many  researchers  to 
postulate  that  a  loss  of  roll  damping  is  an  indication  of 
wing  rock.  (2; 5) 

Figure  9  shows  the  F-15B  simulator  mathematical  model 
of  roll  danping  versus  angle  of  attack.  This  figure 
represents  a  curve  fit  of  tabular  data  at  M=0.6,  20,000  ft 
PA.  As  evident  by  this  figure  the  idea  that  roll  damping 
hits  a  critical  low  value  at  the  onset  of  wing  rock  does 
not  apply  to  the  F-15B,  as  roll  damping  never  changes  sign 
(loses  stability)  and  decreases  significantly  long  before 
the  onset  of  wing  rock.  Following  the  hypothesis  that  the 
roll  darping  hits  a  critical  value  upon  onset  of  wing  rock, 
would  lead  one  to  believe  that  onset  would  occur  at  an  AOA 
of  15  degrees  (15  units) .  Flight  test  (Table  El)  indicates 
onset  at  16  degrees  (26  units)  AOA  for  a  1  g  stall  at 
M=0.6,  20,000  ft  PA.  Similar  mathematical  studies  on  the 
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Figur-3  9.  Roll  D-rinping  vs  AOA 

I 

F-4J  "^Iso  exhibit  King  rock  without  drastic  loss  in  roll 
deut^inn  (16:14)  .  In  fact,  the  matnematical  model  for  roll 
danping  in  the  F-4J  is  increasing  throughout  the  wing  rock 
region,  and  while  it  hits  a  minimum  at  onset  (onset 
occurring  at  approximately  18  degrees  AOA)  it  is  still 
fairly  large,  (-.1)  and  since  roll  dan^jing  is  increasing 
with  AOA,  wing  rock  should  die  out .  This  does  not  happen 
in  flight,  and  wing  rock  occurs  until  the  departure  angle 
of  attack  is  reached.  This  concludes  that  roll  dancing 
alone  cannot  be  used  as  an  indication  of  wing  rock,  for 
while  a  loss  of  roll  danping  has  been  associated  with  wing 
rock  in  some  configurations,  a  similar  loss  of  roll  danping 
in  the  F-15B/D  and  F-4J  cannot  be  linked  to  the  wing  rock 
phenomena.  In  the  F-15B/D,  wing  rock  occurs  some  1  unit 
after  a  local  maximum  in  roll  damping  and  occurs  at  values 
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much  lower  than  that  maxima.  In  the  F-4J,  wing  rock  occurs 
at  a  local  maxima  in  the  roll  damping  derivative,  but  the 
increase  in  roll  danping  throughout  the  wing  rock  region 
would  intuitively  lead  to  a  damping  out  of  this  motion, 
something  not  demonstrated  in  flight  test  (16:6). 
Therefore,  aircraft  with  apparently  good  roll  damping  still 
exhibit  wing  rock. 

While  some  aircraft  exhibit  wing  rock  with  apparently 
benign  roll  damping,  it  is  well  documented  that 
artificially  increasing  roll  damping  through  feedback  or 
arbitrarily  increasing  the  derivative,  as  in  reference  16 
and  17,  eliminates  the  wing  rock.  This  aspect  of  wing  rock 
will  be  discussed  later. 

Figure  10  is  the  dihedral  effect,  versus  angle  of 
attack  for  the  F-i5B.  Note  that  at  the  onset  of  wing  rock, 
20  degrees,  dihedral  effect  is  -0.0070. 

Figure  11  shows  versus  angle  of  attack  for  the  F- 
15B.  Figure  11  shows  a  decreasing  trend  and  negative 
directional  stability  beyond  19  degrees  (29  units)  angle  of 
attack. 

High  dihedral  and  a  loss  of  directional  stability  are 
two  of  the  reasons  that  the  wing  rock  is  primarily  a 
rolling  motion.  Reference  16,  which  details  the 
aerodynamic  model  of  the  F-4J,  notes  that  passes  through 
zero  at  about  20  degrees  AOA,  while  at  the  same  point  Cip  is 
still  relatively  large.  This  is  credited  with  making  the 
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Figure  10.  Dihedral  Effect  vs  AOA 

dutch  roll  motion  largely  a  rolling  motion,  namely  wine 
rock  (16:43) .  Reference  1  says  the  same  thing  for  the  F- 
14A.  The  high  level  of  dihedral  effect,  combined  with 
essentially  no  directional  stability  and  further  combinec 
with  a  high  ratio  of  Iz/Ix  results  in  the  oscillation  beinc 
almost  entirely  about  the  roll  axis  with  very  little  yav 
motion  (1:150).  Since  the  F-15B  has  similar 
characteristics,  it  is  expected  that  the  F-15B's  wing  rock 
would  be  primarily  roll,  which  is  what  was  found  in  the 
simulation  runs.  Also,  like  the  F-14,  both  the  F-15  and  th€ 
F-4J  also  have  high  yaw-to-roll  (I,/I,)  moment  of  inertias 
(5  for  the  F-15  and  6.12  for  the  F-4J)  .  These  factors 
combine  to  make  the  wing  rock  that  occurs  in  the  F~15  a 
rolling  motion.  This  high  dihedral  effect,  combined  with 
a  high  yaw-to-roll  moment  of  inertia,  and  a  directionally 
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unstable  aircraft,  combine  to  make  wing  rock  in  the  F-14A, 
F-4J,  and  the  F-15B  primarily  a  rolling  motion. 

Conclusion 

It  has  been  shown  that  the  wing  rock  experienced  by 
the  F-15  is  primarily  a  rolling  motion,  with  much  smaller 
excursions  in  yaw.  This  is  verified  by  dynamic  simulator 
runs  of  the  F-15B  and  through  flight  test.  During  the  wing 
rock  region,  the  F-15B  aerodynamic  model  shows  the  dihedral 
effect  remaining  high,  with  the  aircraft  becoming 
directionally  unstable  at  the  onset  of  wing  rock.  These 
factors,  combined  with  the  aircraft's  high  yaw-to-roll 
moment  of  inertia,  makes  the  wing  rock  primarily  a  rolling 
motion. 

Bifurcation  analysis  predicted  the  onset  point  of  the 
wing  rock  and  mapped  the  branch  of  the  periodic  solutions. 
This  prediction  of  the  wing  rock  onset  was  4  units  off  from 
the  actual  flight  test  data. 

Finally,  important  stability  derivatives  in  the  F-15B 
were  reviewed  to  examine  their  effects  on  wing  rock.  Roll 
damping,  while  important,  remains  high  when  wing  rock  onset 
occurs.  This  leads  to  the  conclusion  that  an  airframe  with 
good  roll  damping  can  still  exhibit  wing  rock,  and  that  an 
aircraft's  wing  rock  susceptibility  cannot  necessarily  be 
explained  by  examining  roll  damping  alone. 
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V.  Wing  Rock  During  Maneuvering  Flight 


Up  to  this  point,  the  wing  rock  that  has  been  studied 
has  concentrated  the  airplane  flying  in  a  steady,  straight 
and  level  flight.  The  elevator  is  then  increased,  and  the 
aircraft  approaches  a  1-g  stall.  This  is  useful  in 
studying  the  wing  rock  phenomena,  but  it  is  a  flight 
condition  that  is  seldom  encountered  in  the  combat  arena, 
as  it  leaves  the  aircraft  vulnerable  to  opposition.  What 
is  most  encountered  in  air  to  air  combat  are  more  dynamic 
maneuvers,  such  as  turns,  and  abrupt  pullups.  It  makes 
some  sense,  therefore,  to  model  the  aircraft  in  these 
maneuvers,  to  see  if  wing  rock  occurs,  and  to  check  at  what 
point  it  occurs.  In  this  regard,  the  airplane  was  studied 
in  three  flight  conditions;  a  rudder  sweep,  in  which  the 
elevator  was  held  fixed  and  the  rudder  was  swept  from  one 
side  to  the  other,  steady  turns,  in  which  the  aircraft  is 
initially  in  a  steady  turn,  and  then  the  elevator  is 
increased,  and  a  pullup,  in  which  the  pitch  rate,  q,  is 
non-zero,  and  the  elevator  is  increased  to  increase  the 
pullup.  Bifurcation  diagrams  are  presented  for  each  case, 
mapping  out  the  stability  of  each  flight  condition.  Again, 
flight  test  is  used  to  verify  the  bifurcation  model  in  the 
steady  turns  and  pullup  maneuver. 
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Rudder  Sweeps 


The  first  case  studied  was  a  rudder  sweep.  In  this 
maneuver  the  aircraft  is  fixed  at  a  certain  angle  of 
attack  hy  fixing  the  elevator  deflection.  The  rudder  is 
then  varied,  both  positive  and  negative,  and  the  point  at 
which  wing  rock  is  noted  and  the  periodic  branch  is  mapped 
out.  It  should  be  noted  that  this  is  not  a  practical 
maneuver,  as  it  is  the  classic  way  to  enter  a  spin. 
However,  it  does  serve  as  a  good  starting  point  for  the 
study  of  wing  rock  in  maneuvering  flight.  Four  runs  were 
made,  at  elevator  deflections  of  -16.72,  -  10.255,  -19.12, 
and  one  at  -14.88  degrees  and  are  shown  in  figures  11,12, 
13,  14.  Despite  the  different  initial  elevator 
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Figure  11.  Rudder  Sweep  with  Elevator  Fixed  at  -16.72 
Degrees 
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Figure  12.  Rudder  Sweep  with  Elevator  Fixed  at  -10.255 
Degrees 


Figure  13.  Rudder  Sweep  with  Elevator  Fixed  at  -19,12 
Degrees 


deflections,  the  onset  point  of  the  wing  rock  remains 
relatively  unchanged,  occurring  at  20+,!  degrees  (30+,l 
units)  AOA-  Also,  all  periodic  branches  extend  to  30 
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Figure  14.  Rudder  Sweep  with  Elevator  Fixed  at  -14.88 
Degrees 

degrees  angle  of  attack,  where  they  tenninate.  (The  stable 
solution  then  •juit^js"  to  the  stable  spin  solutions.  For 
clarity,  these  are  not  shown.)  So,  despite  the  different 
initial  elevator  settings  onset  occurs  at  the  same  point 
(20+,l  degrees),  and  the  periodic  branches  extend  over  the 
S2ume  angle  of  attack  range. 

steady  Turns 

While  pure  rudder  sweeps  are  something  that  a  pilot 
wouldn't  normally  do  unless  he  was  entering  a  spin,  steady 
turns  are  coirononly  encountered,  and  the  ability  to  turn 
quickly  is  inportant  in  gaining  an  advantage  over  an 
adversary  in  air  to  air  combat.  Since  optimum  turning 
flight  utilizes  the  maximum  useable  lift,  wing  rock  is 
often  encountered  in  very  steep  turns,  when  maximum  useable 
lift  is  exceeded. 
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Turning  flight  poses  a  problem  in  that  solving  for  the 
starting  point  is  a  great  deal  more  difficult  than 
rectilinear  flight.  Given  that  the  stability  derivatives 
for  one  flight  condition  are  known,  it  is  possible  to  solve 
for  the  three  control  surface  deflections,  angle  of  attack, 
and  sideslip.  However,  since  the  stability  derivatives  in 
the  bifurcation  model  are  functions  of  the  control  surface 
deflections,  angle  of  attack,  and  sideslip,  it  becomes 
impossible  to  solve  for  these  values  using  traditional 
methods.  But,  through  continuation,  it  is  possible  to  find 
a  steady  turn  solution  by  varying  the  control  surface 
deflections.  This  is  how  it  is  done;  first,  the  aircraft, 
which  is  in  rectilinear  flight,  is  given  an  aileron  sweep. 
The  aileron  sweep  begins  to  roll  the  aircraft.  To  get  the 
aircraft  in  a  steady  turn,  a  rudder  sweep  is  then  done,  and 
the  aircraft  enters  a  turning  spiral,  with  each 
continuation  step  increasing  the  bank  angle.  At  selected 
bank  angles,  the  change  in  altitude  is  calculated,  and  if 
the  change  in  altitude  is  less  than  100  ft  for  the  entire 
turn  (which,  at  M=0.6,  can  cover  several  miles)  then  the 
turn  is  taken  to  be  a  steady  turn.  At  this  point,  the 
elevator  deflection  is  then  increased,  to  simulate  the 
pilot  pulling  the  nose  to  follow  a  target,  v/hile  the 
airplane  was  turning.  Three  turns  were  analyzed,  a  shallow 
turn  of  10  degrees  bank  angle,  a  moderate  turn  of  30 
degrees,  and  a  steep  turn  of  70  degrees. 
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One  other  parameter  that  was  varied  was  thrust.  It  was 
found  early  in  the  investigation  that  too  much  thrust  with 
too  little  bank  angle  would  produce  a  climbing  turn,  and 
without  '■nough  thrust,  the  aircraft  would  not  be  able  to 
pull  the  nose  up  to  a  position  where  a  Hopf  point  could  be 
detected.  Thrust  was  not  varied  like  the  other  parameters, 
but  instead  was  set  as  a  constant  before  the  initial 
aileron  sweep.  Using  these  thrust  values,  and  the  initial 
values  of  the  states  at  these  thrust  values,  it  was 
possible  to  achieve  steady  turns  with  different  bank 
angles.  This  allowed  more  of  the  turning  envelope  to  be 
explored,  for  only  one  steady  turn  was  available  at  the 
value  of  thrust  fixed  at  8300  lbs. 

The  next  three  figures  (15,  16,  17)  represent  the 
bifurcation  diagrams  for  the  three  different  turns 
previously  mentioned.  What  is  interesting  is  that  all  three 
onset  points  are  identical,  and  all  are  identical  to  the 
onset  point  for  the  1-g  approach  to  stall.  (20  degrees,  or 
30  units  AOA)  Also,  all  periodic  branches  extend  to  about 
the  same 
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Figure  15.  Bifurcation  Diagram  of  F-15  Starting  from  a 
10  Degree  Banked  Turn 


Figure  16.  Bifurcation  Diagram  of  F-15  Starting  from  a  30 
Degree  Banked  Turn 
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Figure  17.  Bifurcation  Diagram  of  the  F-15  Starting  from 
a  68  Degree  Banked  Turn 

maximum  point,  about  33.5  degrees  (43.5  units),  and  none  of 
the  periodic  branches  reattach  to  the  equilibrium  branch. 
Despite  the  difference  in  the  initial  starting  maneuvers, 
all  wing  rock  branches  are  remarkably  similar,  with  neither 
the  different  bank  angles  or  thrust  values  having  a 
significant  effect  on  the  periodic  branch. 

Flight  test  results  of  the  same  bank  angles  modeled 
through  bifurcation  analysis  revealed  onset  points  4  units 
lower  than  predicted  by  bifurcation.  Like  the  1  g  approach 
to  stall,  wing  rock  onset  for  these  three  points  occurred 
at  a  lower  mach  number  than  was  strictly  valid  for  the 
computer  model.  Table  El  shows  the  flight  test  onset 
points  for  the  four  maneuvers,  and  the  mach  numbers  at 
which  onset  occurred.  The  mach  numbers  varied  from  0.32  to 
0.44,  different  from  the  0.6  mach  modeled  in  the  F-15 
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aerodtynainics .  This  reveals  that  it  may  not  be  possible  to 
extrapolate  the  aerodynamics  of  the  model  down  to  lower 
mach  numbers.  Even  though  the  onset  points  are  identical 
at  each  point,  like  the  bifurcation  model,  no  conclusions 
can  be  drawn  from  this  because  of  the  different  mach 
numbers  that  occurred  during  flight  test  of  the  wing  rock 
onset  points. 

Steady  Pullup 

Another  coirenonly  encountered  flight  condition  is  the 
steady  pullup,  in  which  the  aircraft  pitch  rate,  q  ,  is 
non-zero.  This  flight  profile  is  encountered  in  a  split-s 
maneuver,  common  in  air  to  air  combat.  Also,  since  it  is 
a  quasi-steady  condition,  it  easily  lends  itself  to 
bifurcation  analysis. 

In  this  study,  only  one  pullup  was  studied,  for  it  was 
assumed  that  the  pitch  rate  would  have  no  effect  on  the 
onset  or  development  of  the  wing  rock  cycles .  To  perform 
the  analysis,  though,  a  modification  to  the  12  state  model 
had  to  be  made  for  the  steady  pullup  case  to  work.  This 
change  was  to  neglect  gravity,  for  it  was  found  that 
gravity  had  no  effect  on  the  onset  point  of  the  wing  rock, 
although  it  did  have  an  effect  on  the  amplitudes  in  the 
developed  periodic  branch.  Gravity  was  neglected  because 
with  it  the  equilibriiim  assun^tion  in  the  theta  dot 
equation  is  violated,  making  the  regular  12  state  model 
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unsolvable.  By  removing  gravity,  the  theta  dot  and  phi  dot 
equations  become  uncoupled,  and  theta  becomes  independent, 
since  it  is  multiplied  by  gravity  in  all  the  equations. 
This  approximation  makes  the  equations  solvable  with  the 
non-zero  pitch  rate  because  the  remaining  equations  no 
longer  rely  on  pitch  attitude,  which  is  transient  during 
this  maneuver.  As  a  result,  the  order  of  the  model  is 
reduced  by  two  states. 


Figure  18.  Bifurcation  Diagram  of  Pullup  Maneuver 

Figure  18  shows  the  results  of  the  bifurcation 
analysis  for  the  steady  pullup.  As  seen  on  the  diagram, 
onset  occurs  at  21  degrees  (31  units)  AOA.  The  fact  that 
q,  the  pitch  rate  is  non-zero  does  not  effect  the  onset  of 
wing  rock.  The  amplitude  of  the  periodic  solutions  are  not 
as  pronounced  as  the  1-g  approach  to  stall  or  turn  cases, 
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but  this  is  because  gravity  is  neglected.  In  the  pullup 
case,  the  maximum  AOA  achievable  is  30  degrees  (40  units), 
slight  lower  than  both  turning  flight  and  1  g  stall. 

This  was  the  one  maneuver  that  could  be  flight  tested 
at  M=0.6,  and  as  a  result,  was  the  one  maneuver  in  which 
flight  test  results  for  wing  roc)c  onset  correlated  well 
with  the  prediction  of  bifurcation  analysis.  Flight  test 
results  revealed  onset  at  30  ^  1  units  which  is  what  was 
predicted  within  the  fidelity  of  the  AOA  gauge. 

Conclusions 

It  was  found  from  the  flight  test  that  there  is  a  four 
unit  difference  between  the  flight  test  and  the  bifurcation 
model  for  the  wing  roc)c  onset  point.  This  discrepancy  from 
the  bifurcation  model  to  the  flight  test  is  believed  to  be 
caused  by  the  differing  aerodynamics  at  the  mach  numbers 
modeled  (M=0.6)  and  the  mach  numbers  actually  flown. 

It  should  be  fully  determined  what  is  causing  the 
discrepancy  between  the  flight  test  results  and  the  model. 
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VI.  Flight  Control  Solutions  to  Wing  Rock 


In  previous  research,  two  solutions  to  the  wing  rock 
problem  have  been  explored;  aerotfynamic  solutions  and 
flight  control  solutions.  The  aerodynamic  solutions 
involve  shaping  either  the  forebody  or  leading  edge 
extensions  of  the  aircraft.  Although  this  method  has  proven 
successful  on  such  aircraft  as  the  F-20  Tigershark,  it 
involves  extensive  fluid  dynamics  research,  as  well  as 
potentially  extensive  modifications  to  the  airframe.  Both 
of  these  methods  are  expensive  to  implement,  and  do  not 
easily  lend  themselves  to  computer  simulations.  As  a  result 
this  research  did  not  investigate  aerodynamic  solutions  to 
the  problem,  concentrating  on  the  flight  control  solutions. 

This  study  is  theoretical,  with  no  flight  test  to 
supplement  the  theoretical  research.  This  is  primarily  due 
to  the  inherent  safety  risks  involved  in  implementing  the 
single  control  structures  studied;  there  are  solutions  of 
such  limited  stability  that  once  wing  rock  was  reached,  the 
aircraft  may  depart  with  only  the  slightest  increase  in 
angle  of  attack.  In  actuality,  such  control  systems  would 
employ  a  complex  series  of  washouts,  which  would  prevent 
the  aircraft  from  departing. 
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F-14A  High  AOA  Flight  Control  System 


Before  examining  the  F-15B,  this  chapter  will  examine 
the  research  done  on  the  F-14A  high  angle  of  attack  flight 
control  system  (1)  .  It  serves  as  a  good  starting  point  for 
the  F-15B  in  that  it  incorporates  a  single  feedback 
strategy,  while  also  con^ensating  for  the  control  induced 
departures  that  this  type  of  feedback  strategy  presents. 
Furthermore,  this  system  was  successfully  flight  tested. 

Wing  rock  suppression  was  accomplished  by  feeding  back 
roll  rate  to  the  differential  tail.  This  augmented  roll 
dairping,  which  decreased  significantly  prior  to  the  wing 
rock  region  (1:152)  .  This  has  a  side  benefit  of  quickening 
the  roll  mode,  which  helps  with  the  roll  reversal  problem. 
However,  this  feedback  architecture,  which  employs  a  high 
gain  roll  damper,  has  been  shown  to  be  very  detrimental  to 
departure/spin  characteristics  under  certain  conditions 
(1:153).  To  prevent  departure/spin,  a  series  of  washout 
circuits  were  incorporated  to  prevent  the  roll  damper  being 
employed  in  a  departure  situation.  The  first  washout 
employed  was  Mach  limited,  washing  out  the  roll  danper 
where  it  was  not  needed  (For  M^v.77). 

Another  problem  encountered  was  the  use  of  the  roll 
damper  during  maneuvering  flight  at  high  angles  of 
attack.  A  common  way  to  roll  an  aircraft  at  a  high  angle 
of  attack  is  with  rudder.  With  the  roll  danper  employed, 
however,  the  flight  system  will  apply  opposite  roll  control 
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to  danp  the  motion.  As  a  result,  the  classical  crossed, 
prospin  controls  are  obtained  which  can  cause  a  departure 
to  occur.  To  avoid  this  problem,  a  pedal  fadeout  was 
employed  such  that  the  roll  danper  is  deactivated  for  large 
pedal  inputs  (1:153).  Finally,  a  fadeout  was  enployed  that 
deactivated  the  roll  danger  at  high  angle  of  attack,  («*>35 
degrees),  to  prevent  the  system  from  applying  inputs  which 
could  aggravate  departure  should  one  occur.  As  it  is, 
precise  tracking  is  no  longer  possible  in  this  region,  and 
the  roll  dait?>er  is  not  needed  (1:153). 

F-15B  Flight  Control  System 

Like  the  F-14,  the  F~15  also  has  roll  rate  feedback  to 
the  differential  tail,  but  it  is  washed  out  prior  to  the 
onset  of  wing  rock.  This  is  not  without  merit,  as  it  was 
found  in  the  initial  flight  testing  of  the  F-15  that  the 
roll  danper  was  applying  adverse  (prospin)  control  inputs 
during  both  wing  rock  and  departure;  i.e.  the  differential 
stabilator  feeds  inputs  that  are  opposite  the  roll  rate 
(7:10).  As  a  result,  one  of  the  initial  recommendations 
implemented  was  a  washout  of  the  roll  danger  to  prevent 
such  inputs  that  would  lead  to  spin. 

Also  in  the  F-15  is  an  aileron-rudder  interconnect, 
which  is  capable  of  controlling  the  aircraft  during  wing 
rock,  because  of  the  rudder's  ability  to  roll  the  aircraft 
in  this  region.  Unfortunately,  The  aileron  rudder 
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interconnect  is  connected  to  a  washout  circuit  that  filters 
out  low  frequency  oscillations,  such  as  wing  rock.  As  a 
result  of  the  two  washout  circuits,  the  F-15  has  minimal 
flight  control  system  enployed  at  the  angle  of  attack  where 
the  wing  rock  occurs . 

Overview  of  Flight  Control  Study 

To  study  the  feedback  effects  on  wing  rock,  two 
approaches  are  taken.  The  first  is  to  study  the  effects  of 
siir^jle  feedback  strategies  on  the  12  state  model .  Three 
feedback  strategies  are  en^loyed:  feedback  of  roll  rate  to 
the  aileron,  feedback  of  roll  rate  to  the  differential 
tail,  and  feedback  of  angle  of  attack  and  roll  rate  to  the 
rudder.  The  first  tv/o  methods  are  to  supplement  the  roll 
daitping  for  the  F-15,  and  the  last  feedback  strategy  is  to 
enploy  the  effectiveness  of  the  rudder  at  moderately  high 
angles  of  attack. 

The  other  approach  is  to  study  the  effects  of  a 
reduced  order  model  of  the  actual  control  system  on  the  F- 
15  developed  in  (13).  This  model  is  presented  with  and 
without  the  roll  danper  limiter  engaged,  to  show  what  the 
current  flight  control  system  can  do  and  to  show  what 
realistic  modifications  can  do  to  solve  the  wing  rock 
problem. 
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The  first  sinple  feedback  strategy  studied  was  the 
roll  rate  to  aileron  feedback.  Feedback  to  the  aileron 
makes  sense;  it  is  the  classical  lateral  control,  and  is 
successfully  employed  to  danp  wing  rock  in  the  X-29  (8:6)  . 
The  feedback  strategy  employed  is  -60*K*p,  with  four 
different  gains  studied.  Figure  19  and  20 
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Figure  19.  Roll  Rate  to  Aileron  Bifurcation  Diagram, 
Showing  the  Periodic  (Wing  Rock)  Solutions 


shows  the  effect  of  this  feedback  on  the  periodic  branch. 
(The  equilibrium  branch  leading  up  to  the  Hopf  point  is  not 
shown,  and  for  simplicity  the  traditional  representation  of 
the  bifurcation  analysis  is  abandoned.  In  the  following 
cases  studied,  a  solid  line  indicates  a  stable  periodic 
solution,  a  dashed  line  an  unstable  one.) 
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Figure  20.  Detail  of  Roll  Rate  Feedback  to  Aileron 
Bifurcation  Diagreun 


Two  things  are  immediately  evident:  first,  the  Hopf 
point  is  initially  delayed  with  gain,  and  secondly,  and 
perhaps  more  inportantly,  the  length  of  the  periodic  branch 
is  greatly  reduced  as  gain  is  increased.  For  exaitple,  in 
the  no  feedback  case,  the  stable  wing  rock  region  extends 
from  an  angle  of  attack  of  20  degrees  to  an  angle  of  attack 
of  33  degrees.  But  for  the  case  where  the  feedback  gain 
egi^als  0.25  the  stable  wing  rock  region  extends  only  one 
degree  before  departure.  As  the  gain  is  increased,  the 
wing  rock  region  is  increased,  but  the  onset  point  also 
decreases,  until  the  onset  point  is  back  to  the  no  gain 
case  when  K=2  .  Thus,  this  simple  feedback  strategy  exposes 
the  danger  of  feeding  back  roll  rate  to  suppress  wing  rock; 
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it  greatly  increases  the  chance  of  departure.  This  yields 
mathematical  support  to  the  results  that  a  high  gain  roll 
danper  aggravates  the  divergence  susceptibility  of  the 
aircraft.  It  also  shows  that  for  a  successful  design, 
there  has  to  be  a  tradeoff  between  the  onset  of  wing  rock 
and  the  stability  cf  the  solution.  As  clearly  shown  in 
figure  21,  the  onset  can  be  delayed  up  to  26  degrees  AOA, 


Figure  21.  Position  of  Hopf  (Onset)  Point  versus  Gain 
Aileron  Feedback  Case 

but  the  stable  branch  is  so  small  that  an  effective  washout 
circuit  has  to  be  employed  to  prevent  the  aircraft  from 
diverging.  The  area  of  effective  gain  is  small,  as  beyond 
K=2  the  onset  point  decreases  beyond  the  case  of  no 
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feedback,  with  a  substantial  decrease  in  the  length  of  the 
stable  region  of  the  branch. 

Roll  Rate  to  Differential  Tail 

Figure  22  shows  the  analysis  of  the  roll  rate  feedback 
to  the  differential  tail.  Like  the  roll  rate  to  aileron 
feedback  case,  the  increase  in  gain  initially  delays  the 
onset  of  wing  rock,  but  also  decreases  the  stable  wing  rock 
region.  (For  clarity,  the  zero  feedback  case  is  not  shown.) 
The  feedback  delays  the  onset  point  in  a  similar  fashion  to 
the  roll  rate  to  aileron  feedback  case,  with  the  longest 
delay  occurring  at  K=0.5  (see  Figure  23). 


Differential  Tail 
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Figure  23.  Gain  versus  Onset  Point,  Roll  Rate  to 
Differential  Tail  Feedback 

Likewise,  as  gain  is  increased  beyond  the  optimum  point, 
the  onset  point  decreases  back  to,  and  beyond  the  onset 
point  for  no  feedback. 

It  is  known  that  roll  rate  feedback  delays  the  onset 
point,  but  what  does  the  feedback  do  when  the  aircraft 
enters  wing  rock?  The  answer  is  ttiat  the  lateral  motion 
that  is  characteristic  of  wing  rock  gets  danped  out. 
Figure  24  provides  an  exanplc  of  this,  as  roll  angle  is 
plotted  against  angle  of  attack,  for  a  typical  limit  cycle. 
As  gain  is  increased  the  excursions  in  roll  angle  go  from 
50  degrees  in  the  zero  feedback  case  to  5  degrees  in  the 
K=0.2  feedback  case.  Thus  as  gain  is  increased  the 
excursions  on  both  angle  of  attack  and  roll  angle  are 
greatly  damped.  Therefore,  not  only  is  it  possible  to.  move 


62 


the  onset  point  almost  five  degrees  AOA  as  gain  is 
increased,  but  the  wing  rock  that  does  occur  is  greatly 
damped  over  the  zero  feedback  case. 


Figure  24.  Roll  Angle  versus  Angle  of  Attack  for  Various 
Gains  in  to  the  Differential  Elevator 

Augmenting  the  roll  danqping  in  the  F-15  greu.ly 
suppresses  the  wing  rock  motion,  and  can  give  the  F-15  six 
to  seven  additional  degrees  of  angle  of  attack.  However, 
as  pointed  out  in  (7),  care  must  be  taken  so  that  the  roll 
d2ui^)er  does  not  induce  prospin  inputs  in  the  high  AOA 
region.  As  a  result,  any  flight  control  scheme  to  give  the 
F-15  a  wider  angle  of  attack  envelope  would  have  to  include 
the  necessary  washouts  to  prevent  pro  spin  inputs  and  to 
allow  the  pilot  to  control  the  aircraft. 
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Angle  of  Attack  and  Roll  Rate  Feedback  to  the  Rudder 

As  stated  previously,  the  rudder  is  often  used  to  roll 
the  aircraft  at  high  angles  of  attack.  Thus,  some  sort  of 
feedback  to  the  rudder  is  logical  to  inprove  control  of  the 
aircraft  at  high  AOA.  In  this  analysis  the  feedback 
architecture  used  was  This  was  chosen  as  it 

represents  the  aileron-to-rudder  interconnect  that  is  found 
in  the  F-15.  Unfortunately,  in  the  F-15  this  interconnect 
is  attached  to  a  washout  circuit  (s/s+1)  that  filters  out 
all  low  frequency  inputs,  such  as  wing  rock.  As  a  result, 
the  aileron-to-rudder  interconnect,  which  is  capable  of 
suppressing  wing  rock,  is  not  effective  in  the  wing  rock 
region. 

Figure  25  shows  the  periodic  solutions  for  this  type 
of  feedback  structure.  Like  all  of  the  other  feedback 
structures,  an  increase  in  gain  widens  the  equilibrium 
angle  of  attack  envelope,  but  also  lowers  the  stable  range 
of  solutions.  Figure  26  shows  the  position  of  the  Hopf 

(onset)  point  with  different  gains.  Unlike  the  two 

I 

previous  solutions,  this  feedback  structure  can  accept 
higher  gains. 

I 

I 

I 

I 

j 

i 

Despite  the  increased  gain,  this  feedback  strategy 
basically  gives  the  same  performance  in  moving  the  onset 
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Figure  25.  Wing  Rock  Solutions  for  Angle  of  Attack  and 
Roll  Rate  Feedback  to  the  Rudder 


Figure  26.  Onset  Point  versus  Gain  for  the  Rudder  Feedback 
case 

point  and  daitping  the  wing  rock  motion,  and  like  the  other 
two  strategies,  it  requires  a  series  of  washout,  circuitry 


to  prevent  the  flight  control  system  from  entering  any 
departure  inducing  controls. 

35  State  Model 

In  an  effort  to  model  the  d^amics  of  the  F-15  control 
system,  models  were  developed  (13)  that  included  the 
effects  of  the  F-15  control  system  on  the  airplane 
dynamics.  Three  additional  aircraft  states,  forward 
acceleration,  normal  load  factor,  and  lateral  load  factor 
were  added,  as  well  states  for  the  longitudinal  stick 
dynamics,  the  lateral  stick  (^namics,  rudder  pedal 
dynamics,  aileron-rudder  interconnect,  and  c^amics  for  the 
aircraft  control  system.  Originally,  these  models  included 
all  of  the  servo  and  high  frequency  filters  that  were 
incorporated  in  the  F-15  control  system,  rasing  the  number 
of  states  of  the  model  to  50.  Since  most  of  the  high 
frequency  effects  had  little  effect  on  the  aircraft 
dynamics,  they  were  eliminated,  which  reduced  the  aircraft 
model  to  35  states.  Thus,  a  workable  model  of  the  F-15 
flight  control  system  was  developed,  and  this  model  was 
used  to  study  the  effects  of  the  flight  control  system  on 
the  aircraft  dynamics. 

In  developing  the  model  of  the  control  augmentation 
system,  one  of  the  assunptions  made  was  to  remove  all  the 
limiters  from  the  flight  control  system.  This  was  done 
because  in  limiting  the  continuation,  in^)ortant,  physically 
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realizable  solutions  may  exist  in  a  different  area  of  the 
bifurcation  diagram  and  attached  by  unstable  solutions 
that  exist  outside  the  physical  limits  of  the  airplane.  For 
exarr^le,  if  the  continuation  of  an  elevator  sweep  stopped 
at  the  physical  limit  of  elevator  travel,  -29  degrees,  the 
stable  spin  regions  which  exist  at  extremely  high  AOA,  (and 
are  physically  realizable)  would  be  missed.  This  would 
severely  limit  the  usefulness  of  the  bifurcation  analysis. 

Some  of  the  limiters,  however,  have  an  important 
effect  on  the  aircraft  dynamics.  One  of  the  limiters 
removed  in  both  the  50  and  35  state  models  developed  by 
Beck  was  the  roll  daii?>er  washout.  While  this  limiter  has 
little  effect  on  the  equilibrium  solutions,  (17; 6)  it  has 
a  profound  effect  on  the  periodic  solutions.  When  this 
limiter  is  removed,  the  flight  control  system  provides 
excellent  danping  of  the  wing  rock,  and  possible 
inprovement  in  the  onset  point.  Unfortunately,  this  is  not 
realistic,  for  as  recommended  in  (7),  the  roll  danper  is 
scheduled  out  by  the  onset  of  wing  rock. 

Figure  27  shows  the  washout  schedule  for  the  roll 
damper  in  the  F-15  and  in  the  conputer  model.  To  modify 
the  35  state  model,  this  limiter  was  fit  to  a  fourth  order 
polynomial  for  smoothness  considerations,  and  placed  in  the 
model,  to  wash  out  the  gain  as  angle  of  attack  is 
increased.  After  an  angle  of  attack  of  20.2  degrees,  the 
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gain  is  set  to  zero,  identical  to  the  schedule  found  in  the 
F-15. 


Figure  27.  Washout  Schedule  for  F-15  Roll  Damper 

Figure  28  shows  the  results  of  the  35  state  model, 
with  and  without  the  limiter.  The  results  are  a  bit 
truncated,  due  to  difficulty  with  the  model,  but  there  is 
enough  to  show  the  trends  of  the  analysis.  For  a  more 
complete  run  of  the  35  state  model  with  ae  limiter 
disengaged,  the  reader  is  referred  to  (17)  .  The  results 
from  the  previous  analysis  of  the  simple  feedback 
structures  become  helpful  here.  Similar  to  those  cases, 
when  feedback  is  present,  the  onset  point  is  increased,  and 
the  maximum  amplitude  AOA  excursions  are  greatly  reduced 
for  the  same  longitudinal  stick  force.  With  the  limiter, 
the  periodic  branch  greatly  resembles  the  no  feedback  case. 


68 


with  the  amplitude  staying  initially  flat  and  then  greatly 
increasing.  Also,  wing  rock  occurs  at  a  lesser  stick 
force,  and  onset  occurs  at  a  slightly  lower  angle  of 


attack.  The  inprovement  in  the  onset  point  is  to  be 
expected,  as  in  the  sinple  feedback  cases.  The  difference 
in  shape  and  onset  of  the  periodic  branch  with  the  limiter 
is  due  to  additional  control  system  dynamics  in  the  model. 


Figure  28.  Periodic  Branches  Derived  from  the  35  State 


Model 
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Conclusions 

As  seen  by  the  three  sinple  feedback  cases,  a  fairly 
sinple  feedback  structure  can  expand  the  angle  of  attack 
range  of  the  F-15.  In  all  three  cases  explored,  an  optimum 
gain  was  found  that  raised  the  onset  of  wing  rock  from  20 
degrees  AOA  to  26  degrees  AOA.  Increasing  gain  was  also 
shown  to  damp  the  wing  rock  when  the  oscillations  do  occur. 

The  most  profitable  feedback  structures  involve 
augmenting  the  roll  daitping  derivative.  Feedback 
strategies  involving  augmenting  the  roll  d2UTiping  by  feeding 
back  roll  rate  to  the  differential  tail  have  been 
successfully  flight  tested  on  the  F-14A  (1) .  And,  despite 
the  significantly  higher  roll  damping  the  F-15B  has  over 
the  F-14A  in  the  wing  rock  region,  augmenting  the  roll 
danping  in  the  F-15B  also  proves  beneficial  in  delaying  and 
suppressing  the  wing  rock. 

Despite  the  inprovements  that  these  simple  feedback 
strategies  offer  in  expanding  the  angle  of  attack  envelope, 
a  number  of  washout  circuits  would  have  to  be  employed  in 
order  to  keep  the  roll  damper  from  providing  departure 
inducing  control  inputs  at  higher  angles  of  attack,  and  to 
control  the  aircraft  at  these  flight  conditions.  In  the 
current  condition,  the  F-15B  has  a  circuit  that  washes  out 
the  roll  damper  by  the  onset  of  wing  rock  to  prevent  the 
departure  inducing  control  inputs  from  the  roll  damper. 
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The  35  state  model  of  the  flight  control  system 
employs  a  roll  danper  that  is  modeled  after  the  roll  danper 
found  in  the  F-15.  When  this  is  not  washed  out,  it 
increases  the  onset  point  of  the  wing  rock .  When  the 
washout  circuit  is  engaged,  the  wing  rock  solution  moves 
towards  the  unaugmented  solution. 

There  is  good  reason  for  washing  out  the  roll  damper, 
however,  and  the  bifurcation  diagrams  of  the  siitple 
feedback  structures  bear  this  out.  As  the  gain  is  increased 
from  the  no  feedback  case,  the  stable  periodic  solution 
branch  length  is  greatly  decreased.  This  lends  mathematical 
support  to  the  hypothesis  that  high  gain  roll  dampers  often 
aggravate  divergence.  Thus,  despite  the  advantages  that 
roll  danpers  give,  care  must  be  taken  that  the  feedback 
strategy  eirployed  to  suppress  the  wing  rock  does  not 
aggravate  the  other  non-linear  phenomena  that  occur  in  this 
flight  regime. 
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VII.  Conclusions  and  Recommendations 


Both  the  flic  it  test  analysis  and  computer  model 
revealed  that  wing  rock  in  the  F-15  is  primarily  a  rolling 
motion.  The  aircraft's  high  dihedral  effect,  combined  a 
high  yaw  to  roll  inertia  ratio  and  a  loss  of  directional 
stability  at  wing  rock  onset  make  the  resultant  motion 
primarily  roll. 

Flight  test  correlated  well  with  the  computer  model  of 
the  symmetric  pullup,  with  onset  occurring  at  30  units  AOA 
in  both  cases.  However,  in  all  other  maneuvers  modeled 
using  bifurcation,  the  onset  point  differed  by  4  units  AOA 
from  the  flight  test  results.  This  puts  into  question  the 
assumption  that  the  aerodynamics  modeled  at  M=0.6,  20,000 
ft  PA  is  valid  for  all  mach  numbers  below  M=0.6. 
Furthermore,  prior  research  and  the  conputer  generated  wing 
rock  limit  cycles  showed  symmetric,  periodic  behavior,  none 
of  which  was  observed  in  flight  test.  Wing  rock  was 
characterized  by  non -symmetric,  almost  random  behavior, 
with  constantly  changing  bank  angle  oscillations  and 
equilibrium  bank  angles. 

Computer  modeling  and  analysis  showed  that  a  loss  of 
roll  damping  did  not  necessarily  correspond  to  the  onset  of 
wing  rock,  as  the  F-15B  still  had  significant  roll  damping 
when  onset  occurred  in  the  coitputer  model. 
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Also  examined  was  the  effects  of  a  flight  control 
system  on  the  suppression  of  wing  rock.  It  was  seen  in  the 
three  simple  feedback  cases  studied  that  an  feedback  of 
roll  rate  was  highly  beneficial  in  delaying  the  onset  of 
wing  rock,  as  well  as  damping  out  the  resultant  wing  rock 
oscillations.  However,  the  resultant  stability  of  each 
wing  rock  branch  is  greatly  reduced,  increasing  the 
aircraft's  departure  susceptibility.  In  order  to  be 
implemented,  any  feedback  strategy  studied  must  employ  a 
coirplex  series  of  washout  circuits  to  reduce  this  departure 
susceptibility.  A  reduced  order  model  of  the  F-15  flight 
control  system  was  also  studied,  with  and  without  the  ^roll 
daitper  washout  circuit.  The  35  state  model's  performance 
was  comparable  with  the  unaugmented  modol  with  the  roll 
damper  engaged.  With  the  roll  damper  washout  disengaged 
the  performance  of  the  system  is  greatly  improved. 

Despite  the  failure  of  the  analysis  to  accurately 
predict  wing  rock  onset  in  the  majority  of  cases,  and  that 
the  coioputer  model  and  the  flight  test  occurred  at 
different  points  in  the  envelope,  making  any  results 
difficult  to  compare,  bifurcation  does  hold  the  power  to 
predict  non-linear  phenomena  such  as  wing  lock,  and  show 
the  regions  where  non  linear  behavior  occurs  in  the 
aircraft  flight  envelope.  With  the  following 
recommendations,  it  would  be  possible  to  predict  the  onset 
of  wing  -ock  and  the  developed  limit  cycle  as  done  in  this 
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research,  but  more  in^ortantly,  it  will  be  possible  to 
verify  the  model  at  more  than  one  point  through  flight 
test. 

Recommendations 

1.  Determine  why  there  is  a  4  degree  difference  in 
onset  point  between  the  flight  test  results  and  the 
bifurcation  model  for  the  1  g  stall  and  the  constant  bank 
turns . 

2.  Examine  more  closely  the  effects  of  various 
parameters  on  wing  rock,  such  as  pitch  rate,  g  and  mach 
number  to  get  an  idea  as  to  what  pareuneters  are  important 
in  the  developed  wing  rock  limit  cycle. 

3.  Oetejmine  why  the  wing  rock  motion  in  both  the  T- 
38*1  and  F-15D  is  so  random. 

4.  Examine  the  developed  wing  rock  limit  cycle,  and 
determine  if  the  change  in  equilibrium  bank  angle  during 
wing  rock  has  any  effect  on  the  resultant  motion. 

5-  Rebuild  the  model  with  the  Pitch  CAS  on,  to 
reflect  the  actual  flight  condition.  The  AOA  achieved  in 
this  study  required  the  pitch  CAS  to  be  on  in  flight,  which 
was  not  reflected  in  the  12  state  model.  For  a  more 
accurate,  realistic  model,  the  effects  of  the  pitch  CAS 
should  be  modeled. 
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APPENDIX  A 


test  item  description 


The  F-15B/D  aircraft  was  a  two-place,  supersonic  air 
superiority  fighter.  The  aircraft  was  powered  by  two  Pratt 
and  Whitney  FlOO-PW-lOO  turbofan  engines  rated  at  14,670  lb 
thrust  military  and  23,830  lbs  thrust  in  afterburner.  The 
F'lSB/D  flight  control  system  incorporated  irreversible, 
hydraulically  powered  ailerons,  horizontal  stabilizers,  and 
rudders.  The  horizontal  stabilizers  could  also  be  moved 
differentially,  to  aid  in  rolling  the  aircraft.  Spring 
cartridges  provided  simulated  aerodyneunic  feel  to  the 
pilot.  All  inputs  to  the  actuators  were  provided  by  the 
CAS  which  augmented  the  basic  stability  of  the  airframe  in 
all  three  axes.  This  CAS  system  also  included  an 
aileron-rudder  interconnect  which  served  to  coordinate  the 
ailerons  and  rudder.  The  test  F-15D  flown  was  not  equipped 
with  any  special  instrumentation  and  was  outfitted  with 
standard  production  avionics.  A  full  description  of  the 
aircraft  is  contained  in  the  F-15  A/B/C/D  Flight  Manual 
(22)  . 
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The  physical  dimensions  and  weight  and  balance  data 
for  the  F-15B/0  is  listed  below.  This  data  was  obtained 


from  Beck  (13)  and  (23). 

Physical  Characteristics  of  the  F-15B/D 


Wing 

Area  (Theoretical) 

Aspect  Ratio 

Airfoil 

Root 

Xw  155 

Tip 

Span 

Taper  Ratio 

Root  Chord  (Theoretical) 
Tip  Chord 

Mean  Aerodynamic  Chord 

Leading  Edge  Sweep  Angle 

25%  Chord  Sweep  Angle 

Dihedral 

Incidence 

Twist  at  Tip 

Aileron  Area 

Flap  Area 

Speed  Brake  -  Area 

Control  Surface  Movement 

Aileron 

Speedbrake 

Flap 

Horizontal  Tail 
up 

Rudder 

Vertical  Tail 

Area  (Theoretical  Each) 

Rudder  Area  (Each) 

Span 

Aspect  Ratio 
Root  Chord 
Tip  Chord 
Airfoil  -  Root 
-  Tip 

Taper  Ratio 

Leading  Edge  Sweep  Angle 


61/8  sg  ft 
3.01 

NACA64006.6 

NACA64A(x) 046  (a=0.8  Mod) 
NACA64A203  (a^O.E  Mod) 

42.8  ft 
0.25 

273.3  in 

68.3  in 

191.3 

45  degrees 

38.6  degrees 
1  degree 
None 
None 

26.5  sq  ft 

35.8  sq  ft 

31.5  sq  ft 


■t/~  20  degrees 

45  degrees  up 

30  degreer.  down 

29  degrees  down,  degieas 

f/-  30  degrees 


62.6  fq  ft 
10.0  sq  ft 

10.3  ft 
1.70 

115.0  in 

30.6  in 
NACA0005-64 
NACA0003.5-64 
0.27 

36.6  degrees 


25%  Chord  Sweep  Angle 
Mean  Aerodynaoiic  Chord 
Cant 

Length  (.25c»  to  .2SCh) 


29.7  degrees 
81.0  in 
2  degrees  out 
241.0  in 


Wetted  Area 

Fuselage 

Nozzles 

Horizontal  Tail 
Vertical  Tail 
Wing 

Total  Area 


1405  sq  ft 
51  sq  ft 
216  sq  ft 
257  sq  ft 
698  sq  ft 
2629  sq  ft 


Engine  Data  (each) 

Non  Afterburning  Thrust 
Afterburning  Thrust 
Y  Direction  C.G.  Offset 
Z  Direction  C.G.  Offset 
Nozzle  Pivot  C.G.  Offset 


14,871  lb 
21,810  lb 
♦/-  25.5  in 
0.25  in 
-20.219  ft 


Miscellaneous  Data 

Aircraft  Length  63.8  ft 

Aircraft  Height  18.6  in 

Aircraft  Volume  1996  cu  ft 

Aircraft  Gross  Weight  (Talceoff)  38400  lbs 

C.G.  Station  X  Direction  557,173 

y  Direction  0.0 

Z  Direction  116.173 


Inertial  Data 

I. 

ly 

I. 


25480  slug-ft’ 
166620  slug-ft’ 
186930  slug-ft’ 
-1000  slug-ft’ 


The  inertia  values  are  for  a  basic,  clean  F-15B  with 
ammo,  50%  fuel  and  gear  up. 

Differences  Between  F-15B  and  F-15D 

F-15B  F-15D 

Engines  PW  FlOO-PW-lOO  PW  FlOO-PW-200 


Internal  Fuel  11,500  lbs 


13,455  lbs 


APPSNDIZ  B 


BIPORCXTIOM  DRIVSR  l?RO<3RAM 
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PROGRAM  STATE12 

LAST  EDITED  ON  22  July  1990 


^••••••••*****>**«******* . 

c 

C  Revised  22  July  1990  -(M. Davison)  This  is  a  revision  of  Capt. 

C  Jeff  Beclc's  Model  5  which  consisted  of  the  eight  aircraft 
C  dynamic  states  plus  three  additional  states  for  the 
C  stabilator, rudder,  and  aileron  control  surfaces.  It  also 
C  consisted  of  3  parameters  for  connutnded  control  surface 
C  deflections.  Furthermore,  Model  5  accounted  for  the 
C  differential  tail  deflection,  and  this  was  the  final  state  in 
C  the  12  state  model.  The  main  diffrence  between  this  progreun  an 
C  the  Beck  model  S  is  the  introduction  of  the  in^roved 
C  aerodynamic  model  devloped  by  Capt .  Dan  Baumann 
C  in  1989 .  This  model  is  a  further  refinement  of  the 
C  aerodynamics  in  the  original  aircraft  model,  and  allows  for  a 
C  more  accurate  simulation  at  high  angles  of  attack.  Finally, 

C  like  the  original  model,  it  does  not  consider  the  Control 
C  Augumentation  System  (CAS) ,  the  Aileron-Rudder  Interconnect 
C  (ARI) ,  or  the  stick  dynamics. 

C 


IMPLICIT  DOUBLE  PRECISI0N(A-H,0-Z) 
DIMENSION  W( 150000 ) ,  IW(lOO) 


C  THIS  WAS  ADDED  SO  THAT  AUll  CAN  READ  FROM  UNIT 
C  3  ON  A  RESTART 


OPEN {UNIT»3, FILES ' f ort.3 ' ) 
OPEN (UNIT-4 , FILE- ' fort . 4 ' ) 
OPEN (UNIT-7,  FILE- ' fort.7 ' ) 
OPEN (UNIT-8 , FILE- ' fort . 8 ' ) 
OPEN (UNIT-9 , FILE- ' fort . 9 ' ) 


REWIND  7 
REWIND  8 
REWIND  9 
REWIND  3 
REWIND  4 
CALL  AUTO(W,IW) 
CLOSE (3) 

CLOSE (4) 

CLOSE (7) 

CLOSE (8) 

CLOSE ( 9 ) 

STOP 

END 
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oon>r)  nnonno  oo  on 


SUBROUTINE  FUNC  (NDIM,  NPAR,  U.  ICP,  PAR,  IJAC,  F,  DFDU,  DFDP) 


C  Revised  13  Aug  89  -  Moved  all  calls  to  subroutine  COEFF  to 
C  the  start  of  subroutine  FUMX.  Deleted  COMMON  block  SEIZE. 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 

CC»WON  /KS/  K1,K5,K7,K8,K9,K10,K12,X13,K14,K15,K16,K17 
COMMON  /ACDATA/  BWING, OWING, SREF, RHO,  RMASS, THRUST 
DOUBLE  PRECISION  Kl, K5, K7 , K8, K9, KIO, K12 , K13 , K14, K15 , K16, K17 
C  DOUBLE  PRECISION  K2 , K3, K4, K6, Kll,  IX, lY, IZ, IXZ 


C  Revised  13  Aug  89  to  dimension  DBLFl,  DELF2,  and  DX  to  15  (was 
8)  . 

C  Revised  5  Sep  89  -  Changed  dimension  DELFl,  DELF2,  and  DX  to 
50  (was  15) . 

DIMENSION  DFDU(NDIM, NDIM) ,DFDP(NDIM, NPAR) , DELFl (50) , 

+  DELF2(50) ,U(ND1M) .PAR(IO) ,F(NDIM) ,DX(50) 

C  WRITE(6,*)  'ENTERED  FUNC' 


INITIALIZE  SOME  CONSTANTS  THAT  ARE  PASSED  THROUGH 
THE  COMMON  BLOCK  ACDATA 

DATA  IS  FROM  MCAIR  REPORT#  A4172  AND  AFFTC-TR-75-32 
F-15A  APPROACH-TO-STALL/ STALL/ POST-STALL  EVALUATION 

BWING  -  A/C  WINGSPAN,  FT 

CWING  -  A/C  MEAN  AERODKMAMIC  CHORD,  FT 

SREF  -  A/C  WING  REFERENCE  AREA,  SQ  FT 

RHO  -  AIR  DENSITY  AT  20000  FT  ALTITUDE,  SLUG/FT''3 

RMASS  -  A/C  MASS,  SLUGS 

THRUST  -  TOTAL  A/C  THRUST,  LBS 

BWING»42 . 8 

CWING«15.94 

SREF«608. 

RHO». 0012673 

RMASS-37000./32.174 

THRUST-8300. 


DETERMINE  CONSTANTS  Kl  THROUGH  K17 ,  SOME  ARE  MADE  COMMON 

PASSED  TO  SUBROUTINE  FUNX  AND  USED  IN  THE  EQUATIONS 
OF  MOTION  THERE 
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INERTIAS  HAVE  UNITS  OF  SLUG-FT''2 
K1  HAS  UNITS  OF  1/FT 

K6,  K8,  Kll,  K14,  AND  K17  HAVE  UNITS  OF  l/FT^'I 


IX>  25480. 

IY»  166620. 

IZ»  186930. 

IXZ»  -1000. 

K1 • . 5 *RHO*SREF/RMASS 
K2*(IZ-IY)/IX 
K3=1XZ*IXZ/ (IX*IZ) 
K4*(IY-IX)/IZ 
K5=1XZ/IX 

K6«.5*RHO*BWING*SREF/IX 

K7=IXZ/IZ 

K8a .5*RHO*SREF*CWING/IY 

K9*{IZ-IX)/IY 

K10»IXZ/IY 

Kll* . 5*RHO*SREF*BWING/IZ 

K12a(K2+K3)/(l.-K3) 

K13»(l.-K4) *KS/ (1.-K3) 

K14«K6/(1.-K3) 

K15*(K3-K4)/(1.-K3) 

K16*(1.-kK2)*K7/(1.-K3) 

K17aKll/(l.-K3) 


K1  .  3.350088890D-04 
K5  *-3,9246467810-02 
K7  =-5.3495961050-03 
K8  *  3.6856509710-05 
K9  *  .96897131196 
KIO  --6.0016804710-03 
K12  -  .79747314581 
K13  «-9. 6157553410-03 
K14  *  6.4727458470-04 
K15  *-.754990553922 
K16  -  K13 

K17  K  8.8228515580-05 


FIND  THE  VALUES  OF  F(l)  THROUGH  F(NDIM)  .  SUBROUTINES 
COEFF  AND  FUNX  ARE  CALLED  ONCE. 


TEST  ALPHA  AND  BETA  FOR  OUT  OF  BOUNDS  CONDITION, 

-8  .GE.  ALPHA  .LE.  50  DEGREES  AND  -30  .GE.  BETA  .LE.  30 
DEGREES.  IF  A  VARIABLE  PASSES  OUTSIDE  A  BOUNDARY,  SET  THAT 
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VARIABLE  EQUAL  TO  THAT  LIMIT  AND  ISSUE  A  MESSAGE  TO  THE 
SCREEN  AND  CONTINUE  TO  EXECUTE. 

IF(U{1)  .GT.  50.)  WRITE(6.*)  'ALPHA  >  50  DEGREES.  CONTINUE 
TO  RUN' 

IF(U(1)  .LT.  -8.)  WRITE(6,*)  'ALPHA  <  -8  DEGREES,  CONTINUE 
TO  RUN' 

IF(U(2)  -GT.  30.)  WRITE(6,*)  'BETA  >  30  DEGREES,  CONTINUE 
0  RUN' 

IF(U(2)  .LT.  -30.)  WRITE ( 6, *) 'BETA  <  -30  DEGREES,  CONTINUE 
TO  RUN' 


WRITE  (6,*)  '  ALPHA.  BETA.  P.Q='.U(1)  .U(2)  .U(3)  ,U  (4) 

WRITE  (6.  •)  'R,  THETA,  PHI,  VELOCITY*'  ,U(5)  ,U{6),U(7),U(8) 

CALL  FUNX(NDIM,U,PAR,F) 


WRITE(6, *) 
WRI-TEdl,  *) 
WRITE(6, •) 
WRITE (11, *) 
WRITE(6, *) 
WRITEdl,  *) 
WRITE (6,*) 
WRITEdl,  *) 


'F1,F2=  ',Fd),F(2) 
'F1,F2=  '.Fd),F(2) 
'F3,F4»  ',F(3).F(4) 
'F3,F4=  ',F(3),F(4) 
'F5.F6*  ',F(5),F(6) 
'F5,F6=  ',F(5).F(6) 
'F7,F8»  ',F(7),F(8) 
'F7,F8x  ',F(7),F(8) 


IF(IJAC.EQ.O)  RETURN 
C  SET  THE  VALUES  OF  DX 

C  MODIFIED  TO  SCALE  DX  ACCORDING  TO  VARIABLE 
C  13  JUN  88 


C  Revised  13  Aug  89  -  Added  DX(9)  to  DXdl)  for  the  control 
C  surface  states.  Used  DX0*50.  since  this  was  already  being 
C  used  for  DPAR  when  these  were  parameters  rather  than  states. 
C  Revised  5  Sep  89  -  Added  DXd2)  for  differential  tail. 
DX0=1.0D-9 
DXd)  =DX0*50. 

DX(2) =DX0*10 . 

DX(3)*DX0*.5 

DX(4)=DX0*,25 

DX{5)=DX0*.5 

DX(6)sDX0*50. 

DX(7)*DX0*50. 

DX(8) *DX0» .5 
DX(9) *DX0*50. 

DXd0)«DX0*50. 

DXdl) -0X0*50. 

DXd2)»DX0*50. 


C  NEXT  THE  PARTIAL  OF  F  W.R.T.  A  GIVEN  PARAMETER  ARE  FINITE 
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DIFFERENCED 


PTEMP=PAR(ICP) 

PAR { ICP) -PTEMP+DX { 1 ) 

CALL  FUNX(NDIM.U.PAR.DELFl) 

PAR ( ICP) «PTEMP-DX { 1 ) 

CALL  FUNX(MDIM.U,PAR.DELF2) 

DO  13  I=1.NDIM 

DFDP(I.ICP)=(DELF1(I)-DELF2(I) )/{2.*DX(l) ) 

13  CONTINUE 

PAR(ICP)«PTEMP 

THE  NEXT  DO  LOOP  CALCULATES  THE  PARTIAL  DERIVATIVE  OF  F 
W.R.T.TO  U  USING  FINITE  DIFFERENCES. 

SET  U(J)  EQUAL  TO  U+DU,  THEN  CALL  COEFF  WITH  THIS  UPDATED 
STATE  VECTOR.  THIS  IS  DONE  SIMILARLY  WITH  U-DU 

DO  20  J=1,NDIM 

UTEMP*U(J) 

U(J)  =UTEMP4-DX(J) 

CALL  FUNX(NDIM,U,PAR,DELF1) 

U(J)aUTEMP-DX(J) 

CALL  FUNX(NDIM,U,PAR,DELF2) 

DO  16  Ixl,NDIM 

DFDU(I,J)x(DBLFl(I)-DELF2(I) )/(2.*DX(J) ) 

16  CONTINUE 

0(J) -UTEMP 

20  CONTINUE 


WRITE(6,30)  (  (DFDUd.J) ,  J*!, NDIM)  ,  1  =  1 , NDIM) 
30  F0RMAT(8{1X,E8.2)/) 

WRITE (6,*)  'LEAVING  FUNC' 

RETURN 

END 

SUBROUTINE  FUNX {NDIM,U, PAR,  F) 


SUBROUTINE  FUNX  EVALUATES  THE  NDIM  EQUATIONS  GIVEN  THE 
STATE  VECTOR  U. 

Revised  13  Aug  89  -  Added  PAR  to  the  argument  list 

for  FUNX.  Moved  the  call  for  subroutine  COEFF  from  several 
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C  places  in  FUNC  to  the  start  of  FUNX.  Added  PAR  to  DIMENSION 
C  statement . 

C  NDIM~  THE  DIMENSION  OF  THE  PROBLEM 

C  U  -  THE  VECTOR  OF  STATES  ALPHA.  BETA.  ...  (INPUT) 

C  PAR  -  THE  VECTOR  OF  PARAMETERS 

C  F  -  THE  VECTOR  RESULT  OF  FUNCTION  EVALUATIONS  (OUTPUT) 

IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

CO;4MON  /SEIZE/  CX.CY.CZ.CLM.CMM.CNM 

COMMON  /KS/  K1.K5.K7.K8.K9.K10,K12.K13.K14,K15.K16.K17 
DOUBLE  PRECISION  K1.K5.K7.K8.K9.K10.K12.K13.K14.K15.K16.K17 
DIMENSION  U(NDIM) .F(NDIM) .PAR(IO) 

C  CALL  SUBROUTINE  COEFF  TO  SET  VALUES. 

CALL  COEFF (U.NDIM) 

C  SET  TRIGONOMETRIC  RELATIONSHIPS  OF  THE  STATES  ALPHA,  BETA. 

C  THETA,  AND  PHI  AND  THEN  SET  P,  Q,  R,  AND  VTRFPS 

C  WRITE (6.*)  'ENTERED  FUNX' 

IWRITE=1 

C 

C  IF(IWRITE.EQ.1)WRITE(6. *) ' K1 , 5, 7 , 8 , 9 , 10 , 12 -17* ' , K1 , K5, K7 , K8, 

C  K9,K10.K12,K13,K14.K15,K16,K17 

C  IF(IWRITE.EQ.1)WRITE(6,*) 'FUl«  U ( 1 ) -U ( 8)  *  ' ,  U  ( 1 )  , U (2 )  , U (3 )  , 

C  U(4)  .U(5)  ,U(6)  ,U(7)  ,U(8) 

DEGRAD*57. 29577951 

CA*COS(U(l) /DEGRAD) 

SA=SIN(U(1) /DEGRAD) 

CB=COS(U(2) /DEGRAD) 

SB»SIN(U(2) /DEGRAD) 

C'rHE*COS(U(S)  /DEGRAD) 

S-THEsSIN  (U  ( 6 )  /DEGRAD) 

CPHI=COS(U (7) /DEGRAD) 

SPHI=SIN(U(7) /DEGRAD) 

P«U(3) 

Q»U(4) 

R=U(5) 

VTRFPS=1000.*U(8) 

C  Revised  13  Aug  89  -  Added  3  new  states  for  the  control 
surfaces . 

C  Deflections  are  in  degrees.  Also  added  variables  for 
commanded 

C  control  surface  deflections  in  degrees. 

C  Revised  5  Sep  89  -  Added  new  state  for  differential  tail 
control 
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C  surface  deflection. 

DSTBD«U(9) 

DRUDD»U(10) 

DAILD-U(ll) 

DDTD»U(12) 


CDSTBDsPAR(l) 

CDRUDDsPAR|2) 

CDAILDsPARO) 

C  SET  THE  GRAVITATIONAL  CONSTANT,  FT/SEC 

0=32.174 

C  THE  FOLLOWING  SYSTEM  OF  NONLINEAR  DIFFERENTIAL  EQUATIONS 
C  GOVERN  AIRCRAFT  MOTION 

C  UPDATED  FOR  PROPER  DEGREE-RADIAN  UNITS  AND  PROPERLY 
C  SCALED  VELOCITY  EQUATION:  7  JUN  88 

C  F(l) =ALPHA-DOT 

1  F  (1)  *Q-*-  ( -  (K1*VTRFPS*CX-G*STHE/VTRFPS+R*SB)  *SA+  {K1*VTRFPS 
+  *CZ+ (G*CTHE*CPHI/VTRFPS) ~P*SB) *CA) /CB 

F(1)*F(1)*DEGRAD 


C  F(2)=BETA-D0T 

2  F(2)  *-  {  (K1*VTRFPS*CX-G*STHE/VTRFPS)  ‘SBi-R)  *CA-k  (Kl*VTRFPS*CY 
-*•  ♦G*CTHE*SPHI/VTRFPS)  *CB-  ( (Kl*VTRFPS*CZfG*CTHE*CPHI/VTRFPS) 

*SB-P)*SA 
F(2) =F(2) *OEGRAD 

C 

WRITE(6, *)  'K1,CX,CY,CZ,F(2) ,C, SB, CB, CA, SA,  SPHI ,CPHI, CTHE, STHE ' , 

C  + 

K1,CX, CY,CZ,F (2 ) ,G,SB,CB,CA,SA, SPHI, CPHI, CTHE, STHE 
C  F(2)=P-D0T 

3  F(3)  =-K12*Q*R+Kl3*P*Q+K14*  (CLM-^K7*CNM)  *VTRFPS*VTRFPS 


C  F(4)=Q-D0T 

4  F(4)  =K8*VTRFPS*VTRFPS*CMM+K9*P*Rt-K10*  (R*R-P*P) 


C  F(5)=R-D0T 
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5  F(5)»K15*P*Q-K16*Q*RfK17*VTRFPS*VTRFPS*(K5*CLM+CNM) 


C  F(6)*THETA-D0T 

6  F(6)«Q*CPHI-R*SPHI 
F(6)«.F(6)*DEGRAD 


C  F (7). PHI-DOT 

7  F  (7 )  «P+Q*  (STHE/CTHE)  *SPHIi-R*  (STHE/CTHE)  *CPHI 
F(7)nF{7) *DEGRAD 


C  F(8)»VTRFP3-DOT  (SCALED  BY  A  FACTOR  OF  1000) 

8  F(8)  =U(8)  *  (  (K1*VTRFPS*CX-G«STHE/VTRFPS)  *CA*CB+  {K1*VTRFPS*CY 

♦G*CTHE*SPHI/VTRFPS)  ‘SB-*- (K1  *”rRFPS*CZ+G*CTHE*CPHI/VTRFPS)  *SA*CB) 


C  Revised  13  Aug  89  -  Added  differential  equations  governing 
control  surfaces.  Equations  assume  CAS  off  and 
Aileron-Rudder  Interconnect  off.  Deflections  are  in  degrees. 

C  F(9)»DSTBD-DOT 

9  F(9)=20.*(CDSTBD-DSTBD) 


C  F(10) =DRUDD-OOT 

10  F(10)s28.*(CDRUDD-DRUDD) 


C  F(ll) =DAILD-DOT 

11  F(ll)e20.*(CDAILD-DAILD) 


C  F;12) =DDTD-DOT 

C  Revised  5  Sep  89  -  Added  new  state  for  differential  tail 
C  deflection.  Neglecting  CAS  inputs,  the  differential  tail 
C  deflection  is  6/20  of  the  commanded  aileron  deflection, 

C  acting  through  the  stabilator  actuators. 

12  F(12)«20.*( .3*CDAILD-DDTD) 
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WRITr.(6,*)  'LEAVING  FUNX' 

WRITE(6,*)'lv  funx 
=  '  ,F(1)  ,F{2)  .F(3)  ,F{4)  .F(5).F(6)  .F(7)  ,F(8) 

RETURN 

END 

SUBROUTINE  STPNT{NDIM, U, NPAR, ICP, PAR) 


THIS  SUBROUTINE  SETS  THE  VALUES  OF  THE  STATES  AND 
PARAMETERS  AT  THE  START  OF  THE  ANALYSIS.  THE  STATES  AND 
CONTROL  SURFACE  SETTINGS  REPRESENT  AN  EQUILIBRIUM  STATE  OF 
THE  AIRCRAFT 

Revised  13  Aug  89  -  Changed  to  incorporate  control  surface 
deflections  as  states  and  commanded  control  surface 
deflections  as  parameters.  Changed  READ  statements  to  DO 
loop  format . 

Revised  5  Sep  89  -  Added  new  state  for  differential  tail 
control  surface  deflection. 

IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

DIMENSION  U(NDIM) ,PAR(10) 

WRITE ( 6 , * )  ' ENTERED  STPNT' 

U(l)  -  ALPHA.  DEG 
U(2)  -  BETA,  DEG 
U(3)  -  P,  RAD/SEC 
U{4)  -  Q,  RAD/SEC 
U{5)  -  R,  RAD/SEC 
U(6)  -  THETA,  DEG 
U{7)  -  PHI,  DEG 

U(8)  -  TRUE  VELOCITY,  IN  THOUSANDS  OF  FT/ SEC 

U(9)  -  DSTBD,  STABILATOR  DEFLECTION,  DBG 

U(10)  -  DRUDD,  RUDDER  DEFLECTION,  DEG 

U(ll)  -  DAILD,  AILERON  DEFLECTION,  DEG 

Ua2)  -  DDTD,  DIFFERENTIAL  TAIL  DEFLECTION,  DEG 

THE  STARTING  POINT  (VECTOR)  IS  READ  IN  THIS  WAY 
SO  THAT  THE  DATA  FILE  (UNIT  15)  HAS  A  COLUMN  OF 
NUMBERS.  I  BELIEVE  THAT  USING  A  COLUMN  WILL  MAKE  IT 
EASIER  TO  GET  THE  NUMBERS  RIGHT  THE  FIRST  TIME. 

OPEN (UNIT-15,  FILE- 'fort. 15 ' ) 

REWIND  (15) 

DO  10  I»1,NDIM 
READ(15,*)  U(I) 

10  CONTINUE 

U(8)»U(8)/1000. 
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PAR(1)*CDSTBD 

PAR(2)«CDRUDD  THE  PARAMETERS,  IN  DEGREES 
PAR(3)-CDAZLD 


DO  20  Zsl,3 

READdS,*)  PAR(I) 
20  CONTINUE 


RETURN 

END 


SUBROUTINE  INIT 
C  . 


IMPLICIT  DOUBLE  PRECISION (A-H,  0-2) 


COMMON  /BLCSS/ 
COMMON  /BLCPS/ 
COMMON  /BLDLS/ 
COMMON  /BLLIM/ 
COMMON  /BLOPT/ 
COMMON  /BLEPS/ 


NDIM,  ITMX,NPAR.  ICP,  IID,NMX,  IPS,  IRS 
NTST,NCOL,  IANCH,NMXPS,  lAD,  NPR,  NWTN,  ISP,  ISWl 
DS, DSMIN, DSMAX, IADS 
RLO.RL1,AO,A1,PAR(10) 

ITNW,MXBF, IPLT, ICP2, ILP 
EPSU,EPSL,EPSS,EPSR 


IN  THIS  SUBROUTINE  THE  USER  SHOULD  SET  THOSE  CONSTANTS  THAT 
REQUIRE  VALUES  DIFFERENT  FROM  THE  DEFAULT  VALUES  ASSIGNED  IN 
THE  LIBRARY  SUBROUTINE  DFINIT.  FOR  A  DESCRIPTION  OF  THESE 
CONSTANTS  SEE  THE  DOCUMENTATION  CONTAINED  IN  THJ  LIBRARY. 

COMMON  BLOCKS  CORRESPONDING  TO  CONSTANTS  THAT  THE  USER  V/ANTS  TO 
CHANGE  MUST  BE  INSERTED  ABOVE.  THESE  COMMON  BLOCKS  SHOULD  OF 


C 

C 


COURSE  BE 

IDENTICAL 

DSMAX 

U 

10. 

DSMIN 

M 

.0000001 

EPSU 

S 

l.OD-07 

EPSL 

C 

l.OD-07 

EPSS 

9 

1.00-05 

EPSR 

9 

l.OD-07 

lAD 

U 

1 

ILP 

8 

1 

ITMX 

8 

10 

ITNW 

8 

5 

MXBF 

8 

5 

Revised 

13 

Aug  89 

Revised 

5 

Sep  89  - 

NDIM 

a 

12 

NPAR 

a 

3 

-  Changed  NDIM  from  8  to  11. 
Changed  NDIM  from  11  to  12, 
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OPEN(UNITa25,FILE«'fort.25' ) 

REWIND  (25) 

READ(25,*)  RLO.RLl 
READ(25.*)  AO.Al 
READ(25.*)  DS 
R£AD(25,*)  NMX 

READ(25,*)  NTST,NCOL,NMXPS,NPR 
READ(25, •)  ISP. IRS. ICP. ICP2. IPLT.  IPS 
READ{25.*)  ISWl 


C  WRITE (6.*)  'LEAVING  INIT' 

RETURN 

END 

SUBROUTINE  BOND 


ftETURN 

END 

SUBROUTINE  ICND 


RETURN 

END 


SUBROUTINE  COEFF (U.NDIM) 
Last  Edited:  30  May  90  by  J.Beck 


C 


C  This  subroutine  computes  the  aerodynamic  force  and  moment 
C  coefficients  to  be  used  in  subroutine  FUNX.  It  was  taken 
C  directly  from  Dan  Baumann's  program  D2ICC.  The  only  changes 
C  made  were  to  incorporate  the  CONSTANT  common  block,  to  change 
C  the  control  surface  deflections  from  parameters  to  states 
C  (i.e.,  DSTBDbU(13)  rather  than  D3TBD>PAR(1) ) ,  to  remove  PAR 
C  from  the  argument  list,  to  use  the  1988  data  base  control 
C  surface  deflection  names  (DSTBD,  DAILD,  DTALD,  and 
C  DRUOD),  and  to  use  the  differential  tail  deflection  state 
C  rather  than  .3*DAILD. 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-2) 

DIMENSION  U(NDIM) 

COMMON  /ACDATA/  BWING.CWING, SREF,  RHO.  RMASS, THRUST 
COMMON  /SEIZE/  CX.CV.CZ.CLM.CMM.CNM 


90 


C  THE  PRIMARY  SOURCE  OF  THESE  COEFFICIENT  EQUATIONS  IS 
C  SUBROUTINE  AROlO  FROM  MCAIR  CODE  USED  IN  THE  FI 5  BASELINE 

C  SIMULATOR.  MOST  OF  THE  COEFFICIENTS  USED  IN  THE  EQUATIONS 

C  WERE  COMPUTED  USING  SAS  WITH  RAW  DATA  FROM  THE  F15  SIMULATOR 

C  DATA  TABLES.  THIS  SUBnOUTINE  IS  CALLED  BY  THE  DRIVER  PROGRAM 

C  FOR  THE  AUTO  SOFTWARE.  IT  MERELY  TAKES  INPUTS  ON  THE  A/C 
C  STATE  AND  CONTROL  SURFACE  POSITIONS  AND  RETURNS  THE 
C  APPROPRIATE  AERO  COEFFICIENTS  CX,  CY,  CZ,  CL,  CM,  AND  CM. 

C  22  July  1990-  Due  Co  Che  face  chac  Chis  is  a  merger  of  a 

C  lacer  subroucine  Co  an  earlier  version,  some  of  Che 

C  complexicy  of  Che  lacer  version  has  co  be  caken  ouC .  In 
C  chis  case  ic  involves  rescacing  cercian  conscancs. 

C  GRAVITATIONAL  CONSTANT 

G=32.174 

C  DEGREE  TO  RADIAN  CONVERSION 

DEGRAD=57. 2957795131 


C  PI 

PI=3. 1415926536 


C  INPUTS  TO  THIS  SUBROUTINE 


C  AL 
C  BETA 
C  DDA 
C  DELEDD 
C  DELESD 
C  DRUDD 


ANGLE  OF  ATTACK,  DEG 

SIDESLIP  ANGLE,  DEG 

AILERON  DEFLECTION  ANGLE,  DEG 

DIFFERENTIAL  TAIL  DEFLECTION  ANGLE,  DEG 

SYMMETRICAL  TAIL  DEFLECTION  ANGLE,  DEG 

RUDDER  DEFLECTION,  POSITIVE  TRAILING  EDGE  LEFT,  DEG 


C  P 
C  Q 
C  R 
C  THRUST 
C  VTRFPS 


ROLL  RATE,  RAD/SEC 

PITCH  RATE,  RAD/SEC 

YAW  RATE,  RAD/ SEC 

TOTAL  ENGINE  THRUST,  POUNDS 

TRUE  AIRSPEED,  FT/SEC 


C  INTERMEDIATE  VARIABLES  USED  IN  THIS  SUBROUTINE 


C  ABET 
C  ARUD 
C  BWING 
C  CA 
C  CD 
C  CL 
C  ewING 
C  DAHD 
C  DAHLD 
C  DAHRD 
C  DELEDR 
C  DELESR 


ABSOLUTE  VALUE  OF  BETA,  DEG 

ABSOLUTE  VALUE  OF  RUDDER  DEFLECTION,  DEG 

WING  SPAN,  FEET 

COSINE  RAL  (RAL  IN  RADIANS) 

COEFFICIENT  OF  DRAG 
BASIC  LIFT  COEFFICIENT 
MEAN  AERODYNAMIC  CHORD,  FEET 
DIFFERENTIAL  ELEVATOR  DEFLECTION,  DEG 
LEFT  AILERON  DEFLECTION,  DEG 
RIGHT  AILERON  DELFECTION,  DEG 
DIFFERENTIAL  TAIL  DEFLECTION  ANGLE,  RAD 
SYMMETRIC  TAIL  DEFLECTION  ANGLE,  RAD 
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C  QBARS  -  DYNAMIC  PRESSURE  TIMES  WING  REFERENCE  AREA,  LBF 
C  RABET  -  ABSOLUTE  VALUE  OF  BETA,  RADIANS 

C  RAL  -  ABSOLUTE  VALUE  OF  ALPHA,  RADIANS 

C  RARUD  -  ABSOLUTE  VALUE  OF  RUDDER,  RADIANS 

C  SA  -  SINE  RAL  (RAL  IN  RADIANS) 

C  OUTPUTS  FROM  THIS  SUBROUTINE 

C  CX  *  BASIC  AXIAL  FORCE  COEFFICIENT,  BODY  AXIS,  f  FORWARD 

C  CY  -  BASIC  SIDE  FORCE  COEFFICIENT,  BODY  AXIS,  ♦  RIGHT 

C  CZ  -  BASIC  NORMAL  FORCE  COEFFICIENT,  BODY  AXIS,  +  DOWN 

C  CLM  -  BASIC  ROLLING  MOMENT  COEFFICIENT,  BODY  AXIS,  +  R  WING 
C  DOWN 

C  CMM  -  BASIC  PITCHING  MOMENT  COEFFICIENT,  BODY  AXIS,  NOSE  UP 

C  CNM  -  BASIC  YAWING  MOMENT  COEFFICIENT,  BODY  AXIS,  ♦  NOSE 

C  RIGHT 

C  ANGLES  USED  IN  CALCULATING  CL,  CLLDB,  ...,  ARE  IN  RADIANS. 

C  THIS  IS  BECAUSE  RADIANS  WERE  USED  IN  THE  CURVE  FITTING 
C  PROGRAM  TO  OBTAIN  THE  COEFFICIENTS  OF  THE  ALPHA,  BETA,  .  .  . , 
C  TERMS  IN  THE  FOLLOWING  EQUATIONS. 

C  MOMENT  REFERENCE  CENTER  WAS  SET  IN  AROlO  PROGRAM  AS: 

C  DATA  CMCGR  /.2565/,  CNCGR  /.2565/ 

C  THE  AERO  STABILITY  DATA  WAS  TAKEN  REFERENCED  TO  THESE  CG 
C  LOCATIONS.  THE  MOMENTS  OF  INERTIA  AND  OTHER  AIRCRAFT  DATA 

C  ARE  FOR  A  CLEAN  CONFIGURATION  TEST  AIRCRAFT  WITH  A  CG  AT 

C  TKF  SAME  CG.  AS  A  RESULT,  THERE  IS  NO  'CG  OFFSET'  TO  BE 
C  COMPUTED. 

AL»U(1) 

BETA»U(2) 

P«U(3) 

Q-U(4) 

R-U(5) 

THETA»U ( 6 ) 

PHI=U(7) 

VTRFPS»U<8)*1000. 

DSTBD-U ( 9 ) 

DAILD-U(IO) 

DTALD»U(11) 

DRUDD>U(12) 

QBARS- . 5 'RHO'VTRFPS’VTRFPS'SREF 
C02V-CWING/ ( 2 . *VTRFPS) 

B02V-BWING/ (2 . *VTRFPS) 

QSB«BWING*QBARS 
ARUDsABS (DRUDD) 

RARUD-ARUD/ DEGRAD 
RAL-AL/DEGRAD 
ABET- ABS (BETA) 
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RABETsAfiET/DEGRAD 
RBETAbBETA/DEGRAD 
DAILA-ABS (DAILD) 
DSTBR-DSTBD/ DEGRAD 


C 


NEW  SECTION  OF  CODE  -  1)  ALL  THE  AERODYNAMIC  COEFFICIENTS 

IN  THIS  VERSION  OF  THE  DRIVER 
PROGRAM  ARE  TAKEN  DIRECTLY  FROM 
THE  1988  F15  AEROBASE  (0.6  MACH, 

20000  FEET) 


2)  THIS  SECTION  SUMMARIZES  THE 
AERODYNAMIC  COEFFICIENTS  AS  TO 

WHAT  THEY  ARE  AND  HOW  THEY  ARE 
USED.  THE  FIRST  ACCRONYM  IS  THE 
JOVIAL  NAME  OF  THE  AERODYNAMIC 
COEFFICIENT  (CFXl,  ETC),  THE 
SECOND  ACCRONYM  IS  THE  FI 5 
,SE  CODE  OR  CTAB  NAME 

(ATAB15,  ETC).  A  BRIEF  DEFINITION 
OF  THE  AERODYNAMIC  COEFFICIENT  IS 
ALSO  PROVIDED. 

3)  THERE  IS  ALSO  A  SECTION  THAT 
PROVIDES  A  TABLE  OF  CONVERSIONS 

BETWEEN  WHAT  THE  VARIABLE  IS 
CALLED  IN  THE  ORIGINAL  SECTION  OF 
THIS  PROGRAM  AND  ITS  NAME  IN  THE 
1988  FI 5  AEROBASE. 

FOR  THE  SAKE  OF  CONTINUITY  THE 
ORIGINAL  PROGRAM  NAME  IS  USED  AND 
THE  1988  F15  AEROBASE  NAME 
IS  PROVIDED  AS  BOOK  KEEPING 
INFORMATION. 


CFX  -  FORCE  IN  STABILITY  AXIS  X  DIRECTION  (CD  IN  BODY  AXIS) 
(FUNCTION  OF  CL  OR  CFZl) 

CFX  •  CFXl  +  CXRB  ♦  STORE  INCREMENTS  ♦  CXDSPD  +  DCXLG  +  DCD 

CFXl  ■  ATAB15  »  PERFORMANCE  DRAG  COEFFICIENT  -  CD 
CXRB  *  ATAB22  -  DELTA  CD  DUE  TO  CG  (»0.0) 

CXDSPD  *  ATAB27  «  DELTA  CD  DUE  TO  SPEEDBRAKE  (NORMALLY  » 
0.0436) 

C  SET  TO  0  SINCE  THIS  STUDY  IS  CONCERNED 
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WITH  HIGH  ANGLES 

OF  ATTACK  PHENOMENON  (>40  DEGREES) 
AND  BECAUSE  THE  SPEEDBRAKE  WILL  NOT  DEPLOY 
AT  ANGLES  OF  ATTACK  GREATER  THAN  15  DEGREES. 

DCXLG  =  ATAB19  =  DELTA  CD  DUE  TO  REYNOLD'S  NUMBER  (=-0.0005) 
DCD  =  BTAB03  =  DELTA  CD  DUE  TO  2-PLACE  CANOPY  (F15B)  (=0.0005) 

*******  NOTE  THAT  DCXLG  AND  DCD  CANCEL  EACH  OTHER 


C 


CFY  =  FORCE  IN  BODY  AXIS  Y  DIRECTION 

CFY  =  CFY1*EPA02  ♦  CYDAD*DAILD  t-  (CYDRD*DRUDD*DRFLX5 )  *EPA43 
t- (CYDTD*DTFLX5  -*•  DTFLX6 )  *DTALD  +  CFYP'PB  t  CFYR*RB 
■►CYRB  *  STORE  INCREMENTS  +  DCYB*BETA 

CFYl  =  ATAB16  *  BASIC  SIDE  FORCE  COEFFICIENT  -  CY(BETA) 

EPA02  =  ATAB21  =  BETA  MULTIPLIER  TABLE 
CYDAD  =  ATAB75  =  SIDE  FORCE  COEFFICIENT  DUE  TO  AILERON 
DEFLECTION 

DAILD  =  AILERON  DEFLECTION  (DEG) 

CYDRD  =  ATAB69  =  SIDE  FORCE  COEFFICIENT  DUE  TO  RUDDER 
DEFLECTION 

DRUDD  =  RUDDER  DEFLECTION  (DEG) 

DRFLX5  =  ATAB88  =  FLEX  MULTIPLIER  ON  CYDRD  (=0.89) 

EPA43  =  ATAB30  =  MULTIPLIER  ON  CNDR,  CLDR,  CYDR  DUB  TO 
SPEEDBRAKE 
(=1.0) 

CYDTD  =  ATAB72  »  SIDE  FORCE  COEFFICIENT  DUE  TO  DIFFERENTIAL  C 
TAIL 

DEFLECTION  -  CYDDT 

DTFLX5  =  ATABIO  =  FLEX  MULTIPLIER  ON  CYDTD  (=0.975) 

DTFLX6  »  ATAB77  =  FLEX  INCREMENT  TO  CYDTD  (=0.0) 

DTALD  =  DIFFERENTIAL  TAIL  DEFLECTION  (DEG)  WHICH  IS 

DIRECTLY  PROPORTIONAL  TO  AILERON  DEFLECTION 
AND  IS  PRIMARILY  USED  TO  ASSIST  IN  ROLLING 
THE  F-15B  (DTALD=0.3*DAILD) 

CFYP  =  ATAB13  =  SIDE  FORCE  COEFFICIENT  DUE  TO  ROLL  RATE  (CYP) 
PB  ■  (PEOBB*SPAN)/ (2*VILWF) 

PEOBB  =  ROLL  RATE  IN  RAD/ SEC  =  P 
SPAN  »  WING  SPAN  =42.8  FEET  =  BWING 
VILWF  *  VELOCITY  IN  FT/SEC  =  VTRFPS 
CFYR  =  ATAB07  =  SIDE  FORCE  COEFFICIENT  DUE  TO  YAW  RATE  (CYR) 

RB  *  (REOBB*SPAN)/(2*VILWF) 

REOBB  «  YAW  RATE  IN  RAD/ SEC  =  R 
CYRB  «  ATAB93  =  ASSYMETRIC  CY  AT  HIGH  ALPHA  (ANGLE  OF  ATTACK) 

DCYB  =0.0  THERE  IS  NO  INCREMENT  DELTA  CYB  (SIDE  FORCE) 

DUE  TO  A  2 -PLACE  CANOPY  ON  THE  F15B.  THIS  IS 
BECAUSE  THE  SAME  CANOPY  IS  USED  ON  BOTH  THE 
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BASELINE  F15A  AND  THE  FISB.  THE  SIDEFORCE  IS  THE 
SAME  FOR  BOTH  VERSIONS  OF  THE  F15  AND  ALREADY 
INCLUDED  IN  THE  BASIC  SIDE  FORCE  (CFYl) .  THE  TWO 
PLACE  CANOPY  IS  MOUNTED  DIFFERENTLY  HOWEVER,  SO 
THERE  IS  A  DIFFERENCE  IN  YAWING  AND  ROLLING  MOMENT. 
(SEE  DCNB  AND  DCLB) 


CFZ  =  FORCE  IN  STABILITY  AXIS  2  DIRECTION  (CL  IN  BODY  AXIS) 
CF2  =  CFZl  +  CZDSPD  +  STORE  INCREMENTS  +  DCL*BETA 


CFZl  =  ATAB17  =  BASIC  LIFT  COEFFICIENT  -  CL 
CZDSPD  =  ATAB26  =  DELTA  CL  DUE  TO  SPEEDBRAKE 

SET  TO  0  DUE  TO  THE  REASONS  GIVEN  ABOVE  IN  CXDSPD 
DCL  =  BTABOl  =  DELTA  CL  DUE  TO  2 -PLACE  CANOPY  (F15B)  (=0.0) 


CML  »  TOTAL  ROLLING  MOMENT  COEFFICIENT  IN  BODY  AXIS 
CML  *  CML1*EPA02  +  CLDAD*DAILD  (CLDRD*DRUDD*DRFLX1 )  *EPA43 
(CLDTD'DTFLXl  +  DTFLX2 )  *DTALD  CMLP*PB  ♦  CMLR*RB  + 
STORE  INCREMENTS  +  CLDSPD  ♦  DCLB ‘BETA 


CMLl  =  ATABOl  =  BASIC  ROLLING  MOMENT  COEFFICIENT  -  CL (BETA) 
EPA02  =  ATAB21  =  BETA  MULTIPLIER  TABLE 

CLDAD  =  ATAB73  =  ROLL  MOMENT  COEFFICIENT  DUE  TO  AILERON  C 
DEFLECTION 
- (CLDA) 

DAILD  =  AILERON  DEFLECTION  (DEG) 

CLDRD  =  ATAB67  =  ROLLING  MOMENT  COEFFICIENT  DUE  TO  RUDDER 
DEFLECTION  - (CLD) 

DRUDD  *  RUDDER  DEFLECTION  (DEG) 

DRFLXl  =  ATAB80  =  FLEX  MULTIPLIER  ON  CLDRD  (=0.85) 

EPA43  =  ATAB30  =  MULTIPLIER  ON  CNDR,  CLDR,  CYDR  DUE  TO 
SPEEDBRAKE 
(=1.0) 

CLDTD  =  ATAB70  =  ROLL  MOMENT  COEFFICIENT  DUE  TO  DIFFERENTIAL 
TAIL 

DEFLECTION  -  CLDD 

DTFLXl  =  ATAB04  =  FLEX  MULTIPLIER  ON  CLDTD  (=0.975) 

DTFLX2  =  ATAB84  =  FLEX  INCREMENT  TO  CLDTD  (=0.0) 

DTALD  =  DIFFERENTIAL  TAIL  DEFLECTION  (DEG)  WHICH  IS 
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DIRECTLY  PROPORTIONAL  TO  AILERON  DEFLECTION 
AND  IS  PRIMARILY  USED  TO  ASSIST  IN  ROLLING 
THE  F-15B  (DTALD  =  0.3*DAILD) 

CMLP  =  ATAB02  *  ROLL  DAMPING  DERIVATIVE  -CLP 

PB  =  (PEOBB*SPAN> / (2*VILWF) 

PEOBB  =  ROLL  RATE  IN  RAD/ SEC  =  P 
SPAN  =  WING  SPAN  =42.8  FEET  =  BWING 
VILWF  =  VELOCITY  IN  FT/SEC  =  VTRFPS 
CMLR  =  ATABll  =  ROLLING  MOMENT  COEFFICIENT  DUE  TO  YAW  RATE  - 

CLR 

RB  =  (REOBB*SPAN) / (2*VILWF) 

REOBB  =  YAW  RATE  IN  RAD/SEC  =  R 
CLDSPD  =  ATAB29  =  DELTA  CL  DUE  TO  SPEEDBRAKE 

SET  TO  0  DUE  TO  THE  REASONS  GIVEN  ABOVE  IN 
CXDSPD 

DCLB  =  BTAB04  =  INCREMENT  DELTA  CLB  (ROLLING  MOMENT)  DUE  TO  C 
2 -PLACE  CANOPY  FROM  PSWT  499 


C 


CMM  =  TOTAL  PITCHING  MOMENT  COEFFICIENT  IN  STABILITY  AXIS 
(BODY  AXIS  -  AS  WELL) 

CMM  =  CMMl  CMMQ*QB  +  STORE  INCREMENTS  +  CMDSPD  ♦  DCM 

CMMl  =  ATAB03  s  BASIC  PITCHING  MOMENT  COEFFICIENT  -  CM 
CMMQ  =  ATAB05  »  PITCH  DAMPING  DERIVATIVE  -  CMQ 
QB  *  (QEOBB*MAC)/(2*VILWF) 

QEOBB  =  PITCH  RATE  IN  RAD/ SEC  =  Q 
MAC  *  MEAN  AERODYNAMIC  CHORD  a  15.94  FEET  =  C 
OWING 

VILWF  =  VELCXTITY  IN  FT/SEC  =  VTRFPS 
CMDSPD  =  ATAB25  «  DELTA  CM  DUE  TO  SPEEDBRAKE 

SET  TO  0  DUE  THE  REASONS  GIVEN  ABOVE  IN  C 
CXDSPD 

DCM  =  BTAB02  =  DELTA  CM  DUE  TO  2 -PLACE  CANOPY  (F15B)  (=0.0) 


’.*C 


CMN  =  TOTAL  YAWING  MOMENT  COEFFICIENT  IN  BODY  AXIS 

CMN  =  CMN1*EPA02  +  CNDAD'DAILD  ♦  [CNDRD*DRUDD*DRFLX3 ] *EPA43 

+[CNDTD*DTLX3  +  DTFLX4 ] ‘DTALD  +  CMNP*PB  +  CMNR*RB  +  CNRB 
+DCNB2*EPA36  +  STORE  INCREMENTS  +  CNDSPD  DCNB*BETA 


CMNl  =  ATAB12  =  BASIC  YAWING  MOMENT  COEFFICIENT  -  CN  (BETA) 

EPA02  =  ATAB21  =  BETA  MULTIPLIER  TABLE 

CNDAD  =  ATAB74  =  YAW  MOMENT  COEFFICIENT  DUE  TO  AILERON 

DEFLECTION  -CNDA 
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C  DAILD  =  =  AILERON  DEFLECTION  (DEG) 

C  CNDRD  =  ATAB68  =  YAWING  MOMENT  COEFFICIENT  DUE  TO  RUDDER 
C  DEFLECTION  -CNDR 

C  DRUDD  *  RUDDER  DEFLECTION  (DEG) 

C  DRFLX3  *  ATAB85  *  FLEX  MULTIPLIER  ON  CNDRD 
C  EPA43  =  ATAB30  *  MULTIPLIER  ON  CNDR,  CLDR,  CYDR  DUE  TO  C 
SPEEDBRAKE 

CNDTD  =  ATAB71  a  YAWING  MOMENT  COEFFICIENT  DUE  TO  DIFFERENTIAL 
TAIL 

DEFLECTION  -  CNDDT 

DTFLX3  a  ATAB08  «  FLEX  MULTIPLIER  ON  CNDTD 
DTFLX4  a  ATAB09  a  FLEX  INCREMENT  ON  CNDTD  (=0.0) 

DTALD  a  a  DIFFERENTIAL  TAIL  DEFLECTION  (DEG)  WHICH  IS 

DIRECTLY  PROPORTIONAL  TO  AILERON  DEFLECTION 
AND  IS  PRIMARILY  USED  TO  ASSIST  IN  ROLLING 
THE  F-15B  (DTALD  =  0.3*DAILD) 

CMNP  a  ATAB06  =  YAWING  MOMENT  COEFFICIENT  DUE  TO  ROLL  RATE  - 
CNP 

PB  =  (PEOBB*SPAN) / (2*VILWF) 

PEOBBaROLL  RATE  IN  RAD/ SEC  a  P 
SPAN  a  WING  SPAN  a  42.8  FT  =  SWING 
VILWF  a  VELOCITY  IN  FT/ SEC  a  VTRFPS 
CMNR  a  ATAB14  =  YAW  DAMPING  DERIVATIVE  -  CNR 
RB  a  (REOBB*SPAN) / (2*VILWF) 

REOBB  a  YAW  RATE  IN  RAD/ SEC  =  R 
CI^RB  a  ATAB86  =  ASSYMETRIC  CN  AT  HIGH  ALPHA 

DCNB2  a  ATAB44  a  DELTA  CNB  WITH  STABILATOR  EFFECT  -  DELCNB  C 
(aO.O) 

EPA36  a  ATAB94  *  MULTIPLIER  ON  DCNB2  (aBETA) 

CNDSPD  a  ATAB28  =  DELTA  CN  DUE  TO  SPEEDBRAKE 

SET  TO  0  DUE  TO  THE  REASONS  GIVEN  ABOVE  IN 
CXDSPD 

DCNB  a  BTAB05  =  INCREMENT  DELTA  CNB  (YAWING  MOMENT)  DUE  TO 
2 -PLACE  CANOPY  (F15B) 


C 


MISCELLANEOUS  COEFFICIENTS  AND  NAME  CONVERSION  TABLE 


1988  F15 
AEROBASE  NAME 

AL77D 

BE77D 

BE77D 


ORIGINAL 
PROGRAM  NAME 

AL 

BETA 

RBETA 


DEFINITION 


ANGLE  OF  ATTACK 
(DEG) 

SIDESLIP  ANGLE 
(DEG) 

SIDESLIP  ANGLE 
(RAD) 
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B077D 

ABET 

ABSOLUTE  VALUE  OF 
SIDESLIP  ANGLE 
(DEG) 

DAI  LA 

DAILA 

ABSOLUTE  VALUE  OF 
AILERON  DEFLEC¬ 
TION  (DEG) 

DAILD 

DDA 

AILERON  DEFLEC¬ 
TION  (DEG) 

DRUABS 

ARUD 

ABSOLUTE  VALUE  OF 
RUDDER  DEFLEC¬ 
TION  (DEG) 

DRUABS 

RARUD 

ABSOLUTE  VALUE  OF 
RUDDER  DEFLEC¬ 
TION  (RAD) 

DRUDD 

DRUDD 

RUDDER  DEFLECTION 
(DEG) 

DST3D 

DSLESD(R) 

AVERAGE 

STABILATOR 

DEFLECTION 

DEG  (RAD) 

DTALD 

DELEDD(R) 

DIFFERENTIAL  TAIL 
DEFLECTION 

DEG  (RAD) 

PB= (P'BWING) / (2*VTRFPS) 

QB=(Q*CWING) / {2*VTRFPS) 

RBa {R*BWING) / (2*VTRFPS) 

THE  F-15B  AERO  DATA  TABLES  DO  NOT  CONTAIN  STABILITY 
COEFFICIENT 

DATA  FOR  BETA  AND  RUDDER  DEFLECTION  ,DRUDD,  LESS  THAN  0 

DEGREES.  THE  ABSOLUTE  VALUE  OF  BETA,  ABET,  AND  THE  ABSOLUTE 

VALUE  OF  RUDDER  DEFLECTION,  ARUDD,  ARE  USED  IN  THE  FOLLOWING 

EQUATIONS.  IN  RADIANS  THESE  PARAMETERS  ARE  RABET  AND  RARUD, 

RESPECTIVELY.  IN  SOME  CASES  THE  COEFFICIENT  IS  MULTIPLIED  BY 
A 

-1  FOR  PARAMETER  VALUES  LESS  THAN  ZERO. 

EPA02  IS  A  MULTIPLIER  THAT  ADJUSTS  THE  PARTICULAR  COEFFICIENT 

C  IT  IS  WORKING  ON  (CFY1,CML1,CMN1)  BY  CHANGING  THAT  PARTICULAR 
C  COEFFICIENTS  SIGN  (POSITIVE  OR  NEGATIVE)  DEPENDENT  ON  THE  SIGN 
C  OF  THE  SIDESLIP  ANGLE  (BETA) .  IF  BETA  IS  NEGATIVE  THEN 
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C  EPA02=-1.0.  IF  BETA  IS  POSITIVE  THEN  EPA02=1.0.  SINCE  THIS 

C  FUNCTION  IS  DISCONTINUOUS  AT  THE  ORIGIN  A  CUBIC  SPLINE  HAS 

C  BEEN  EMPLOYED  TO  REPRESENT  THIS  FUNCTION  IN  ORDER  THAT 

C  AUTO  CAN  RUN. 

C 

C  Revised  15  Jun  90  -  Combined  three  IF  statements  into 

IF-THEN-ELSE. 

IF ( BETA. LE.- 1.0) THEN 
EPA02S=-1 .00 
ELSEIF{BETA.GE.i.O)THEN 
EPA02S=1.00 
ELSE 

EPA02S  =  -1.00f  (1.5*(  (EETA+1.0)**2)  )  -  (0.5*  (  (BETA-^1 .0)  **3)  ) 
ENDIF 
C 

C  Revised  ?.6  Jun  90  -  Combined  three  IF  statements  into 

IF-THEN-ELSE. 

IF ( BETA. LE .-5.0) THEN 
EPA02L=-1 .00 
ELSEIF { BETA .GE .5.0) THEN 
EPA02Lrl .00 
ELSE 

EPA02L=-1.00-^(0 .06*  (  (BETA+5.  O'  **2)  )  -(0.004*  (  (BETA+5 . 0)  **3  J  ) 

ENDIF 

C 

C 

C 

c 

c 

c 

C*  ****** . **** . . . . . 

c 

c 

c 

CFZl=-0.00369376+(3 .78028702*FJL)+ i0.6921459*RAL*RAL) - (5.0005867 

-h*  (RAL**3)  )♦  (1 . 94478199*  {RAL**4)  )  + (0 .40781955*DSTBR)  i- (0 . 10114579 
+• (DSTBR*DSTBR) ) 

C 

CFZ*CFZ1 

C 

C 

c*  ******** . . . . . . . . . 

c 

c 

c 

c 

CL=CFZl/57 .29578 
C 
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THIS  CONVERSION  OF  CFZl  TO  CL  IS  AN  ARTIFACT  FROM  THE 
CURVE  FITTING  PROCESS  WHERE  ALL  THE  INDEPENDENT  VARIABLES 
WERE  ANGLES  THAT  WERE  CONVERTED  FROM  DEGREES  TO  RADIANS. 
IT  JUST  SO  HAPPENED  THAT  FOR  CFXl  ONE  OF  THE  VARIABLES 
WAS  NOT  AN  ANGLE  BUT  A  DIMENSIONLESS  COEFFICIENT. 


CFXl-0. 018068211- (0.01556573*CL)  + (498. 96208868*CL*CL) 

1-- (14451.56518396*  (CL**3)  )♦  (2132344 . 6184755*  (CL**4 )  ) 

TRANSITIONING  FROM  LOW  AOA  DRAG  TABLE  TO  HIGH  AOA  DRAG  TABLE 

CFX2=0. 0267297- {0.10646919*RAL)i- (5. 39836337*RAL*RAL) 
-*•-(5.0086893*  (RAL* *3 ))♦(  1 . 34148193*  (RAL**4 )  )  + 

1- (0.20978902*DSTBR)i-(0. 30604211*  (DSTBR**2)  >1-0.09833617 

Als20.0/DEGRAD 

A2s30.0/DEGRAD 

A12sAli-A2 

BA=2.0/(-Al**3i-3  .*A1*A2*(A1-A2)i-A2**3) 

BB«-.3  ,*BA*(Ali-A2)/2. 

BCs3.*BA*A1*A2 

BDsBA*A2  *  *2  *  ( A2  -3  .  'AD  /2  . 

F1sBA*RAL**3i-BB*RAL**2i-BC*RALi-BD 
F2»-BA*RAL**3i-(3  .*A12*BAi-BB)  *RAL**2- 
1-  (BCi-2  .  •A12*BBi-3  .*A12**2*BA)  *RAL* 

♦  BDi-A12*BCi-A12**2*BB-*A12**3*SA 

IF  (RAL  .LT.  Al)  THEN 

CFXsCFXl 

ELSEIF  (RAL  .GT.  A2)  THEN 
CFX»CFX2 
ELSE 

CFX=CFXl*Fli-CFX2*F2 

ENDIF 


C 


C 

DTFLX5«0.975 
DRFLXSsO .89 
C 

CFY1  =  -W.  05060386- (0.12342073*RAL)i-(  1. 04501 136*RAL*RAL) 
1-- (0.17239516*  (RAL*  *3)  >-(2.90979277*  {RAL**4)  ) 

1-1- (3. 06782935*  (RAL*  •5>  )- (0 . 88422116*  (RAL**6 )  ) 
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♦- ,0. 06578812 *RAL*RABET) - (0 .71521788*RABET) - ( 0 . 00000475273 
(i^ABET**2 )  )  -  ( 0 , 04856168*RAL*DSTBR)  -  ( 0 . 05943607 *RABET»DSTBR)  ♦ 
+(0.02018534*DSTBR) 

C 

IF  (RAL  .LT.  .52359998)  THEN 
C 

CFYP.O. 014606188^ (2. 52405055*RAL)- (5. 02687473* (RAL**2) ) 

♦- (106. 43222962* (PAL**3) )♦ (256.80215423* (RAL**4) ) 

♦•►(  1256. 39636248*  (RAL**5) ) 

V- (3837 .92878173* (RAL**6) )- (2863 .16083460* (RAL**7) ) + 

+ (17382. 72226362* (RAL**8) )- (13731. 65408408* (RAL**9) ) 

ENDIF 

C 

IF  ((RAL  .GE.  .52359998)  .AND.  (RAL  .LE.  .610865))  THEN 
C 

CFyPa0.00236511f (0.52044678* (RAL-0. 5235 9998)  )- (12 . 8597002* (RAL- 
+0.52359998) **2) + (75.46138* (RAL-0 . 52359998) * *3) 

ENDIF 

C 

IF  (RAL  .GT.  0.610865)  THEN 
C 

CFYPsO.O 

ENDIF 

C 

IF  (RAL  .LT.  -0.06981)  THEN 
C 

CFYR=0.35 

ENDIF 

C 

IF  ((RAL  .GE.  -0.06981)  .AND.  (FAL  .LT.  0.0))  THEN 
C 

CFYRsO. 34999999+ (35. 4012-.13*(RAL+0. 06981)* *2) -(493  .33441162* 
+(RAL+0. 06981) **3) 

ENDIF 

C 

IF  ((RAL  .GE.  0.0)  .AND.  (RAL  .LE.  0.523599))  THEN 
C 

CFYRsO . 35468605- (2 .26998141 *RAL) ♦ (51 . 82178387 *RAL*RAL) 
+-(718.55069823*(RAL**3) ) 

++ (4570. 00492172* (RAL**4) )- (14471. 88028351* (RAL*‘5) )+ 

+ (22026. 58930662* (RAL**6) )- (12795. 99029404* (RAL* *7) ) 

ENDIF 

C 

IF  ((RAL  .GT.  0.523599)  .AND.  (RAL  .LE.  0.61087))  THEN 
C 

CFyRs0.00193787+( 1.78332496*(RAL-0. 52359903)  )  -(41.63198853* (RAL- 
+0.52359903) * *2) +(239. 97909546* (RAL-0 . 52359903 ) **3 ) 

ENDIF 

C 
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IF  (RAL  .GT.  0.61087)  THEN 
C 

CFYR=0 . 0 

ENDIF 

C 

IF  (RAL  .LT.  0.55851)  THEN 
C 

CYDAD--0 .00020812+ (0 .00062122*RAL) ♦ (0 .00260729*RAL*RAL) 

++ (0. 00745739* (RAL**3) )- (0.0365611* (RAL**4) ) 
+-(0.04532683*(RAL**5) )+(0.20674845*(RAL**6) ) 

T- (0. 13264434* (RAL**7) )- (0.00193383* (RAL**8) ) 

ENDIF 

C 

IF  ((RAL  .GE.  0.55851)  .AND.  (RAL  .LT.  0.61087))  THEN 
C 

CYD AD«0. 00023894+ (0.00195121* (RAL-0. 55851001) )+ (0.02459273 
+*(RAL-0. 55851001) **2) - (0 . 1202244* ( (RAL-0 .55851001) **3 ) ) 

ENDIF 

C 

IF  (RAL  .GE.  0.61087)  THEN 
C 

CYDAD=0 . 27681285- (2 . 02305395*RAL)  +  (6 . 01180715*RAL*RAL) 
+-(9.24292188* (RAL**3) )+ (7.59857819* (RAL* *4) ) 

+-(2.8565527* (RAL**5) )+ (0.25460503* (RAL**7) ) 

+- (0 .01819815* (RAL**9) ) 

ENDIF 


IF  (RAL  .LE.  0  C  THEN 

EPA43-1.0 

ENDIF 

IF  (RAL  .GT.  0.0  AND  .LE.  0.6283185)  THEN 
0.6283185  RADIANS  «  36  DEGREES 

EPA43-0 . 9584809+ (4 . 13369452*RAL) -( 18 . 31288396*RAL*RAL) + 

+ (19 . 5511466* (RAL**3 ))- (1 .09295946*RAL*DSPBD) + (0 . 17441033* 
+DSPBD*DSPBD) 

ENDIF 

IF  (RAL  .GT.  0.6283185)  THEN 

EPA43sl.O 

ENDIF 


C  *  NOTE  -  THE  PARAMETER  EPA43  IS  A  MULTIPLIER  ON  RUDDER 

C  *  EFFECTIVENESS  DUE  TO  SPEEDBRAKE.  THIS  TABLE  IS  ALSO 

* 

C  *  LIMITED  TO  36  DEG  AOA.  HOWEVER,  THERE  IS  NO  AERODY 

# 

C  *  NAMIC  EFFECT  FOR  ANGLES  OF  ATTACK  LESS  THAN  16  DEG, 

# 

C  *  AND  THE  SPEEDBRAKE  IS  AUIOMATICALLY  RETRACTED  AT  AOA 
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C  *  GREATER  THAN  15  DEG.  THEREFORE,  THIS  TABLE  SHOULD 

* 

C  *  NOT  BE  NECESSARY  FOR  THE  ORDINARY  OPERATION  OF  THE 

* 

C  *  AIRCRAFT 

C 


c 

CYDRD-0 . 00310199+ (0 . 00119963*RAL)  ^ ( 0 . 02806933*RAL«RAL) 

+- (0 . 12408447* (RAL**3) ) - (0 . 12032121* (RAL**4) ) 

++ (0.79150279* (PAL**5) >-(0.86544347 *(RAL**6) ) 

++ (0.27845115* (RAL**7) )+(0.00122999*RAL*RARUD) +(0.00145943 
♦  *RARUD)-(0.01211427*RARUD*RARUD)  + (0.00977937*  (RARUD**3)  ) 

C 

CYDTD»-0. 00157745- {0.0020881*RAL)+(O.OOS57239*RAL*RAL) 

+- (0. 00139886* (RAL**3) >+(0.04956247* (RAL**4) ) 

+  - (0. 0135353* (RAL**5>  >-(0.11552397* (RAL**6> ) 

++ (0. 11443452* (RAL**7>  >-(0.03072189* (RAL**8> )- (C .01061113* 

+  (RAL**3)  *DSTBR)  -(0.00010529*RAL*RAL*DSTBR*DSTBR) 

+  - { 0 . 0 0 5 7 2 4 6 3  * RAL *DSTBR * DSTBR > 

++(0.01885361*RAL*RAL*DSTBR) -(0.01412258*RAI.*(DSTBR**3) ) 

+  - (0 .00081776*DSTBR)  +  (O . 00404354* {DSTBR**2 >  > - 
+ (0. 00212189* (DSTBR**3) >+(0.00655063* (DSTBR**4) > 

++(0. 03341584* (DSTBR**5)) 

C 

RALYl-0. 6108652 
RALY2- 9 0,0 /DEGRAD 
RBETYl«-0. 0872665 
RBETY2«0. 1745329 
C 

AY-0.164 
ASTARY-0. 95993 
BSTARY-0, 087266 
C 

ZETAY- (2 .ODO* ASTARY- (RALY1+RALY2 ) > / (RALY2-RALY1 > 

ETAY- ( 2 . 0D0*BSTARY- (RBETY1+R3ETY2 ) ) / (RBETY2 -RBETYl ) 

C 

X- ( 2 . ODO *RAL- ( RALY1+RALY2 ) > / ( RALY2 -RALYl ) 

Y- (2 . ODO *RBETA- (RBETY1+RBETY2 ))/ (RBETY2 -RBETYl) 

C 

FY-( {5.0D0*(ZETAY**2> )-(4.0D0*ZETAY*X)-1.0D0)*( ( {X**2) -l.ODO) 
+**2) *(1.0D0/ ( ( (ZETAY**2) -l.ODO) **3) ) 

C 

GY-( (5.0D0*(ETAY**2> )-(4.0D0*BTAY*Y) -l.ODO) * ( ( ( Y**2 ) -1 . ODO ) **2) 
+*(1.0DC/( ( (ETAY**2) -l.ODO) **3) ) 

C 

CYRB»AY*FY*GY 

C 

IF  (RAL  .LT,  0.6108652)  THEN 
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c 

CYRBaO.O 
GOTO  500 
ENDIF 
C 

IF  ((RBETA  .LT.  -0.0872665)  .OR.  (RBETA  .GT.  0.1745329))  THEN 
C 

CYRB-0.0 
GOTO  500 
ENDIF 
C 

500  CFY« (CFY1*EPA02L) ♦ (CYDAD*DAILD) ♦ (CYDRD*DRUDD*DRFLX5*EPA43 ) + 

•►{  (CYDTD*DTFLX5)*DTALD)  +  (CFy..  'PB)  ♦  (CFYR*RB) 

♦♦CYRB 


«  •  «  « 


DTFLXl-0.975 

DRFLXl-0.85 

C 

CMLl.-0 . 00238235- (0 .04616235*RAL)  •!•( 0 . 10553168*RAL*RAI-) 

♦♦ (0 . 10541585* (RAL**3) ) - (0.40254765* (RAL**4) ) 

♦♦ (0. 32530491* (RAL**5) )- (0.08496121* (RAL**6) ) 
♦♦(0.00112288*(RAL**7))-(0.05940477*RABET*RAL)- 
♦  (0.07356236*RABET)  -(0.00550119*RABET*RABET) +  (0.00326191 
♦*(RABET**3) ) 

C 

IF  (RAL  .LT.  0.29671)  THEN 
C 

CMLP.-O . 24963201- (0 .03106297*RAL) + (0 . 12430631*RAL*RAL) 
+-(8.95274618* (RAL**3) )♦ (100. 33109929* (RAL**4) ) 

++ <275. 70069578* (RAL**5) )- (1178.83425699- (RAL**6) ) 

+- (2102. 66811522* (RAL**7) )♦ (2274. 89785551* (RAL**8) ) 

ENDIF 

C 

IF  ((RAL  .GE.  0.29671)  .AND.  (RAL  .LT.  0.34907))  THEN 
C 

CMLP--0. 1635261- (3 .77847099* (RAL-0 . 29671001 ))+( 147 . 47639465 
+* (RAL-0. 29671001 ) **2)- (1295- 94799805* (RAL-0. 29671001) **3) 

ENDIF 

C 

IF  (RAL  .GE.  0.34907)  THEN 
C 

CMLP--1. 37120291+ (7. 06112181*RAL) -(13. 57010422*RAL*RAL) 

++( 11. 21323850* (RAL**3) ) 

+- (4. 26789425* (RAL**4) )+ (0.6237381* (RAL**5) ) 

ENDIF 
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IF  (RAL  .LT.  0.7854)  THEN 
C 

CMLR.0 .03515391<*-(0.59296381*RAL>  +  (2 .27456302*RAL*RAL) 
t-(3.8097803*(RAL**3) ) 

(45. 83162842*  (RAL**4)  )•►( 55. 31669213*  (RAL**5)  )♦ 

♦ (194.29237485* (RAL**6))- (393.22969953* (RAL**7) )+ (192. 20860739* 
♦(RAL**8) ) 

ENDIP 

C 

IP  ((RAL  .GE.  0.7854)  .AND.  (RAL  .LE.  0.87266))  THEN 
C 

Ca4LR-0. 0925579071- (0. 6000000238* (RAL-0. 7853999734) ) 

(1. 3515939713*  ( (RAL-0. 7853999734 )**2) ) 

++(29.0733299255* ((RAL-0. 7853999734) **3) ) 

ENDIP 

C 

IF  (RAL  .GT.  0.87266)  THEN 
C 

CMLR«-311 . 126041+ (1457 .23391042*RAL) - (2680 . 19461944*RAL*RAL) + 

+ (2361.44914738* (RAL**3)) -(893.83567263* {RAL**4) >+(68.23501924* 
+(RAL**6) )-(1.72572994*(RAL**9)) 

ENDIP 

C 

CLDAD.0.00057626+(0.00038479*RAL)-(0.00502091*RAL*RAL) 
++(0.00161407* (RAL**3) >+(0.02268829* (RAL**4) ) 

+- (0. 03935269* (RAL**5) >+(0.02472827* (RAL**6) ) 

+- (0. 00543345* (RAL**7) )+(0.0000007520348*DSTBR*RAL)+ 
+(0.000000390773*DSTBR) 

C 

CLDRD.0.00013713-(0.00035439*RAL)-(0.00227912*RAL*RAL) 

++ (0. 00742636* (RAL**3) >+(0.00991839* (RAL* *4) ) 

+- (0.04711846* (RAL* *5) ) +(0.046124* (RAL**6) ) 

+- (0. 01379021* (RAL**7) )+(0.00003678685*RARUD*RAL)+ 

+ (0 . 00001043751*RARUD) - (0 . 00015866*RARUD*RARUD) + (0 . 00016133 
+*(RARUD**3) ) 

C 

CLDTD-0. 00066663+ (0.00074174*RAL)+(0.00285735*RAL*RAL) 
+-(0.02030692* (RAL**3))-(0. 00352997* (RAL**4) ) 

++ (0. 0997962* (RAL**5) >-(0.14591227* 

+(RAL**6) >+(0.08282004* (RAL**7) ) 

+- (0 . 0168667* (RAL**8) )+ (0.00306142* (RAL**3) *DSTBR) 

+- (0 . 00110266*RAL*RAL* (DSTBR**2 ) ) + (0 . 00088031*RAL* 
+(DSTBR**2) )-(0.00432594*RAL*RAL*DSTBR)- 
+(0.00720141*RAL*(DSTBR**3) ) 

+- (0.00034325*DSTBR)+(0. 00033433* (DSTBR**2) >+(0.00800183 
+*(DSTBR**3) ) -(0.00555986*(DSTBR**4) )-(0.01841172*(DSTBR**5) ) 

n 

IF  (RAL  .LT.  0.0)  THEN 
C 
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DCLB>i-0. 00006 
ENDIF 

IF  ((RAL  .GE.  0.0)  .AND.  (RAL  .LE.  0.209434))  THEN 
DCLB— 0.00006+ (0.0041035078*RAL*RAL)- (0.0130618699* (RAL**3) ) 

ENDIF 

IF  (RAL  .GT.  0.209434)  THEN 

DCLB>0 . 0 
ENDIF 


CML-  (CML1*EPA02S)  ♦  (CLDAD*DAILD)  +  (CLDRD*DRUDD*DRFLX1*EPA43 )  + 
♦  (  (CLDTD*DTFLX1)  •DTALD)  +  (CMLP*PB)  +  (CMLR  RB)  +  (DCLB*BETA) 


C 


CMMl-0. 00501496- (0. 08004901*RAL) -(1.03486675*RAL*RAL) 
+-(0.68580677* (RAL**3) )+ (6,46858488* (RAL**4) ) 

+- (10. 15574108* (RAL**5))+ 

+ (6. 44350808* (RAL**6) )- (1.46175188* (RAL**7) ) 
++(0.24050902*RAL*DSTBR) 

+- (0 . 42629958*DSTBR) - (0 .03337449*DSTBR*DSTBR) 

+-(0.53951733* (DSTBR**3) ) 

C 

IF  (RAL  .LE.  0.25307)  THEN 
C 

CMMQ— 3. 8386262+ (13. 54661297*RAL)+ (402. 53011559*RAL*RAL) 

♦- (6660. 95327122* (RAL**3) )- (62257.89908743* (RAL**4) ) 
++(261526.10242329* (RAL**5) ) 

++ (2177190. 33155227* (RAL**6) )- (703575. 13709062* (RAL**7) ) - 
+(20725000. 34643C54*(RAL**8) )-(27829700.53333649*(RAL**9) ) 
ENDIF 
C 

IF  ((RAL  .GT.  0.25307)  .AND.  (RAL  .LT.  0.29671))  THEN 
C 

CMMQ--8. 4926528931- (2705. 3000488281* (RAL-0 .2530699968) ) 

++( 123801. 5* (RAL-0. 2530699968) **2) 

+-(1414377* (RAL-0. 2530699968 ) **3) 

ENDIF 

C 

IF  (RAL  .GE.  .29671)  THEN 
C 

CMMQ-47 .24676075- (709 .60757056*RAL) + (3359 . 08807193 *RAL*RAL) - 
+  (7565. 32017266* (RAL**3) )  + (8695. 1858091* (RAL**4) ) 

+- (4891. 77183313* (RAL**5) )+ (1061.55915089* {RAL**6) ) 
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c 


ENDIF 

CMM-CMM1+ (CMMQ*QB) 


C 


DTFLX3. 0.975 
DRFLX3a0.89 


CMNl.O. 01441512f(0-02242944*RAL)-(0. 30472558* (RAL**2) ) 
♦♦(0.14475549* {RAL**3)) 

♦♦ (0. 93140112* (RAL**4) >-(1.52168677* (RAL**5) )♦ 

♦ (0. 90743413* (RAL**6) >-(0.16510989* (RAL**7) > 

♦  - (0.0461968* (RAL**8>  > 

♦  ♦ (0 . 01754292* {RAL**9>  > - (0 . 17553807 *RAL*RABET> ♦ 

♦ ( 0 . 15415649*RAL*RABET*DSTBR) 

♦♦ (0. 14829547* (RAL**2>*{RABET**2>  > 

♦  - (0 . 11605031* (RAL**2>  *RABET*DSTBR> 

♦  - (0. 06290678* (RAL**2 >* (DSTBR**2>  > 

♦- (0 . 01404857* (RAL**2 >* (DSTBR**2 > > 

♦♦ (0.07225609*RABET>-(0. 08567087* (RABET**2 >  > 

♦♦ (0.01184674* (RABET**3>  > 

♦-(0.00519152*RAL*OSTBR>^(0.03865177*RABET*D3TBR) 
♦♦(0.00062918*DSTBR> 

C 

CNDRD«-0 .00153402^ (0 .00184982*RAL> - (0 .0068693*RAL*RAL> 

♦♦ (0. 01772037* (RAL**3 >  > 

♦  ♦(0. 03263787* (RAL**4>  >- (0 . 15157163* (RAL**5>  >♦ (0 . 18562888 

♦  *(RAL**6> )- (0.0966163* (RAL**7>  >♦ (0.01859168* (RAL**8>  >♦ (0.0002587 
♦*RAL*DSTBR> - (0 .00018546*RAL*DSTBR*RBETA> - (0 .00000517304*RBETA> 

♦  - (0 .00102718*RAL*RBETA> - (0 . 0000689379*RBETA*DST8R> - (0 . 00040536 

♦*RBETA*RARUD> - {0.00000480484*OSTBR*FARUD> 

♦- (0 . 00041786*RAL*RARUD> 

♦♦(0.0000461872*RBETA>^{0. 00434094* (RBETA**2>  > 

♦  - (0.00490777* (RBETA**3>  > 

♦♦ (0 . 000005157867*RARUD> ♦ (0 . 00225169*RARUD*RARUD> - (0 . 00208072 

♦  *(RARUD**3>  > 

C 

IF  (RAL  .LT.  0.55851>  THEN 
C 

CMNP--0 . 00635409- (1 . 14153932*RAL) ♦ (2 . 82119027* (RAL**2 )> ♦ 

♦  (54.4739579*  (RAL**3>  >-(140.89527667* (RAL**4>  >-(676.73746128* 

♦  (RAL**5>  >♦ (2059.18263976* (RAL**6>  >♦ (1579. 41664748*  (RAL**7>  > 
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♦- (8933. 08535712* (RAL**8 ))♦ (6806. 54761267* (RAL**9) ) 

ENDIF 

C 

IF  ((HAL  .GE.  C. 55851001)  .AND.  (RAL  .LE.  0.61087)}  THEN 
C 

CMNP«-.07023239-*-(l  .085815*(RAL  -0.55851)  ) 

♦♦ (8. 852651* ( (RAL-. 55851) **2) )- (192.6093* ( (RAL-0. 55851) **3) ) 
ENDIF 
C 

IF  (RAL  .GT.  0.61087)  THEN 
C 

CMNP=-71 . 03693533+ (491 . 32506715*RAL) 

♦-(1388.11177979* (RAL**2) )♦ 

♦ (2033. 48621905* (RAL**3 ) ) 

*-(1590.91322362* (RAL**4) )♦ (567. 38432316* (RAL**5) ) 
♦-(44.97702536* (RAL**7) )♦ (2. 8140669* (RAL* *9) ) 

ENDIF 


IF  (RAL  .LE.  -.069813)  THEN 
C 

CMNR-  -0.2805 

ENDIF 

C 

IF  ((RAL  .GT.  -.069813)  .AND.  (RAL  .LT.  0.0))  THEN 
C 

CMNR--0. 2804999948* (35. 9903717041* (RAL* . 0698129982 ) **2 ) 
♦-(516.1574707031* (RAL*. 0698129982 ) **3) 

ENDIF 

C 

IF  ((RAL  .GE.  0.0)  .AND.  (RAL  .LE.  0.78539801))  THEN 
C 

CMNR--. 28071511- ( 2. 52183924*RAL)*( 68. 90860031* (RAL**2) ) 

♦- (573.23100511* (RAL**3) )+ (2009.08725005* (RAL**4) ) 
♦-(3385.15675307* (RAL**5) ) 

♦♦ (2730. 49473149* (RAL**6))- (848. 12322034* (RAL**7 ) ) 

ENDIF 

C 

IF  ((RAL  .GT.  0.78539801)  .AND.  (RAL  .LT.  0.95993102))  THEN 


C 

<n4NR=-0.1096954+(0.52893072*(RAL-0.78539801))-(6.09109497*(RAL- 
♦0 .78539801) **2> ♦ (17 .47834015* (RAL-0 .78539801) **3) 

ENDIF 

C 

IF  (RAL  .GE.  0.95993102)  THEN 
C 

CMNR»-0.11 

ENDIF 

C 
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CNDTD-0 . 00058286+ (0 . 0007341*RAL) - (0 . 00746113 *RAL*RAL) 

+- (0.00685223* (RAL**3) ) 

++ (0.03277271* (RAL**4) )- (0.02791456* (RAL**5) ) 

++ (0.00732915* (RAL**6) ) 

♦+ (0. 00120456*RAL*DSTBR)-(0. 00168102 *DSTBR)+( 0.0006462* 
+DSTBR*DSTBR) 

C 

CNDAD=0 . 00008228887- (0 .00014015*RAL) - (0 . 0013493 *RAL*RAL) + 
♦(0.00020487* (RAL**3) )+ (0.00561241* (RAL**4) ) 

♦-(0. 00634392* (RAL**5) ) 

+♦(0.00193323* (RAL*  6) ) - (2 . 05815E-17* (RAL*DAILA) ) + (3 .794816E-17* 
+(DAILA**3) ) 

C 

DCNB=t-2.500E-4 

C 

RALNlsO. 69813 
RALN2= 9 0.0 /DEGRAD 
RBETNl=-0. 174532 
RBETN2=0. 34906 
C 

ANsO.034 
ASTARNal.0472 
BSTARN=0. 087266 
C 

ZETAN=(2 .0D0*ASTARN-(RALN1+RALN2)  )  /  ( RALN2 -RALNl ) 

ETANs(2 .0D0*BSTARN-(RBETN1+RBETN2)  )  /  (RBETN2 -RBETNl ) 

C 

XN=  (2  .  ODO*RAL- (RALN1+RALN2 )  )  /  (RALN2-RALN1 ) 

YN=(2 .0D0*RBETA-(RBETN1+RBETN2) ) / ( RBETN2 ‘RBETNl ) 

C 

FN»( (5.0D0* (ZETAN**2) )-(4.0D0*ZETAN*XN) -l.ODO) * 

+( ( (XN**2) -l.ODO) **2)/( ( (ZETAN**2) -l.ODO) **3) 

C 

GN»( (5.0D0*(ETAN**2) ) -(4.0D0*ETAN*YN) -l.ODO) * 

+(( (YN**2) -l.ODO) **2)/( ( (ETAN**2) -l.ODO) **3) 

C 

CNRB=AN*FN*GN 

C 

IF  (RAL  .LT.  0.69813)  THEN 
C 

CNRBsO.O 
GOTO  1000 
ENDIF 
C 

IF  ((RBETA  .LT.  -0.174532)  .OR.  (RBETA  .GT.  0-34906))  THEN 
C 

CNRBsO.O 
GOTO  1000 
ENDIF 
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C 

1000  CMN=(CMNl*EPA02S)  +  (CNDAD*DAILD>  +  (  (CNDRD*DRUDD*DRFLX3)  *EPA43)  + 

+  {  (CNDTD*DTFLX3 )  *DTALD)  (CMNP*PB)  +  {CMNR*RB)  (DCNB*BETA) 
++CNRB 


C 


CX*CFZ*SIN (RAL) -CFX*COS (RAL) +THRUST/QBARS 
CYsCFY 

CZ=-(CF2*COS{RAL)+CFX*SIN(RAL) ) 

CLM=CML 

CMMrCMMfTHRUST* (0. 25/12. 0)/( QBARS  *CW1NG ) 

THE  (0.25/12.0)  IS  THE  OFFSET  OF  THE  THRUST  VECTOR  FROM  THE 
CG 

CNM=CMN 

RETURN  CX,  CY,  CZ.  CLM,  CMM,  CNM  TO  CALLING  PROGRAM. 


IWRITEsO 

IF(IWRITE.EQ.l)  WRITE(6,*)'  CX, CY, CZ, CLM, CMM, CNM= ' , 

CX , C Y , CZ , CLM , CMM , CNM 


WRITE { 6 , * )  '  END  OF  SUB  COEFF  U(1,2,3)*',U(1),U(2),U(3) 
WRITE ( 6 , * ) ' LEAVING  COEFF ' 

RETURN 

END 
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n  o 


c 

program  integ 
c 

c  ordinary  differential  equations  propagator 
c  using  routine  'haming*  as  supplied  by  w.  wiesel 
c 

implicit  double  precision (a-h, 0-2) 
common  /ham/  t,x(12,4) , f (12,4) ,errest (12) ,n.h 
common  /bpars/  PAR (10) 
open (15, files'  fort .15 ' ) 
open (7 , files ' fort .7 ' ) 
write  (*,*)  'help  me' 
c 

c  read  in  max  printout,  number  of 

c  steps  between  printing 

c 

n  «  8 
npar  »  3 

read  (15,*)  nprint,nstep 
write  (*,*)  'I  M  trapped' 
c 

c  read  timestep 

c 

t  =  O.D  00 
read  (15,*)  h 
write  (*,*)  'in  the  body' 
write  {*,*)  'hs',h 
c 

c  read  initial  conditions  and  parameter  values 

c 

read  (15,*)  (x(i,l) , isl,n) 
write  (*,*)  'of  an  aeronautical  engineer' 
read  (15,*)  ( PAR ( i) , isi, npar) 
write  (*,*)  'I  m  really' 
c 

c  initialize  haming 

c 

nxt  3  0 

call  haming (nxt) 
if (nxt  .ne.  0)  go  to  50 
write  (*,1) 

1  format (2x,'  haming  did  not  initialize') 

stop 
c 

50  continue 


write  (*,*)  'a  wily  coyote' 
c 

c  integrate  ode.... two  nested  loops 
c 


112 


n o o n o o o o o o n n o n  on  o  o o o o o  o  no 


do  200  ipr  s  l,n*'.rint 
c 

do  100  istp  >  l.nstep 
c 

c  each  call  to  haming  advances  one  step... 
c 

call  haming (nxt) 
c 

100  continue 
c 

c  after  nstp  integration  steps,  print  current  values 
c 

write(6,2)  t, (x{ii,nxt) , ii=l,n) 

2  format ( lx, 9ei4 .7) 
c 

200  continue 


scop 

end 

subroutine  rhs{nxc) 

rhs  is  the  right  -  hand  -  side  subroutine,  customized 
to  the  particular  sec  of  odes  being  integrated 


IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 

COMMON  /KS/  K1,K5,K7,K8,K9,K10,K12,K13,K14,K15,K16,K17 

COMMON  /ACDATA/  SWING, CWING,  SREF,  RHO,  W  AS S,  THRUST 

DOUBLE  PRECISION  Kl, K5, K7 , K8 , K9 , KIO , K12 , K13 , K14 , K15, K16 ,  K17 

DOUBLE  PPECISION  K2, K3 , K4 , K6, Kll , IX, lY, IZ, IXZ 

COMMON  /SEIZE/  CX,Cy,CZ,CLM,CMM.CNM 

common  /ham/  t,xin(12,4) , f (12,4) ,err(12)  ,n,h 

common  /jpars/  PAR (10) 


DIMENSION  U(20),FF(20) 
WRITE{6,*)  'ENTERED  FUNC' 


INITIALIZE  SOME  CONSTANTS  THAT  ARE  PASSED  THROUGH 
THE  COMMON  BLOCK  ACDATA 

DATA  IS  FROM  MCAIR  REPORT#  A4172  AND  AFFTC-TR-75-32 
F-15A  APPROACH-TO-STALL/ STALL/ POST-STALL  EVALUATION 

BWING  -  A/C  WINGSPAN,  FT 

CrWING  -  A/C  MEAN  AERODYNAMIC  CHORD,  FT 

SREF  -  A/C  WING  REFERENCE  AREA,  SQ  FT 

RHO  -  AIR  DENSITY  AT  20000  FT  ALTITUDE,  SLUG/FT''3 
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RMASS  -  A/C  MASS,  SLUGS 
THRUST  -  TOTAL  A/C  THRUST,  LBS 

BWINGb42.8D  00 
CWING>15.94D  00 
SREF«608.D  00 
KHOs. 00126730  00 
RMASS-37000.D  00/32. 174D  00 
THRUST.8300.0  00 


DETERMINE  CONSTANTS  Kl  THROUGH  K17 .  SOME  ARE  MADE  COMMON 

PASSED  TO  SUBROUTINE  FUNC  AND  USED  IN  THE  EQUATIONS 
OF  MOTION  THERE 

INERTIAS  HAVE  UNITS  OF  SLUG-FT"2 
Kl  HAS  UNITS  OF  1/FT 

K6,  K8,  Kll,  K14,  AND  K17  HAVE  UNITS  OF  1/FT^2 


IX«  25480. 

IY=  166620. 

I2«  186930. 

IXZ»  -1000. 

Kl= . 5*RHO*SREF/RMASS 
K2»(IZ-IY) /IX 
K3«IXZ*IX2/(1X*1Z) 
K4»(IY-IX) /IZ 
KS-IXZ/IX 

K6o.5*RHO*BWlNG*SREF/IX 

K7.IXZ/IZ 

K8».5*RHC  .-REF*CWING/iy 
K9«(IZ-IX) /lY 
KIO-IXZ/IY 

Kll*.5*RHO'SREF*BWING 

K12=(K2fK3)/(l.-K3) 

K13=(1.-K4)*K5/(1.-K3) 

K14.K6/ (1.-K3) 

K15-(K3-K4)./|1.-K3) 

K16-(l.+K2) *K7/(1.-K3) 

K17-K11/(1.-K3) 


Kl  «  3 .3500888900-04 
K5  «-3.924646781D-02 
K7  .-5.349596105D-03 
K8  .  3 .6856509710-05 
K9  ■  .96897131196D  00 
KIO  --6.001680471D-03 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 
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c 

c 

c 

c 

c 

c 


K12  .  .79747314581D  00 
K13  S-9.615755341D-03 
K14  «  6.472745847D-04 
K15  S-.754990553922D  00 
K16  »  K13 

K17  s  8.822851S58D-05 


FIND  THE  VALUES  OF  F(l)  THROUGH  F(NDIM)  .  SUBROUTINES 
COEFF  AND  FUNX  ARE  CALLED  ONCE. 

U(l)  -  ALPHA,  DEG 
U(2)  -  BETA,  DEG 
U(3)  -  P,  RAD/SEC 
U(4)  -  Q,  RAD/SEC 
U{5)  -  R,  RAD/SEC 
U(6)  -  THETA,  DEG 
U(7)  -  PHI,  DEG 

U(8)  -  TRUE  VELOCITY,  IN  THOUSANDS  OF  FT/SEC 
PAR ( 1 ) -DELESD 

PAR(2)«ORUDD  THE  PARAMETERS,  IN  DEGREES 

PAR(3)«DDA 


do  12  isl,n 
u ( i ) Bxin ( i , nxt ) 
continue 

TEST  ALPHA  AND  BETA  FOR  OUT  OF  BOUNDS  CONDITION. 

•8  .GE.  ALPHA  .LE.  50  DEGREES  AND  -30  .GE.  BETA  .LE.  30 
DEGREES.  IF  A  VARIABLE  PASSES  OUTSIDE  A  BOUNDARY,  SET  THAT 
VARIABLE  EQUAL  TO  THAT  LIMIT  AND  ISSUE  A  MESSAGE  TO  THE 
SCREEN  AND  CONTINUE  TO  EXECUTE. 


c 

c  IF(U(1)  .GT.  50.)  WRITE(6,*)  'ALPHA  >  50  DEGREES,  CONTINUE 

TO  RUN' 

C  IF{U(1)  .LT.  -8.)  WRITE(6,*)  'ALPHA  <  -8  DEGREES,  CONTINUE 

TO  RUN' 

C  IF(U{2)  .GT.  30.)  WRITE;6,*)  'BETA  >  30  DEGREES.  CONTINUE  TO 

RUN' 

c  IF(U(2)  .LT.  -30.)  WRITE(6,*) 'BETA  <  -30  DEGREES,  CONTINUE 

TO  RUN' 

C 

C 

c  WRITE (*,  *)  'ALPHA, BETA,  P.Q- '  ,U(1)  ,U (2)  ,U{3)  ,Ui4) 

C  WRI'TEi*,*)  'R,'rHETA,PHI,VELOCITY»',U{5)  ,U(6;  ,U(7)  ,U(8) 

c 

CALL  COEFF (U, PAR, n) 

C 

CALL  FUNX(n,U,FF) 
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c 


13 


c 

c 

c 

C 

c 

C 

c 

c 

c 

C 

C 

C 


do  13  isl.n 
f (i.nxt) sFF{i) 
contim^e 


WRITE (*.*) 
WRITE (11,*) 
WRITE ( * , * ) 
WRITEdl,*) 
WRITE {*. *) 
WRITEdl,  *) 
WRITE (*,*) 
WRITEdl,  *) 


'Fl,F2a  ',FFd),FF(2) 
'F1,F2=  ',FFd),FF(2) 
'F3,F4*  ',FF(3) ,FF(4) 
'F3,F4=  ',FF(3) ,FF(4) 
'F5,F6«  ',FF(5) ,rF(6) 
'F5,F6*  ',FF(5) ,FF{6) 
'F7,F8»  ' ,FF(7) ,FF(8) 
'F7,F8=  ',FF(7) ,FF(8) 


RETURN 

END 


SUBROUTINE  FUNX (NDIM, U, F) 


SUBROUTINE  FUNX  EVALUATES  THE  NDIM  EQUATIONS  GIVEN  THE 
STATE  VECTOR  U. 

NDIM-  THE  DIMENSION  OF  THE  PROBLEM 
U  -  THE  VECTOR  OF  STATES  ALPHA,  BETA,  ...  (INPUT) 

F  -  THE  VECTOR  RESULT  OF  FUNCTION  EVALUATIONS  (OUTPUT) 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-2) 

COMMON  /SEIZE/  CX,CY,CZ,CLM,CMM,CNM 

COMMON  /KS/  K1,K5,K7,K8,K9,K10,K12,K13,K14,K15,K16,K17 
DOUBLE  PRECISION  Kl, K5, K7 , K8, K9 , KIO ,  K12 , K13 ,  K14 , K15 , K16, K17 
DIMENSION  U (NDIM) ,F (NDIM) 

SET  TRIGONOMETRIC  RELATIONSHIPS  OF  THE  STATES  ALPHA,  BETA, 
THETA,  AND  PHI  AND  THEN  SET  P,  Q,  R,  AND  VTRFPS 

WRITE(6,*)  'ENTERED  FUNX' 

IWRITE-1 

IF(IWRITE.EQ.1)WRITE(6, *) ' 
,5,7,8,9,10,12-17.' ,K1,K5,K7,K8, 

+  K9,K10,K12,K13,K14,K15,K16,K17 

IF(IWRITE.EQ.1)WRITE(6,*) 'FUNX  U(l)-U(8).',  1 ) , U ( 2 ) , U( 3 ) , 

♦  U(4) ,U(5) ,U{6) ,U(7) ,U(8) 


DEGRAD-57. 29577951D  00 

CA.DCOS (U (1 ) /DEGRAD) 

SA. DSIN(U(1) /DEGRAD) 
CB-DCOS(U (2) /DEGRAD) 

SB. DSIN(U(2) /DEGPAD) 
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CTHE-DCOS (U ( 6 ) /DEGRAD) 

S‘niE«DSZN  (U ( 6 )  /DEGRAD) 

CPHI-DCOS  (U(7 )  /DEGRAD) 

SPHZ»0SIN(U(7 ) /OEGRAD) 

C 

P-U(3) 

Q-U(4) 

R>U(5) 

\TRFPS«1000.D  00*U(8) 

C 

C  SET  THE  GRAVITATIONAL  CONSTANT.  FT/SEC 
C 

Ga32.174D  00 
C 

C  THE  FOLLOWING  SYSTEM  OF  NONLINEAR  DIFFERENTIAL  EQUATIONS 
C  GOVERN  AIRCRAFT  MOTION 

C 

C  UPDATED  FOR  PROPER  DEGREE-RADIAN  UNITS  AND  PROPERLY 
C  SCALED  VELOCITY  EQUATION:  7  JUN  88 
C 

C  F ( 1 ) .ALPHA-DOT 

C 

1  F  ( 1)  .Qf  ( -  (K1*VTRFPS*CX-G*STHE/VTRFPS+R*SB)  *SAf  (K1*VTRFPS 

-*■  *CZ>  (G*CTHE*CPHI/VTRFPS)  -P*SB)  *CA)  /CB 

F(1)«F(1)*DEGRAD 

C 

C 

C  F (2) -BETA-DOT 
C 

2  F(2)--(  (K1*VTRFPS*CX-G*STHE/VTRFPS)  *SB+R)  •CA+ (K1*VTRFPS*CY 

♦  ♦G*CTHE*SPHI/VTRFPS)  *CB-  ( (K1*VTRFPS*CZ+G*CTHE*CPHI/VTRFPS) 

♦  •SB-P)*SA 
F(2)-F(2)*DEGRAD 

C 

C 

WRITE (6,*)  'K1,CX,CY,CZ,F(2),G,SB,CB,CA,SA,SPHI,CPHI,CTHE.STHE' , 

C  ♦ 

K1,CX,CY,CZ,F(2) .G.SB.CB.CA, SA.SPHI.CPHI.CTHE.STHE 
C 

C  FI3)»P-DOT 
C 

3  F(3)  «-K12*Q*R+K13*P*Q+Kl4*(CLM+K7*CNM)  ‘VTRFPS'VTRFPS 
C 

C 

C  F(4)-Q-DOT 
C 

4  F(4)-K8*VTRFPS*VTRFPS*CMM^K9*P*R<fKlO*(R*R-P*P) 

C 

C 

C  F(5)-R-DOT 
C 

5  F(5)  «K15*P*Q-K16*0*R-»-K17»VTRFPS*VTRFPS*{K5*CLM+CNM) 
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F { 6 ) -THETA-DOT 


6  F(6)-Q*CPHI-R*SPH1 
F(6)-F(6) •DEGRAO 


F (7) -PHI-DOT 

7  F(7)-P+Q*  (STHE/CTHE)*SPHI+R*  (STHE/CTHE)  *CPHI 

F(7)-F(7)*DEGRAD 


F(8) -VTRFPS-OOT  (SCALED  BY  A  FACTOR  OF  1000) 

8  F  ( 8 )  -0  ( 8 )  *  ( { K1  * VTRFPS*CX-G*STHE/VTRFPS )  •CA*CB+  ( K1  *VTRFPS*CY 

♦ 

;*CTHE*SPHI/VTRFPS)  *SB+ (K1*VTRFPS»CZ4.G*CTHE*CPHI/VTRFPS)  *SA*CB) 


WRITE (6,*)  'LEAVING  FUNX' 

WRITE(6.*)'lv  funx 
',F(1)  ,F(2)  ,F|3)  ,F(4)  ,F(5)  .F(6)  .F(7)  ,F(8) 

RETURN 

END 


subroutine  haming(nxt) 
c 

c  haming  is  a  fourth  order  predictor-corrector  algorithm 
c  for  the  integration  of  systems  of  ordinary  differential 

equations 

c  the  common  block  /ham/  contains  most  of  the  variables: 
c  X  is  the  independent  variable,  the  'time' 

c  y  contains  4  copies  of  the  state  vector,  with 

c  n  odes  being  integrated 

c  f  contains  the  calculated  equations  of  motion 

c  errest  is  a  truncation  error  estimate 

c  n  is  the  number  of  odes 

c  h  is  the  integration  timestep 

c  nxt  assumes  the  values  1,2,3,4.1,2,3,4, . . . .  and  points  to 
c  the  current  value  of  the  state  vector 
c 

c  the  user  must  supply  a  subroutine  'rh8(nxt)'  which 
c  calculates  the  equations  of  motion  f(i,nxt)  from  the 
c  state  vector  y{i,nxt) 
c 

c  to  initialize  haming,  the  initial  conditions  must  be  stored 
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c  in  y(i,l).  i«l,n  ;  x,n,  and  h  must  be  initialized,  and 
c  then  haming  is  called  with  nxtaO.  If  haming  returns  with 

c  nxtsl,  initialization  is  successful.  If  nxt«0  still,  haming 

c  did  not  initialize  (h  is  usually  too  big) 
c 

comnvon  /ham/  x.y(12,4) ,  f  (12,4)  ,errest  (12)  ,n,h 
double  precision  x,y, f ,errest,h,xo, tol,hh 
tol  ■  1.0d'-12 

c  branch  on  nxt:  startup  or  propagating? 

if(nxt)  190,10,200 
c 

c  haming  initialization:  4  point  picard  iteration 
c 

10  XO  S  X 

hh  =  h/2.0 
call  rhsd) 
do  40  1  «  2,4 

X  ■  X  hh 

do  20  i  3  i.n 

20  y(i,l)  .  y(i,l-l)  ♦  hh«f(i,l-l) 
call  rhsd) 

X  3  X  ♦  hh 

do  30  i  3  l,n 

30  y(i,l)  3  y(i,l-l)  ♦  h*f(i,l) 

40  call  rhsd) 
jsw  3  -10 
50  isw  3  1 

do  120  i  3  l,n 

hh  3  y(i,l)  ♦  h*(  9.0*f(i,l)  ♦  19.C*f(i,2)  -  5.0*. (i, 3) 

1  ♦  f(i,4)  )  /  24.0 

if(  dab8(  hh-y(i,2))  .It.  tol  )  go  to  70 
isw  3  0 

70  y(i,2)  3  hh 

hh  3  y(i,l)  ♦  h*(  f(i,l)  +  4.0*f{i,2)  ♦  f(i,3))/3.0 
if(  dab8(  hh-y(i,3))  .It.  tol  )  go  to  90 
isw  3  0 

90  y(i,3)  3  hh 

hh  3  y(i,l)  ♦  h*(  3.0*f(i,l)  i-  9.0*f(i,2)  +  9.0*f(i,3) 

1  3.0*f(i,4)  )  /  8.0 

if(  dabs(hh-y(i,4) )  .It.  tol  )  go  to  110 
isw  3  0 

no  y(i,4)  3  hh 
120  contin".e 

X  3  XO 

do  130  1  3  2,4 

X  3  X  +  h 

130  call  rhsd) 

ifdsw)  140,140,150 
140  jsw  3  jsw  *  1 

if(jsw)  50,280,280 
150  X  3  XO 
isw  3  1 
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jsw  >  1 

do  160  i  >  l.n 
160  err«8t(i)  ■  0.0 
nxt  ■  1 

go  to  280 

190  jsw  >  2 

nxt  ■  isbsCnxt) 

c 

c  hsming  propagation  saction 

200  X  ■  X  h 

npl  a  inod(nxt,4)  *  1 
go  to  (210,230) , law 
210  go  to  (270,270,270,220) , nxt 
220  isw  a  2 

pennuta  indicas 

230  nm2  a  mod(npl,4)  -«■  1 

nml  a  mod(nm2,4)  1 

npo  a  mod  ( nml ,  4 )  1 

pradictor 

do  240  i  a  l,n 

£(i,nm2)  a  y(i,npl)  ♦  4.0*h*(  2.0*£(i,npo)  -  f(i,ninl) 

1  ♦  2.0*f (i,nin2)  )  /  3.0 

240  y(i,npl)  a  f(i,nm2)  -  0.925619835*arra8t (i) 
call  rh8(npl) 

corractor 

do  250  i  a  i,n 

y(i,npl)  a  (  9.0*y(i,npo)  -  y(i,nin2)  *  3.0*h*{  f(i,npl) 
1  ♦  2.0*£(i,npo)  -  £(i,nml)  )  )  /  8.0 

arrastd)  a  f(i,nm2)  -  y(i,npl) 

250  y(i,npl)  a  y(i,npl)  ♦  0.0743801653  *  arrastd) 
go  to  (260,270) , jsw 
260  call  rh8(npl) 

270  nxt  a  npl 
280  ratum 
and 


SUBR(3UTINE  COEFF  (U,ND1M) 

C  Last  Editad:  30  Kay  90  by  J.Bac]c 

C  This  subroutine  computes  the  aerodynamic  force  and  moment 
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C  coefficients  to  be  used  in  subroutine  FUNX.  It  was  taken 
directly 

C  from  Dan  Baumann's  program  D2ICC.  The  only  changes  made  were  to 
C  incorporate  the  CONSTANT  common  block,  to  change  the  control 
surface 

C  deflections  from  parameters  to  states  (i.e.,  DSTB0sU(13)  rather 
than 

C  DSTBDaPARd) ) ,  to  remove  PAR  from  the  argument  list,  to  use  the 
1988 

C  data  base  control  surface  deflection  neunes  (OSTBD,  DAILD,  DTALD, 
and 

C  DRUDD) ,  and  to  use  the  differential  tail  deflection  state  rather 
than 

C  .3*DAILD. 

C==ssssassssass=sa«s«sB»z«asssas«sssssaBSKsssssssx«3ss33SSS3sssss 


IMPLICIT  DOUBLE  PRECISION  (A-H,0-2) 

DIMENSION  U(NDIM) 

COMMON  /ACDATA/  SWING, CWING, SREF,RHO,FMASS, THRUST 
COMMON  /SEIZE/  CX,CY,C2,CLM,CMM,CNM 

C  THE  PRIMARY  SOURCE  OF  'THESE  COEFFICIENT  EQUATIONS  IS  SUBROUTINE 

C  AROlO  FROM  MCAIR  CODE  USED  IN  THE  F15  BASELINE  SIMULATOR. 

C  MOST  OF  THE  COEFFICIENTS  USED  IN  THE  EQUATIONS  WERE  COMPUTED 
C  USING  SAS  WITH  RAW  DATA  FROM  THE  F15  SIMUIATOR  DATA  TABLES. 

C  THIS  SUBROUTINE  IS  CALLED  BY  THE  DRIVER  PROGRAM  FOR  THE  AUTO 
C  SOFTWARE.  IT  MERELY  TAKES  INPUTS  ON  THE  A/C  STATE  AND  CONTROL 

C  SURFACE  POSITIONS  AND  RETURNS  THE  APPROPRIATE  AERO  COEFFICIENTS 

C  CX,  CY,  CZ,  CL,  CM,  AND  CM. 

C  22  July  1990-  Due  to  the  fact  that  this  is  a  merger  of  a  later 

C  subroutine  to  an  earlier  version,  some  of  the  complexity  of 

C  the  later  version  has  to  be  taken  out.  In  this  case  it  C 
involves  restating  certian  constants. 

C  GRAVITATIONAL  CONSTANT 
G-32.174 

C  DEGREE  TO  RADIAN  CONVERSION 

DEGRAC-57. 2957795131 

C  PI 

PI.3.14J.5926536 

C  INPUTS  TO  THIS  SUBROUTINE 

C  AL  -  ANGLE  OF  ATTACK,  DEG 
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C  BETA 
C  DDA 
C  DELEOD 
C  DELESD 
C  DRUDD 


SIDESLIP  ANGLE,  DEG 

AILERON  DEFLECTION  ANGLE,  DEG 

DIFFERENTIAL  TAIL  DEFLECTION  ANGLE,  DEG 

SYMMETRICAL  TAIL  DEFLECTION  ANGLE,  DEG 

RUDDER  DEFLECTION,  POSITIVE  TRAILING  EDGE  LEFT,  DEG 


CP  -  ROLL  RATE,  RAD/ SEC 

C  Q  -  PITCH  RATE,  RAD/ SEC 

C  R  -  YAW  RATE.  RAD/SEC 

C  THRUST  -  TOTAL  ENGINE  THRUST,  POUNDS 
C  VTRFPS  -  TRUE  AIRSPEED,  FT/SEC 


C  INTERMEDIATE  VARIABLES  USED  IN  THIS  SUBROUTINE 


C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


ABET  -  ABSOLUTE  VALUE  OF  BETA,  DEG 

ARUD  -  ABSOLUTE  VALUE  OF  RUDDER  DEFLECTION,  DEG 

BWING  -  WING  SPAN,  FEET 

CA  -  COSINE  RAL  (RAL  IN  RADIANS) 

CD  -  COEFFICIENT  OF  DRAG 

CL  -  BASIC  LIFT  COEFFICIENT 

CWING  -  MEAN  AERODYNAMIC  CHORD.  FEET 

DAHD  -  DIFFERENTIAL  ELEVATOR  DEFLECTION,  DEG 

DAHLD  -  LEFT  AILERON  DEFLECTION,  DEG 

DAHRD  -  RIGHT  AILERON  DELFECTION,  DEG 

DELEDR  -  DIFFERENTIAL  TAIL  DEFLECTION  ANGLE,  RAD 

DELESR  -  SYMMETRIC  TAIL  DEFLECTION  ANGLE,  RAD 

QBARS  -  DYNAMIC  PRESSURE  TIMES  WING  REFERENCE  AREA,  LBF 

RABET  -  ABSOLUTE  VALUE  OF  BETA,  RADIANS 

RAL  -  ABSOLUTE  VALUE  OF  ALPHA,  RADIANS 

RARUD  -  ABSOLUTE  VALUE  OF  RUDDER,  RADIANS 

SA  ~  SINE  RAL  (RAL  IN  RADIANS) 


C  OUTPUTS  FROM  THIS  SUBROUTINE 


C  CX  -  BASIC  AXIAL  FORCE  COEFFICIENT,  BODY  AXIS,  <*■  FORWARD 


C  CY 
C  CZ 
C  CLM 
C 

C  CMM 
C  CNM 


-  BASIC  SIDE  FORCE  COEFFICIENT,  BODY  AXIS,  i-  RIGHT 

-  BASIC  NORMAL  FORCE  COEFFICIENT,  BODY  AXIS,  ♦  DOWN 

-  BASIC  ROLLING  MOMENT  COEFFICIENT,  BODY  AXIS.  +  R  WING 

DOWN 

-  BASIC  PITCHING  MOMENT  COEFFICIENT,  BODY  AXIS,  NOSE 
UP 

-  BASIC  YAWING  MOMENT  COEFFICIENT,  BODY  AXIS,  NOSE 
RIGHT 


C  ANGLES  USED  IN  CALCULATING  CL,  CLLDB,  ....  ARE  IN  RADIANS. 
C  THIS  IS  BECAUSE  RADIANS  WERE  USED  IN  THE  CURVE  FITTING  PROGRAM 
C  TO  OBTAIN  THE  COEFFICIENTS  OF  THE  ALPHA,  BETA,  ....  TERMS  IN 
C  THE  FOLLOWING  EQUATIONS. 

C  MOMENT  REFERENCE  CENTER  WAS  SET  IN  AROlO  PROGRAM  AS: 
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C  DATA  CMCGR  / .2S65/ ,  CNCGR  /.2565/ 

C  THE  AERO  STABILITY  DATA  WAS  TAKEN  REFERENCED  TO  THESE  CG 
C  LOCATIONS.  THE  MOMENTS  OF  INERTIA  AND  OTHER  AIRCRAFT  DATA 

C  ARE  FOR  A  CLEAN  CONFIGURATION  TEST  AIRCRAFT  WITH  A  CG  AT 

C  THE  SAME  CG.  AS  A  RESULT,  THERE  IS  NO  'CG  OFFSET'  TO  BE 
C  COMPUTED. 

AL«U(1) 

BETA«U(2) 

P*U(3) 

Q=U(4) 

RaU(5) 

THETA»U(6) 

PHI«U(7) 

VTRFPS=U(8) *1000. 

DSTBD*U ( 9 ) 

DAILDsU(lO) 

DTALDsUdl) 

DRUDD*U{12) 

QBARS* . 5*RHO*VTRFPS*VTRFPS*SREF 
C02 VaCWING/ ( 2 . * VTRFPS ) 

B02V»BWING/ ( 2 . *VTRFPS ) 

QSB=BWING*QBARS 
ARUDsABS (DRUDD) 

RARUD »ARUD / DEGRAD 

RALsAL/DEGRAD 

ABETaABS{BETA) 

RABET« ABET/ DEGRAD 
RBETA»BETA/DEGRAD 
DAILA»ABS(DAILD) 

DSTBRiiDSTBD/ DEGRAD 


NEW  SECTION  OF  CODE  -  1)  ALL  THE  AERODYNAMIC  COEFFICIENTS  IN 

THIS  VERSION  OF  THE  DRIVER  PROGRAM 
ARE  TAKEN  DIRECTLY  FROM  THE  1988 
F15  AEROBASE  (0.6  MACH,  20000  FEET) 

2)  THIS  SECTION  SUMMARIZES  THE 
AERODYNAMIC  COEFFICIENTS  AS  TO  WHAT 
THEY  ARE  AND  HOW  THEY  ARE  USED. 

THE  FIRST  ACCRONYM  IS  THE  JOVIAL 
NAME  OF  THE  AERODYNAMIC  COEFFICIENT 
(CFX1,ETC),  THE  SECOND  ACCRONYM  IS 
THE  FI 5  AEROBASE  CODE  OR  CTAB  NAME 
(ATAB15,  ETC) .  A  BRIEF  DEFINITION 
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OF  THE  AERODYNAMIC  COEFFICIENT  IS 
ALSO  PROVIDED. 


C 
C 
C 

C  3)  THERE  IS  ALSO  A  SECTION  THAT  C 

PROVIDES  A  TABLE  OF  CONVERSIONS  C 
BETWEEN  WHAT  THE  VARIABLE  IS  CALLED  C 
IN  THE  ORIGINAL  SECTION  OF  THIS  C 
PROGRAM  AND  ITS  NAME  IN  THE  1988 

C  F15  AEROBASE.  FOR  THE  SAKE  OF  C 

CONTINUITY  THE  ORIGINAL  PROGRAM  NAME  C 
IS  USED  AND  THE  1988  F15  AEROBASE  C 
NAME  IS  PROVIDED  AS  BOOK  KEEPING 
C  INFORMATION . 

C 

C 

. . ****** . . 

c 

c 

C  CFX  S  FORCE  IN  STABILITY  AXIS  X  DIRECTION  (CD  IN  BODY  AXIS) 

C  (FUNCTION  OF  CL  OR  CFZl) 

C  CFX  =  CFXl  +  CXRB  +  STORE  INCREMENTS  *  CXDSPD  +  DCXLG  +  DCD 
C 

C  CFXl  =  ATAB15  «  PERFORMANCE  DRAG  COEFFICIENT  -  CD 
C  CXRB  =  ATAB22  =  DELTA  CD  DUE  TO  CG  (=0.0) 

C  CXDSPD  =  ATAB27  =  DELTA  CD  DUE  TO  SPEEDBRAKE  (NORMALLY  =  0.0436) 
C  SET  TO  0  SINCE  THIS  STUDY  IS  CONCERNED 

C  WITH  HIGH  ANGLES 

C  OF  ATTACK  PHENOMENON  (>40  DEGREES)  AND  BECAUSE 

C  THE  SPEEDBRAKE  WILL  NOT  DEPLOY  AT  ANGLES  OF 

C  ATTACK  GREATER  THAN  15  DEGREES. 

C  DCXLG  «  ATAB19  =  DELTA  CD  DUE  TO  REYNOLD'S  NUMBER  (=-0.0005) 

C  DCD  =  BTAB03  =  DELTA  CD  DUE  TO  2-PLACE  CANOPY  (F15B)  (=0.0005) 

C  *******  fioTE  THAT  DCXLG  AND  DCD  CANCEL  EACH  OTHER  ******* 

C 
C 

C**  ******** . . . ***••**••****• . 

c 

c 

c 

C  CFY  «  FORCE  IN  BODY  AXIS  Y  DIRECTION 

C  CFY  *  CFY1*EPA02  CYT)AD*DAILD  (CYDRD*DRUDD*DRFLX5)  *EPA43 
C  +[CYDTD*DTFLX5  +  DTFLX6] *DTALD  +  CFYP*PB  +  CFYR*RB 

C  •►CYRB  ♦  STORE  INCREMENTS  +  DCYB*BETA 

C 

C  CFYl  «  ATAB16  =  BASIC  SIDE  FORCE  COEFFICIENT  -  CY(BETA) 

C  EPA02  =  ATAB21  *  BETA  MULTIPLIER  TABLE 

C  CYDAD  =  ATAB75  =  SIDE  FORCE  COEFFICIENT  DUE  TO  AILERON 
C  DEFLECTION 

C  DAILD  =  AILERON  DEFLECTION  (DEG) 

C  CYDRD  =  ATAB69  =  SIDE  FORCE  COEFFICIENT  DUE  TO  RUDDER  DEFLECTION 
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C  DRUDD  =  RUDDER  DEFLECTION  (DEG) 

C  DRFLX5  =  ATAB88  =  FLEX  MULTIPLIER  ON  CYDRD  (=0.89) 

C  EPA43  =  ATAB30  *  MULTIPLIER  ON  CNDR,  CLDR,  CYDR  DUE  TO 
C  SPEEDBRAKE 

C  (=1.0) 

C  CYDTD  =  ATAB72  =  SIDE  FORCE  COEFFICIENT  DUE  TO  DIFFERENTIAL  TAIL 

C  DEFLECTION  -  CYDDT 

C  DTFLX5  »  ATABIO  «  FLEX  MULTIPLIER  ON  CYDTD  (=0.975) 

C  DTFLX6  =  ATAB77  =  FLEX  INCREMENT  TO  CYDTD  (=0.0) 

=  DIFFERENTIAL  TAIL  DEFLECTION  (DEG)  WHICH  IS 
DIRECTLY  PROPORTIONAL  TO  AILERON  DEFLECTION 
AND  IS  PRIMARILY  USED  TO  ASSIST  IN  ROLLING  THE 
F-15B  (DTALD=0.3*DAILD) 

*  SIDE  FORCE  COEFFICIENT  DUE  TO  ROLL  RATE  (CYP) 

=  (PEOBB*SPAN)/(2*VILWF) 

PEOBB  =  ROLL  RATE  IN  RAD/ SEC  =  P 
SPAN  »  WING  SPAN  =42.8  FEET  =  BWING 
VILWF  =  VELOCITY  IN  FT/SEC  =  VTRFPS 
=  SIDE  FORCE  COEFFICIENT  DUE  TO  YAW  RATE  (CYR) 

=  (REOBB*SPAN)/(2*VILWF) 

REOBB  =  YAW  RATE  IN  RAD/ SEC  =  R 
=  ASSYMETRIC  CY  AT  HIGH  ALPHA  (ANGLE  OF  ATTACK) 


C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 


DTALD 


CFYP  =  ATAB13 
PB 


CFYR 

RB 


ATAB07 


CYRB  =  ATAB93 


C  DCYB  *  0,0  THERE  IS  NO  INCREMENT  DELTA  CYB  (SIDE  C 

FORCE)  DUE  TO  A  2-PLACE  CANOPY  ON  THE  F15B.  C 
THIS  IS  BECAUSE  THE  SAME  CANOPY  IS  USED  ON  BOTH  C 
THE  BASELINE  F15A  AND  THE  F15B.  THE  SIDEFORCE  C 
IS  THE  SAME  FOR  BOTH  '/ERSIONS  OF  THE  F15  AND  C 
ALREADY  INCLUDED  IN  THE  BASIC  SIDE  FORCE  (CFYl) . 

C  THE  TWO  PLACE  CANOPY  IS  MOUNTED  DIFFERENTLY  C 

HOWEVER,  SO  THERE  IS  A  DIFFERENCE  IN  YAWING  AND  C 
ROLLING  MOMENT.  (SEE  DCNB  AND  DCLB) 


CFZ  =  FORCE  IN  STABILITY  AXIS  Z  DIRECTION  (CL  IN  BODY  AXIS) 
CFZ  =  CFZl  +  CZDSPD  *  STORE  INCREMENTS  ♦  DCL*BETA 


CFZl  =  ATAB17  =  BASIC  LIFT  COEFFICIENT  -  CL 
CZDSPD  =  ATAB2(5  =  DELTA  CL  DUE  TO  SPEEDBRAKE 

SET  TO  0  DUE  TO  THE  REASONS  GIVEN  ABOVE  IN 
CXDSPD 

DCL  *  BTABCl  =  DELTA  CL  DUE  TO  2 -PLACE  CANOPY  (F15B)  (=0.0) 
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CML  =  TOTAL  ROLLING  MOMENT  COEFFICIENT  IN  BODY  AXIS 
CML  =  CML1*EPA02  +  CLDAD*DA1LD  +  [CLDRD*DRUDD*DRFLX1 ] *EPA43  + 
[CLDTD*DTFLX1  ♦  DTFLX2 ) 'DTALD  +  CMLP*PB  +  CMLR*RB  + 
STORE  INCREMENTS  ♦  CLDSPD  DCLB*BETA 


CMLl  a  ATABOl  »  BASIC  ROLLING  MOMENT  COEFFICIENT  -  CL (BETA) 
EPA02  =  ATAB21  *  BETA  MULTIPLIER  TABLE 

CLDAD  a  ATAB73  =  ROLL  MOMENT  COEFFICIENT  DUE  TO  AILERON  C 
DEFLECTION 
- (CLDA) 

DAILD  =  AILERON  DEFLECTION  (DEG) 

CLDRD  a  ATAB67  a  ROLLING  MOMENT  COEFFICIENT  DUE  TO  RUDDER 
DEFLECTION  - (CLD) 

DRUDD  =  RUDDER  DEFLECTION  (DEG) 

DRFLXl  a  ATAB80  a  FLEX  MULTIPLIER  ON  CLDRD  (=0.85) 

EPA43  a  ATAB30  =  MULTIPLIER  ON  CNBR,  CLDR,  CYDR  DUE  TO 
SPEEOBRAKE 
(=1.0) 

CLDTD  a  ATAB70  =  ROLL  MOMENT  COEFFICIENT  DUE  TO  DIFFERENTIAL 
TAIL  DEFLECTION  -  CLDD 

DTFLXl  a  ATAB04  =  FLEX  MULTIPLIER  ON  CLDTD  (=0.975) 

DTFLX2  a  ATAB84  =  FLEX  INCREMENT  TO  CLDTD  (=0.0) 

DTALD  a  DIFFERENTIAL  TAIL  DEFLECTION  (DEG)  WHICH  IS 

DIRECTLY  PROPORTIONAL  TO  AILERON  DEFLECTION 
AND  IS  PRIMARILY  USED  TO  ASSIST  IN  ROLLING  THE 
F-15B 

(DTALD  a  0.3 ’DAILD) 

CMLP  a  ATAB02  =  ROLL  DAMPING  DERIVATIVE  -CLP 

PB  a  (PEOBB’SPAN) / (2*VILWF) 

PEOBB  a  roll  rate  IN  RAD/ SEC  =  P 
SPAN  a  WING  SPAN  =42.8  FEET  =  SWING 
VILWF  a  VELCXriTY  IN  FT/ SEC  =  VTRFPS 
CMLR  a  ATABll  =  ROLLING  MOMENT  COEFFICIENT  DUE  TO  YAW  RATE  -  CLR 

RB  «  (REOBB*SPAN)/(2*VILWF) 

REOBB  a  YAW  RATE  IN  RAD/ SEC  =  R 
CLDSPD  a  ATAB29  =  DELTA  CL  DUE  TO  SPEEDBRAKE 

SET  TO  0  DUE  TO  THE  REASONS  GIVEN  ABOVE  IN 
CXDSPD 

DCLB  a  BTAB04  =  INCREMENT  DELTA  CLB  (ROLLING  MOMENT)  DUE  TO  C 
2 -PLACE 

C  CANOPY  FROM  PSWT  499 
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CMM  *  TOTAL  PITChING  MOTIEIT  COEFFICIENT  IN  STABILITY  AXIS 
(BODY  AXIS  -  AS  WELL) 

CMM  =  CMMl  CMMQ*QB  ♦  STORE  INCREMENTS  +  CMDSPD  +  DCM 

CMMl  =  ATAB03  «  BASIC  PITCHING  MOMENT  COEFrTCIENT  -  CM 
CM-IQ  =  ATAB05  »  PITCH  DAMPING  DERIVATIVE  -  CMQ 
QB  *  (QEOBB*MAC)/(2*VILWF) 

QEOBB  «  PITCH  RATE  IN  RAD/ SEC  *  Q 

MAC  ■  MEAN  AERODYNAMIC  CHORD  *  15.94  FEET  *  C 

WING 

VILWF  =  VELOCITY  IN  FT/ SEC  =  VTRFPS 
CMDSPD  =  ATAB25  *  DELTA  CM  DUE  TO  SPEEDBRAKE 

SET  TO  0  DUE  THE  REASONS  GIVEN  ABOVE  IN  CXDSPD 
iX:M  =  BTAB02  =  DELTA  CM  DUE  TO  2-PLACE  CANOPY  (F15B)  (=0.0) 


CMN  =  TOTAL  YAWING  MOMENT  COEFFICIENT  IN  BODY  AXIS 

CMN  =  CMN1*EPA02  CNDAD*DAILD  ♦  {CNDRD*DRUDD*DRFLX3  ]  *EPA43 

+  (CNDTD*DTLX3  +  DTFLX4 )  *DTALD  +  CMNP*PB  >  CMNR*RB  +  CN°'- 
+DCNB2*EPA36  STORE  INCREMENTS  ♦  CNDSPD  ♦  DCNB*BETA 


CMNl  =  ATi'Bl?  »  BASIC  YAWING  MOMENT  COEFFICIENT  -  CN  (BETA) 

EPA02  =  ATAB21  =  BETA  MULTIPLIER  TABLE 

CNDAD  =  ATAB74  *  YAW  MOMENT  COEFFICIENT  DUE  TO  AILERON 

DEFLECTION  -CNDA 

DAILD  =  »  AILERON  DEFLECTION  (DEG) 

CNDKD  =  ATAB68  =  YAWING  MOMENT  COEFFICIENT  DUE  TO  RUDDER 

DEFLECTION  -CNDR 

DRUDD  =  RTODER  DEFLECTION  (DEG) 

DRFLX3  s  ATAB85  =  FJ.EX  MULTIPLIER  ON  CNDRD 
EPA43  =  ATAB30  =  MULTIPLIER  ON  CNDR,  CLDR,  CYDR  DUE  TO  C 
SPEEDBRAKE 

CNDTD  =  ATAB71  s  YAWING  MOMENl'  COEFFICIENT  DUE  TO  DIFFERENTIAL 

TAIL 

DEFLECTION  -  CNDDT 

DTFLX3  »  ATAB08  «  FLEX  MULTIPLIER  ON  CNDTD 
DTFLX4  »  ATAB09  »  FLEX  INCREMENT  ON  CNDTD  (=0.0) 

DTALD  «  «  DIFFERENTIAL  TAIL  DEFLECTION  (DEG)  WHICH  IS 

DIRECTLY  PROPORTIONAL  TO  AILERON  DEFLECTION 
AND  IS  PRIMARILY  USED  TO  ASSIST  IN  ROLLING 
THE  F-15B  (DTALD  =  0.3*DAILD) 

CMNP  =  ATAB06  =  YAWING  MOMENT  COEFFICIENT  DUE  TO  ROLL  RATE  -  CNP 
PB  *  (PEOBB*SPAN) / (2*VILWF) 

PEOBB=ROLL  RATE  IN  RAD /SEC  =  P 
SPAN  »  WING  SPAN  =  42.8  FT  *  BWING 
VILWF  «  VELOCITY  IN  FT/SEC  =  VTRFPS 
CMNR  *  ATAB14  =  YAV.’  DAMPING  DERIVATIVE  -  CNR 
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C  RB  *  (RE0BB*SPAN)/{2*VILWF) 

C  REOBB  «  YAW  RATE  IN  RAD/SEC  *  R 

C  CNRB  *  ATAB86  *  ASSYMETRIC  CN  AT  HIGH  ALPHA 

C  DCNB2  s  ATAB44  =  DELTA  CNB  WITH  STABILATOR  EFFECT  -  DELCNB  C 
(»0.0) 

EPA36  »  ATAB94  *  MULTIPLIER  ON  DCNB2  (=BETA) 

CNDSPD  *  ATAB28  =  DELTA  CN  DOE  TO  SPEEDBRAKE 

SET  TO  0  DUE  TO  THE  REASONS  GIVEN  ABOVE  IN 
CXDSPD 

DCNB  *  BTAB05  *  INCREMENT  DELTA  CNB  (YAWING  MOMENT)  DUE  TO 
2 -PLACE  CANOPY  (FI SB) 


MISCELLANEOUS  COEFFICIENTS  AND  NAME  CONVERSION  TABLE 


1988  F15 

AEROBASE  NAME 

ORIGINAL 
PROGRAM  NAME 

DEFINITION 

AL77D 

AL 

ANGLE  OF  ATTACK 
(DEG) 

BE77D 

BETA 

SIDESLIP  ANGLE 
(DEG) 

BE77D 

RBETA 

SIDESLIP  ANGLE 
(RAD) 

B077D 

ABET 

ABSOLUTE  VALUE  OF 
SIDESLIP  ANGLE 
(DEG) 

DAI  LA 

DAI  LA 

ABSOLUTE  VALUE  OF 
r  ’LERON  DEFLEC- 
i  (DEG) 

DAILD 

DDA 

70N  DEFLEC- 

..  .  :deg) 

DRUADS 

ARUD 

ABSOLUTE  VALUE  OF 
RUDDER  DEFLEC-  . 
TION  (DEG) 

DRUABS 

RARUD 

ABSOLUTE  VALUE  OF 
RUDDER  DEFLEC¬ 
TION  (RAD) 

DRUDD 

DRUDD 

RUDDER  DEFLECTION 
(DEG) 

DSTBD 

DELESD(R) 

AVERAGE 

STABILATOR 

DEFLECTION 

DEG  (RAD) 

DTALD 

DELEDDCR) 

DIFFERENTIAL  TAIL 
DEFLECTION 
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C  DEG  (RAD) 

C 

C 

c 

c 

c 

c 

PB« (P*BWING) / (2*VTRFPS) 

QB« (Q-CWING) / (2*VTRFPS) 

RB« (R*BWING) / (2*VTRFPS) 

C 

C  THE  F-15B  AERO  DATA  TABLES  DO  NOT  CONTAIN  STABILITY  COEFFICIENT 
C  DATA  FOR  BETA  AND  RUDDER  DEFLECTION  ,DRUDD,  LESS  THAN  0 

C  DEGREES.  THE  ABSOLUTE  VALUE  OF  BETA,  ABET,  AND  THE  ABSOLUTE 

C  VALUE  OF  RUDDER  DEFLECTION,  ARUDD,  ARE  USED  IN  THE  FOLLOWING 
C  EQUATIONS.  IN  RADIANS  THESE  PARAMETERS  ARE  RABET  AND  RARUD, 

C  RESPECTIVELY.  IN  SCSffi  CASES  THE  COEFFICIENT  IS  MULTIPLIED  BY 
C  A  -1  FOR  PARAMETER  VALUES  LESS  THAN  ZERO. 

C 

C  EPA02  IS  A  MULTIPLIER  THAT  ADJUSTS  THE  PARTICULAR  COEFFICIENT 
C  IT  IS  WORKING  ON  (CFY1,CML1,CMN1)  BY  CHANGING  THAT  PARTICULAR 
C  COEFFICIENTS  SIGN  (POSITIVE  OR  NEGATIVE)  DEPENDENT  ON  THE  SIGN 
C  OF  THE  SIDESLIP  ANGLE  (BETA)  .  IF  BETA  IS  NEGATIVE  THEN 
C  EPA02--1.0.  IF  BETA  IS  POSITIVE  THEN  EPA02-1.0.  SINCE  THIS 
C  FUNCTION  IS  DISCONTINUOUS  AT  THE  ORIGIN  A  CUBIC  SPLINE  HAS 
C  BEEN  EMPLOYED  TO  REPRESENT  THIS  FUNCTION  IN  ORDER  THAT 
C  AUTO  CAN  RUN. 

C 

C  Revised  15  Jun  90  -  Combined  three  IF  statements  into 

IF-THEN-ELSE. 

IF(BETA.LE.-1.0)THEN 
EPA02S>-1.00 
ELSEIF ( BETA . GE . 1 . 0 ) THEN 
EPA02S-1.00 
ELSE 

EPA02S— 1.00+(1.5*((BETAfl.0)**2)  >-(0.5*(  (BETA-1.0)**3) ,' 
ENDIF 
C 

C  Revised  15  Jun  90  -  Combined  three  IF  statements  into 

IF-THEN-ELSE. 

IF(BETA.LE.-5.0)THEN 
EPA02L— 1.00 
ELSEIF ( BETA . GE . 5 . 0 ) THEN 
EPA02L-1.00 
ELSE 

EPA02L--1. 00+ (0.06* ( (BETA^5.0)**2) )- (0.004* ( (BETA^5.0)**3) ) 

ENDIF 

C 

C 

C 

C 
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FZl—0 .00369376+  (3 .78028702*RAL)  ♦  (0 .6921459*RAL*RAL)  -  (5 . 0005867 

•(RAL**3) )+ (1.94478199* (RAL**4) )+(0.40781955*DSTBR)+ (0.10114579 
♦*(DSTBR*DST3R) ) 

CFZ-CF7.1 


CL=CFZl/57. 29578 

THIS  CONVERSION  OF  CFZl  TO  CL  IS  AN  ARTIFACT  FROM  THE 
CURVE  FITTING  PROCESS  WHERE  ALL  THE  INDEPENDENT  VARIABLES 
WERE  ANGLES  THAT  WERE  CONVERTED  FROM  DEGREES  TO  RADIANS. 
IT  JUST  SO  HAPPENED  THAT  FOR  CFXl  ONE  OF  THE  VARIABLES 
WAS  NOT  AN  ANGLE  BUT  A  DIMENSIONLESS  COEFFICIENT. 


CFXl-0 . 01806821+ (0 .01556573*CL) + (498 . 96208868*CL*CL) 
+-(14451. 56518396*(CL**3) >+(2132344. 6184755*(CL**4)) 

TRANSITIONING  FROM  LOW  AOA  DRAG  TABLE  TO  HIGH  AOA  DRAG  TABLE 

CFX2-0 .0267297- (0 . 10646919*RAL) + (5 . 39836337*RAL*RAL) 

+- (5. 0086893* (RAL**3) >+(1.34148193* (RAL**4) )+ 

+  (0.20978902*DSTBR>  +  (0. 30604211* (DSTBR**2>  >+0.09833317 
C 

A1»20.0/DEGRAO 
A2»3 0,0 /DEGRAD 
A12sA1+A2 

BA-2 .0/(-Al**3+3.*Al*A2*(Al-A2)+A2**3> 

BB--3.*BA*(Al+A2>/2. 

BC-3 .*BA*A1*A2 

BD«BA*A2 * *2  * ( A2 -3 . *A1 >  '2 . 

F1»BA*RAL**3+BB*RAL**2«-BC*RAL+BD 
F2»-BA*RAL**3+ (3 .*A12*BA+BB>  *RAL**2- 
+  (BC+2.*A12*BB+3.*A12**2*BA)*RAL+ 

+  B0+A12*BC+A12**2*BB+A12**3*BA 

C 

IF  (RAL  .LT.  Al>  THEN 
C 

CFX-CFXl 


130 


O  O  *000 


c 

c 

c 

n 

c 


ELSEIF  (RAL  .GT.  A2)  THEN 
CFX-CFX2 
ELSE 

CFX«CFX1*F1+CFX2*F2 

ENDZF 


DTFLX5a0.975 

DRFLX5>0.89 

C 

CFYla-O  .  05060386- (0 . 12342073*RAL)  (1 . 04501136*RAL*RAL) 

(0. 17239516*  (RAL**3) )- (2. 90979277*  (RAL**4)  ) 

(3. 06782935* (RAL**5) )- (0.88422116* (RAL**6) ) 

♦-( 0.065 78812 *R>X*RABET)-(0. 7152 1988*RABET) -(0.00000475273 
+  *  (RABET**2 )  )  -  (0 .04856168*RAL*DSTBR)  -  (0 . 05943607*RABET*DSTBR)  ♦ 
+(0.02018534*DSTBR) 

C 

IF  (RAL  .LT.  .52359998)  THEN 
C 

CFYP«0.014606188*(2.52405055*RAL)- (5.02687473* (RAL**2) ) 

+  -(106.43222962*(RAL**3) )♦ (256.80215423* (RAL**4 ) ) 

++ (1256. 39636248* (RAL**5) ) 

(3887 . 92878173* (RAL**6) )- (2863 .16083460* (RAL**7) ) ♦ 

♦ (17382. 72226362* (RAL**8) )- (13731. 65408408* (RAL**9) ) 

ENDIF 

C 

IF  ((RAL  .GE.  .52359998)  .AND.  (RAL  .LE.  .610865))  THEN 
C 

CFYP-0 . 00236511+ (0 . 52044678* (RAL-0 . 52359998 ))-( 12 . 8597002*  (RAL- 
+0.52359998) **2 )+ (75 .46138* (RAL-0 . 52359998) **3 ) 

ENDIF 

C 

IF  (RAL  .GT.  0.610865)  THEN 
C 

CFYP.0.0 

ENDIF 

C 

IF  (RAL  .LT.  -0.06981)  THEN 
C 

CFYR.0.35 

ENDIF 


C 
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IF  ((RAL  .GE.  -0.06981)  .AND.  (RAL  .LT.  0.0))  THEN 
C 

CFYR-O .34999999^ (35.4012413* (RAL+0 .06981) **2 ) - (493.33441162* 
♦ (RAL+0. 06981 ) **3) 

ENDIF 

C 

IF  ((RAL  .GE.  0.0)  .AND.  (RAL  .LE.  0.523599))  THEN 
C 

CFYRsO. 35468605- (2. 26998141*RAL)+ (51. 82178387*RAL*RAL) 

+- (718.55069823* (RAL**3) ) 

++(4570.00492172* (RAL**4) )- (14471 . 88028351* (RAL**5) ) + 
+(22026.58930662* (RAL**6) )- (12795. 99029404* (RAL**7) ) 

ENDIF 

C 

IF  ((RAL  .GT.  0.523599)  .AND.  (RAL  .LE.  0.61087))  THEN 
C 

CFYR=0 . 00193787+ (1 .78332496* (RAL-0 . 52359903 ) ) - (41 . 63198853*(RAL- 
+0. 52359903 )**2)+ (239. 97909546* (RAL-0. 52359903 )**3) 

ENDIF 

C 

IF  (RAL  .GT.  0.61087)  THEN 
C 

CFYR»0.0 

ENDIF 

C 

IF  (RAL  .LT.  0.55851)  THEN 
C 

CYDAD— 0.00020812+ (0.00062122*RAL)+(0.00260729*RAL*RAL) 

++ (0.00745739* (RAL**3) )- (0.0365611* (RAL**4) ) 

+- (0.04532683* (RAL**5) )+ (0.20674845* (RAL**6) ' 

+- (0.13264434* (RAL**7) )-(0. 00193383* (RAL**8) ) 

ENDIF 

C 

IF  ((RAL  .GE.  0.55851)  .AND.  (RAL  .LT.  0.61087))  THEN 
C 

CYDAD-0. 00023894+ (0.00195121* (RAL-0. 55851001) )+ (0.02459273 
+* (RAL-0. 55851001 ) **2)- (0. 1202244* ( (RAL-0. 55851001 ) **3) ) 
ENDIF 
C 

IF  (RAL  .GE.  0.61087)  THEN 
C 

CYDAD.O. 27681285- (2. 02305395*RAL)+ (6. 01180715*RAL*RAL) 

+- (9. 24292188* (RAL**3) ) +(7.59857819* (RAL**4) ) 

+- (2. 8565527* (RAL**5))+(0. 25460503 •(RAL**7) ) 

+-(0.01819815* (RAL* *9) ) 

ENDIF 


IF  (RAL  .LE.  0.0)  THEN 
EPA43«1.0 
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C  ENDIF 

C  IP  (RAL  .GT.  0.0  AND  .LE.  0.6283185)  THEN 

0.6283185  RADIANS  •  36  DEGREES 

EPA43-0.9564809>(4.13369452*RAL} -(18.31288396*RAL*RAL)1> 
•^{19  .55114S6*  (RAL**3) ) - (1 .09295946*RAL*DSPBD) f (0 . 17441033* 
4-DSPBD*0SPBD) 

ENDIF 

IF  (RAL  .GT.  0.6283185)  THEN 
EPA43-1.0 
ENDIF 


C  •  NOTE  -  THE  PARAMETER  EPA43  IS  A  MULTIPLIER  ON  RUDDER 

C  *  EFFECTIVENESS  DUE  TO  SPEEDBRAKE.  THIS  TABLE  IS  ALSO 

C  *  LIMITED  TO  36  DEG  AOA.  HOWEVER,  THERE  IS  NO  AERODY 

# 

C  *  NAMIC  EFFECT  FOR  ANGLES  OF  ATTACK  LESS  THAN  16  DEG, 

C  *  AND  THE  SPEEDBRAKE  IS  AUTOMATICALLY  RETRACTED  AT  AOA 

* 

C  *  GREATER  THAN  15  DEG.  THEREFORE,  THIS  TABLE  SHOULD 

* 

C  *  NOT  BE  NECESSARY  FOR  THE  ORDINARY  OPERATION  OF  THE 

C  *  AIRCRAFT 

C 


C 

CYDRD.O. 00310199+ {0.00119963*RAL)+(0. 02806933 *RAL*RAL) 

+- (0. 12408447* (RAL**3) )- (0.12032121* (RAL**4) ) 

++ (0. 79150279* (RAL**5) )- (0.86544347* (RAL**6) ) 

♦♦ (0 . 27845115* /RAL**7) ) + (0 . 00122999*RAL*RARUD) + (0 . 00145943 

♦  *RARUD)  -(0.01211427*RARUD*RARUD)  + (0.00977937*  (RARUD**3)  ) 

C 

CYDTD.-0 . 00157745- (0 . 0020881*RAL) +  ( 0 . 00557239*RAL*RAL) 

♦- (0.00139886* (RAL**3) )+ (0.04956247* (RAL**4) ) 

♦  - (0. 0135353  * (RAL**5) )- (0.11552397* (RAL* *6) ) 

♦+ (0. 11443452* (RAL**7) )- (0.03072189* (RAL**8) )- (0.01061113* 

♦  ( RAL  *  *  3 )  *  DSTBR )  -  ( 0 . 0 0 0 1 0 5 2  9  *  RAL  *  RAL  *  DSTBR  *  DSTBR ) 

♦- (0 . 00572463 *RAL*DSTBR*DSTBR) 

♦+(0.01885361*RAL*RAL*DSTBR)-(0.01412258*RAL*(DSTBR**3) ) 
+-(0.00081776*DSTBR)+(0.00404354*(DSTBR**2) ) - 
+(0.00212189*(DSTBR**3) )+(0.00655063*(DSTBR**4) ) 
++(0.03341584*(DSTBR**5) ) 

C 

RALYl-0. 6108652 
RAL Y2> 9 0.0 /DEGRAD 
RBETY1--0. 0872665 
RBETY2.0. 1745329 
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c 

AY>0.164 
ASTARY-O. 95993 
BSTARY«0. 087266 
C 

2ETAY«(2  .ODO*ASTARY-(RALY1^RALY2)  )  /  (RALY2-RALY1 ) 

BTAY.  (2 .0D0*BSTARY-  (RBETY1-**RBETY2 ) )  /  (RBETY2-RBETY1) 

C 

X>  (2 .0D0*RAL- (RALYl-»-RALY2) )  /  (RALY2-RALY1) 

Y.  (2 .0D0*RBETA-  (RBETY1+RBETY2) )  /  (RBETY2-EBETY1) 

C 

rY-( (5.0D0*(2ETAY**2) )-(4.0DO*2ETAY*X)-1.0DO)*( ( (X**2) -l.ODO) 
■••**2)*(1.0DO/(  ((2ETAY**2)-1.0D0)**3) ) 

C 

GY*( {5.0DO*(ETAY**2) ) -(4.0D0*ETAY*Y)-1.0D0) *( ( ( Y«*2 ) -1 . ODC ) **2 ) 
♦•(1.0D0/( ( (ETAY**2) -1.000) **3) ) 

C 

CYRB«AY*FY*GY 

C 

IF  (RAL  .LT.  0.6108652)  THEN 
C 

CYRB-0.0 
GOTO  500 
ENDIF 

C 

IF  ((RBETA  .LT.  -0.0872665)  .OR.  (RBBTA  .GT.  0.1745329))  THEN 
C 

CYRB-0 . 0 
GOTO  500 
ENDIF 
C 

500  CFY.  (CFY1*EPA02L)  ♦  (CYDAD*DAILD)  ♦  (CYDRD*DRUDD*DRFLX5*EPA43 ) 

♦  {  (CYDTD*DTFLX5 )  •DTALD)  ♦  (CFYP’PB)  ♦  (CFYR'RB) 

♦+CYRB 


DTFLXl.0.975 

DRFLXl.0.85 

C 

CMLl— 0.00238235- (0.04616235*RAL)^(0.10553168*RAL*RAL) 
♦♦(0. 10541585* (RAL**3))- (0.40254765* (RAL**4) ) 
♦♦(0.32530491* (RAL**5) )- (0.08496121* (RAL**6) ) 

♦♦ (0 . 00112288* (RAL**7) ) - (0.05940477 *RABET*RAL) - 
♦  (0 .07356236*RABET)  -  (0 .00550119*RABET*RABET) ♦ (0.00326191 
♦*(RABET**3) ) 

C 
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IF  (RAL  .LT.  0.29671)  THEN 
C 

CMLP—0 . 24963201- (0 . 03106297*RAL)  K  (0 . 12430631*RAL*RAL) 

(8. 95274618* (RAL**3 ))-*•  (100. 33109929*  (RAL*'*4) ) 

♦♦ (275.70069578* (RAL**5))- (1178.83425699* (RAL**6) ) 
♦-(2102.66811522*(RAL**7) )+(2274.89785551*(RAL**8) ) 

ENDZF 

C 

IF  ((RAL  .GE.  0.29671)  .AND.  (RAL  .LT.  0.34907))  THEN 
C 

CMLP—0. 1635261- (3. 77847099* (RAL-0. 29671001)  )♦  (147.47639465 
♦*(RAL-0. 29671001) **2)- (1295. 94799805* (RAL-0. 29671001) **3) 
ENDIF 
C 

IF  (RAL  .GE.  0.34907)  THEN 
C 

C.<Ll't  -1 .  .m20291>  (7 .06112181*PAL)  -  (13 . 57010422 •RAL*RAL) 
♦+(11.21323950*(RAL**3) ) 

-(4.26789i25*(R*L**4)  )>(0.6237381*(RAL**5)  ) 

ENOI? 

C 

"F  (EJ»1  .LT.  0.7854)  THEN 
C 

CMLR^O  .  03':i5391+ (0 . 59296381*RAL)  >  (2 . 27456302*RAL*RAL) 
••■-(3.8097803*(RAL**3)) 

♦  -  (45  .e316'i842*  (RAL**4)  /  ♦  (55.31669213*  (RAL**5)  ) 

♦ (194. 29237485* (RAL**S)) -(393. 22969953* (RAL**7) )♦ (192. 20860739* 

♦  (R*JU**8)  ) 

ENDIF 

C 

IF  ((RAL  .GE.  0  7854)  .AND.  (RAL  .LE.  0.87266))  THEN 
C 

CMLR«0. 0925579071- (0. 6000000238* (RAL-0. 7853999734) ) 
♦♦(1.3515939713*( (RAL-0. 7853999734) **2) ) 

♦♦(29.0733299255*( (RAL-0. 7853999734) **3) ) 

ENDZF 

C 

IF  (RAL  .GT.  0.87266)  THEN 
C 

CMLR«-311. 126041-^  (1457. 23391042*RAL)- (2680. 19461944*RAL*RAL)-^ 

♦ (2361. 44914738* (RAL**3))- (893. 83567263* (RAL**4) )+ (68. 23501924* 
♦(RAL**6) )- (1. 72572994* (RAL**9)) 

ENDIF 

C 

CLDAD-0 . 00057626> (0 .00038479*RAL) - (0 . 00502091*RAL*RAL) 
♦♦(0.00161407*  (RAr,**3)  )♦  (0.02268829*  (RAL**4)  ) 

♦- (0. 03935269* (RAL**5) )♦ (0.02472827* (RAL**6) ) 

♦- (0. 00543345* (RAL**7) )^(0.0000007520348*DSTBR*RAL)^ 
♦(0.000000390773*DSTBR) 

C 
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CLDRD-0. 00013713- (0.00035439*RAL)-(0.00227912*RAL*RAL) 

♦♦  (0. 00742636*  {RAL**3 ))  (0.00991839*  (RAL**4)  ) 

(0. 04711846*  (RAL**5)  )♦  (0.046124*  (RAL**6)  ) 

(0.01379021* (RAL**7) )+(0.00003678685*RARUD*RAL)+ 

•►( 0 . 000010437S1*RARUD) - (0 . 00015866*RARUD*RARUD) ♦( 0 . 00016133 
+*(RARUD**3) > 

C 

CLDTDsO .00066663+ (0 .00074174*RAL) + (0 . 00285735*RAL*RAL) 

+- (0.02030692* (RAL**3) )- (0.00352997* (RAL**4) ) 

++ (0. 0997962* (RAL**5) >-(0.14591227* 

+(RAL**6) )+(0.08282004*(RAL**7) ) 

+  -(0.0168667* (RAL**8) >  +  (0.00306142* (RAL**3>  *DSTBR) 

+  -(0.00110266*RAL*RAL*(DSTBR**2>  >  +  (0 .00088031*RAL* 

+  (DSTBR**2 >  > - (0 . 00432594*RAL*RAL*DSTBR> - 
+(0.00720141*RAL*(DSTBR**3> > 

+- (0 . 00034325*DSTBR)  +  (0 .00033433 * (DSTBR**2 >  >  +  (0.00800183 
+  *(DSTBR**3>  > -(0.00555986* (DSTBR**4>  > -(0.01841172* (DSTBR**5> ) 
C 

IF  (RAL  .LT.  0.0>  THEN 
C 

DCLB*-0. 00006 
ENDIF 
C 

IF  ((RAL  .GE.  0.0)  .AND.  (RAL  .LE.  0.209434) >  THEN 
C 

DCLB=-0 .00006+ (0 .0041035078*RAL*RAL> - (0 . 0130618699* (RAL**3>  > 
ENDIF 
C 

IF  (RAL  .GT.  0.209434)  THEN 
C 

DCLBsO.O 

ENDIF 


CML=  (CML1*EPA02S)  +  (CLDAD*DAILD)  +  (CLDRD*DRUDD*DRFLX1  *EPA43 )  + 
+  (  (CLDTD*DTFLX1)  *DTALD)  +  (CMLP*PB)  +  (CMLR*RB)  +  (DCLB*BETA) 


CMMl.O . 00501496- (0 .08004901*RAL) -( 1 . 03486675*RAL*RAL) 
+-(0. 68580677* (RAL**3)) +(6.46858488* (RAL**4) ) 

+- (10. 15574108* (RAL* *5) )+ 

+(6.44350808* (RAL**6)) -(1.46175188* (RAL**7) ) 
++(0.24050902*RAL*DSTBR) 

+-(0.42629958*DSTBR)-(0.03337449*DSTBR*DSTBR) 

+- (0 , 53951733* (DSTBR**3) ) 

C 

IF  (RAL  .LE.  0.25307)  THEN 
C 
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CMMQ=-3.8386262+(13.54661297*RAL)+(402.53011559*RAL*RAL) 
-*■-(6660. 95327122*  (RAL**3)  )- (62257 .89908743*  (RAL**4)  ) 
■►♦(261526. 10242329*(RAL**5)  ) 

•►♦(2177190.33155227*(RAL**6))-(703575.13709062*(RAL**7))- 
(20725000. 34643054*  (RAL**8) )- (27829700. 53333649*  (RAL**9)  ) 
ENDIF 
C 

IF  ((RAL  .GT.  0. 25307)  .AND.  (RAL  .LT.  0.29671))  THEN 
C 

CMMQ.-8.4926528931-(2705.3000488281*(RAL-0.2530699968) ) 
♦♦(123801.5* (RAL-0. 2530699968) **2) 

♦-( 1414377* (RAL-0. 2530699968) **3) 

ENDIF 

C 

IF  (RAL  .GE.  .29671)  THEN 
C 

CMMQ«47.24676075-(709.60757056*RAL)^(3359.08807193*RAL*RAL)- 
♦ (7565. 32017266* (RAL**3) )♦ (8695. 1858091* (RAL**4) ) 

♦- (4891.77183313* (RAL**5) ) ♦(1061.55915089* (RAL**6) ) 

ENDIF 

C 

CMM*:CMM1+  (CMMQ*QB) 


DTFLX3=0 .975 
DRFLX3=0.89 


CMNlsO. 01441512^(0.02242944*RAL)-(0. 30472558* (RAL**2 ) ) 

♦♦(0. 14475549* (RAL**3) ) 

♦♦ 10. 93140112* (RAL**4) )- (1.52168677* (RAL**5) )+ 

♦ (0. 90743413* (RAL**6) )- (0.16510989* (RAL* *7) ) 
♦-(0.0461968*(RAL**8) ) 

♦♦ (0.01754292* (RAL**9))-(0.17553807*RAL*RABET)^ 

♦ (0 . 1541 5649 *RAL*RASET*DSTBR) 

♦+(0. 14829547* (RAL**2)*(RABET**2) ) 

♦- (0. 11605031* (RAL**2)*RABET*DSTBR) 
♦-(0.06290678*(RAL**2)*(DSTBR**2) ) 

+- (0 . 01404857* (RAL**2) * (DSTBR**2 ) ) 

♦♦(0.:?225609*RABET)-(0. 08567087* (RABET* *2) ) 

+♦ (0. 01184674* (RABET**3) ) 

♦- (0 . 00519152 *RAL*DSTBR) ♦ (0 . 03865177*RABET*DSTBR) 
♦♦(0.00062918*DSTBR) 

C 

CNDRD*-0. 00153402+ (0.00184982*RAL)-(0.0068693*RAL*RAL) 

♦+ (0.01772037* (RAL**3) ) 

♦♦(0. 03263787* (RAL**4) )-( 0 . 15157163* (RAL**5 ))♦ (0 . 18562888 
♦*(RAL**6) ) -(0.0966163* (RAL**7) )♦ (0.01859168* (RAL**8) ) +(0.0002587 
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♦*RAL*DSTBR) -(0.00018546*RAL*DSTBR*RBETA) -(0.00000517304*RBETA) 

+-{0.00102718*R;X*RBETA) -(0.0000689379*RBETA*DSTBR) -(0.00040536 
♦*RBETA*RARUD) - (0 .00000480484*DSTBR*RARUD) 

+-( 0 . 00041786 ‘RAL'RARUD) 

++ (0.0000461872*RBETA)-c(0. 00434094*  (RBETA**2 ) ) 
+-(0.00490777*(RBETA**3) ) 

++ (0 . 000005157867*RARUD) ♦ (0 . 00225169*RARUD*RARUD) -( 0 . 00208072 
+*(RARUD**3) ) 

C 

IF  (RAL  .LT.  0.55851)  THEN 
C 

CMNP=-0. 00635409- (1.14153932*RAL) +(2. 82119027* (RAL**2) )♦ 

♦ (54.4739579* (RAL**3) )- (140. 89527667* (RAL**4) )- (676 .73746128* 

■►(RAL**5)  )  + (2059.18263976*  (RAL**6)  )♦  (1579.41664748*  (RAL**7)  ) 

♦- (8933. 08535712* (RAL**8) )+ (6806. 54761267* (RAL**9) ) 

ENDIF 

C 

IF  ((RAL  .GE.  0.55851001)  .AND.  (RAL  .LE.  0.61087))  THEN 
C 

CMNP=-.07023239-*-(1.085815*(RAL  -0.55851)  ) 

-f+(8.852651*(  (RAL-.  55851 )  **2)  )- (192 . 6093*  (  (RAL-0 . 55851)  **3) ) 
ENDIF 
C 

IF  (RAL  .GT.  0.61087)  THEN 
C 

CMNPa-71. 03693533+ (491. 32506715*RAL) 

+- (1388. 11177979* (RAL**2) )♦ 

►  (2033 .48621905* (RAL**3) ) 

+- (1590. 91322362* (RAL**4) )+ (567. 38432316* (RAL**5) ) 

+-( 44. 97702536* (RAL**7))+ (2. 8140669* (RAL**9) ) 

ENDIF 


IF  (RAL  .LE.  -.069813)  THEN 
C 

CMNR«  -0.2805 
ENDIF 
C 

IF  ((RAL  .GT.  -.069813)  .AND.  (RAL  .LT.  0.0))  THEN 
C 

CMNR.-O. 2804999948+ (35. 9903717041* (RAL+. 0698129982 )* *2) 
+  -( 516. 1574707031*  (R;iL+. 0698129982 )**3) 

ENDIF 

C 

IF  ((RAL  .GE.  0.0)  .AND.  (RAL  .LE.  0.78539801))  THEN 
C 

CMNR»-. 28071511- (2. 52183924*RAL) +(68. 90860031* (RAL**2) ) 
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+- (573.23100511* (RAL**3) )♦ (2009. 08725005* (RAL**4) ) 
+-(3385.15675307*(RAL**5) ) 

ff (2730. 49473149* (RAL**6) )- (848. 12322034* (RAL**7) ) 

ENDIF 

C 

IF  ((RAL  .GT.  0.78539801)  .AND.  (RAL  .LT.  0.95993102))  THEN 


C 

CMNR=-0. 1096954+ (0. 52893072* (RAL-0. 78539801) )-(6.09109497*(RAL- 
♦0.78539801) **2)+ (17. 47834015* (RAL-0.78539801) **3) 

ENDIF 

C 

IF  (RAL  .GE.  0.95993102)  THEN 
C 

CMNR=-0 .11 
ENDIF 
C 

CNDTDs0.00058286+(0.0007341*RAL)- (0.00746113 *RAL*RAL) 
♦-i0.00685223*(RAL**3) ) 

♦♦(0.03277271* (RAL**4) )-( 0 . 02791456* (RAL**5 ) ) 

♦♦ (0.00732915* (RAL**6)) 

♦♦(0.00120456*RAL*DSTBR)-(O.OC168102*DSTBR)+( 0.0006462* 
♦DSTBR*DSTBR) 

C 

CNDAD*0. 00008228887.(0. 00014015*RAL)-(0.0013493*RAL*RAL)+ 
♦(0.00020487* (RAL**3) )♦ (0.00561241* (RAL**4) ) 

♦- (0.00634392* (RAL**5) ) 

♦♦ (0.00193323* (RAL**6) )- (2. 05815E-17* (RAL*DAILA) )♦ (3 .794816E-17* 
♦  (DAi'LA**3)  ) 

C 

DCNB=-2.500E-4 

C 

RALNlsO. 69813 
RALN2= 9 0.0 /DEGRAD 
RBETNls-0. 174532 
RBETN2rO. 34906 
C 

AN*0.034 

ASTARN=1.0472 

BSTARN=0.08'^266 

C 

2ETAN= (2 . ODO-ASTARN- (RALN1+RALN2 ) ) / (RALN2 -RALNl ) 

ETAN=(2 .0D0*BSTARN'(RBETN1+RBETN2) ) / (RBETN2 -RBETNl ) 

C 

XN= ( 2 . ODO  *RAL - ( RALNl ♦RALN2 ) ) / ( RALN2 -RALNl ) 

YN= (2 .0D0*RBETA- ( RBETNl ♦RBETN2 ) ) / (RBETN2 -RBETNl ) 

C 

FNs( (5.0D0* (ZETAN**2) ) - (4 . ODO*ZETAN*XN) -1 .ODO) * 

♦( ( (XN**2) -1.0D0)**2)/( ( (ZETAN**2)-1.0D0)**3) 

C 


139 


o o o o o n n o n o  non  nnooo  ooo 


GN={ (5.0D0*(ETAN**2) )-{4.0D0*ETAN*YN)-1.0D0) • 

+( ( (YN**2) -l.ODO) **2)/{ ( (ETAN**2) -l.ODO) **3) 

C 

CNRBaAN*FN*GN 

C 

IF  (RAL  .LT.  0.69813)  THEN 
C 

CNRBxO.O 
GOTO  1000 
ENDIF 
C 

IF  ((RBETA  .LT.  -0.174532)  .OR.  (RBETA  .GT.  0.34906))  THEN 
C 

CNRBxO.O 
GOTO  1000 
ENDIF 


1000  CMN={CMNl*EPA02S)  +  (CNDAD*DAILD)  +  (  (CNDRD*DRUDD’'DRFLX3)  *EPA43)  + 

+  (  (CNDTD*DTFLX3)  *DTALD)-*-(CMNP*PB)  +  (CMNR*RB)  +  {DCNB*BETA) 
•i+CNRB 


CXxCFZ*SIN(RAL) -CFX*COS (RAL) +THRUST/QBARS 
CYxCFY 

C2=- (CFZ*COS (RAL) +CFX*SIN(RAL) ) 

CLMr-CML 

CMMxCMMf THRUST* ( 0 . 2 5 / 12 . 0 ) / ( QBARS *CWING ) 

THE  (0.25/12.0)  IS  THE  OFFSET  OF  THE  THRUST  VECTOR  FROM  THE  CG 
CNMxCMN 

RETURN  CX,  CY,  CZ,  CLM,  CMM,  CNM  TO  CALLING  PROGRAM. 


IWRITExO 

IF(IWRITE.EQ.l)  WRITE(6,*)'  CX, CY,CZ,CLM, CMM, CNM=  ' , 

+  CX,CY,CZ,CLM,CMM,CNM 

WRITE(6,*)  '  END  OF  SUB  COEFF  U(l,2, 3)  =  ' ,U(1) ,U(2) ,U(3) 
WRITE ( 6 , * ) ' LEAVING  COEFF ' 

RETURN 

END 
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AP»SNDIX  D 


T-38A  STRIP  CHART  DATA 


-38A  S/M  68-0205  20,000  ft  PA 
1  rj  Stall 
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Figure  29.  Sideslip,  Roll  and  Yaw  Rato  of  Wing  Rock 


143 


-38A  S/N  68-8205  20,000  ft  l»A 
70*  Uankcd  ‘rurn 
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Figure  3i .  Sideslip,  Roll  and  Yaw  Rate  of  Wing  Rock 


8A  S/N  63-0135  20.000  ft  PA 
30*  Hanked  Turn 


Elapsed  Time  (Sec) 

Figure  32  .  Fideslip.  Koll  and  Yaw  Rate  of  Wing  Rock 


APPSMDXZ  I 

r-15D  TABULATED  FLXQHT  TB8T  DATA 
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Tabltt  Bl 

F-15D  Wing  Rock  OiM«t  Su—ary 


MANEUVER 

WING  ROCK 
ONSET  AOA 
(UNITS) 

MACH 

1  C  STALL 

26 

0.32 

IQO  AOB  TURN 

26 

0.32 

30*  AOB  TURN 

26 

0.35 

70*  AOB  TURN 

26 

0.44 

SYMMETRIC 

PULLUP 

30 

o 

o 

B2 

r-lS  D  wing  Rock  Lialt  Cycl>«  "10  flf 11 


UNITS 

AOA 

WING  ROCK  LIMIT  CYCLES  -  DEG  AOB  EXCURSIONS  | 

NO  WING 

NOCK  LIMCT 
CYCLC 

OiSCKVKO 

.1-5° 

m 

m 

B 

_*10® 

♦.15* 

B 

24 

X 

25 

X 

X 

26 

XX 

27 

XX 

26 

X 

X 

B 

30 

a 

XXX 

X 

32 

X 

X 

33 

X 

X 

34 

X 

XX 

35 

X 

36 

X 

X 

37 

X 

X 

38 

X 

40 

X 

X  >  ONE  OBSERVED  WING  ROCK  LIMIT  CYCLE 
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Table  B3 

r-15D  wing  Rock  Llalt  Cycla*  •  10«  AOB  TORN 


UNITS 

AOA 

MO  MIMO 
•OCX 

LXMtT 

CYCU 

otmvcD 

WING  ROCK  LIMIT  CYCLES  -  1 

DEG  AOB  EXCURSIONS  | 

♦,.5® 

n 

bb 

O 

O 

♦1 

25 

X 

26 

XX 

27 

X 

XX 

29 

X 

31 

X 

32 

X 

X  «  ONE  OBSERVED  WING  ROCK  I. CYCLE 


Table  R4 


P-15D  Wing  Rock  Linit  Cycles  -  30o  AOB  TORN 


UNITS 

AOA 

NO  WINS 
•OCK  LIMIT 
CYCLC 
oMtmo 

WINi 

1 

3  ROCK  LIMIT  CYCLES 
DEG  AOB  EXCURSIONS 

- 

BB 

Bi 

o 

wm 

El 

♦,15® 

25 

X 

26 

X 

X 

27 

X 

X 

28 

X 

X 

29 

B 

32 

X 

X 

33 

B 

X 

38 

B 

X 

X  >  ONE  OBSERVED  WING  ROCK  LIMIT  CYCLE 
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TaJ»l«  15 

r-lSD  Wing  Rock  Limit  Cyclos  -  70*  AOB  TORM 


X  •  ONE  OBSERVED  WING  ROCK  LIMIT  CYCLE 


A»riiiDzx  r 

rLXOflT  nST  nCHHIQOBS 


ISO 


FLZOBT  TBST  TBCHMZQOS8 


0«n«ral 

The  aircraft  was  trinned  to  steady  level  flight  at 
20,000  feet  PA  with  29.92  set  in  the  altimeter  and  feet 
off  the  rudder  pedals.  Extreme  care  was  given  to 
eliminate  all  lateral  rolling  tendencies  and  ensure  the 
ball  was  centered  in  the  front  cockpit  turn  and  slip 
indicator.  The  roll  CAS  was  then  disengaged  and  the 
process  repeated.  These  actions  were  taken  to  guarantee 
any  subsequent  aircraft  rolling  motions  were  the  result 
of  only  aerodynamic  forces.  In  the  front  cockpit,  the 
Heads  Up  Display  was  placed  in  the  Attitude  Directional 
Indicator  (ADI)  mode  to  obtain  the  greatest  amount  of 
flight  information. 

1  g  stall 

The  aircraft  was  trimmed  at  20,000  ft  PA  and  25 
units  AOA  as  described  above.  The  throttles  and  pitch 
rate  were  then  smoothly  modulated  to  establish  a  1  kt/sec 
bleed  rate.  This  bleed  rate  was  then  held  until  30  units 
AOA,  This  condition  was  maintained  for  10  seconds  to 
observe  the  F-15D  flight  characteristics.  The  AOA  was 
then  increased  by  two  cockpit  unit  increments  and  the 
observation  process  repeated  until  the  onset  of  wing  rock 
was  observed.  After  the  onset  of  wing  rock,  the  AOA  was 
increased  in  one  to  two  unit  increments,  as  conditions 
allowed,  and  the  observation  process  repeated.  The  AOA 
was  incrementally  increased  in  this  manner  until  a 
central  AOA  of  40  units  was  achieved  or  the  bottom  of  the 
data  band  was  exited.  The  control  stick  was  centered  to 
eliminate  lateral  inputs.  At  each  stabilized  AOA  value, 
the  magnitude  of  the  wing  rock  was  measured  by  observing 
the  bank  angle  oscillations  on  the  ADI.  The  ADI  was  also 
used  to  determine  the  pitch  angles  at  which  wing  rock  was 
occurring.  The  wing  rock  period,  when  possible,  was  also 
measured  from  the  ADI.  The  front  cockpit  AOA  indicator 
was  used  to  measure  the  wing  rock  AOA  and  excursions. 
Airspeed  and  altitude  were  taken  from  the  front  cockpit 
main  instruments. 

10  Degree  Bank  Turn 

The  F-15D  was  trimmed  at  20,000  ft  PA  and  25  units 
AOA  in  stea^  level  flight.  A  10  degree  bank  turn  was 
then  established  with  aileron  and  the  control  stick 
centered  after  the  input.  The  AOA  was  then  increased  in 
the  same  manner  as  during  the  1  g  Stall  test  point  and 
the  aircraft  wing  rock  characteristics  observed.  The 
process  was  repeated  in  both  directions. 
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30  Dttgr**  Bank  Turn 

The  F-15D  was  trimmed  at  21,000  ft  PA  and  25  units 
AOA  in  steady  level  flight.  A  30  degree  bank  was 
established  with  aileron  and  the  control  stick 
subsequently  centered.  AOA  was  increased  as  described 
above  and  the  F-15D  wing  rock  characteristics  observed. 

The  test  was  repeated  in  both  directions. 


70  Degree  Bank  Turn 

The  aircraft  was  trimmed  at  23,000  ft  PA  above  the 
data  band  at  25  units  AOA.  A  70  degree  bank  turn  was 
established  and  AOA  allowed  to  increase.  Power  was 
modulated  to  minimize  sink  rate.  Because  of  the  sink 
rate  and  narrow  data  bank,  the  AOA  was  not  stopped  at 
intermittent  values,  but  rather  increased  up  to  the  onset 
of  wing  rock.  The  wing  rock  characteristics  of  the 
aircraft  were  then  observed  and  recorded  at  this  and 
higher  AOA  values .  The  process  was  terminated  upon 
exiting  the  data  band.  The  maneuver  was  then  repeated  in 
the  opposite  direction. 

SyBB&atric  Pullup 

To  complete  the  symmetric  pullup  maneuver,  the 
aircraft  was  trimmed  to  steady  level  flight  at 
23,000-24,000  ft  PA  and  MachsO.6.  This  corresponded  to 
approximately  275  KCAS.  The  aircraft  was  then  pulled  to 
a  set  AOA  and  airspeed  was  allowed  to  increase  to 
maintain  0.6  mach  until  the  aircraft  exited  the  data 
band.  Any  wing  rock  characteristics  were  recorded.  This 
maneuver  was  repeated  for  each  AOA  tested,  with  a  single 
AOA  captured  for  each  maneuver. 
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APPENDIX  G 


TEST  POINT  SUMMARY 
& 

DATA  TOLERANCE 
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Table 
Teat  Point 

01 

SUSBMXy 

TEST 

PT. 

TRIM 

CONDITION 

MANEUVER 

STABILIZED  AOA‘ 
(COCKPIT  UNITS) 

0 

BLEED 

RATE 

KT/SEC 

REMARK 

1 

20K.  25 

CPU  AOA 

1  C  Stall 

30.32.34,36,38,4 

0 

■ 

1 

Roll  CAS 
Off 

2 

20K.  25 

CPU  AOA 

10  deg 
banked  turn 

30,34,38 

■ 

BB 

Roll  CAS 
Off 

3 

20K.  25 

CPU  AOA 

30  deg 
banked  turn 

30,34.38 

B 

BB 

Roll  CAS 

Off 

1  * 

20K,  25 

CPU  AOA 

70  deg 
banked  turn 

30,34,38 

2.9 

BB 

Roll  CAS  1 
off  B 

1  ^ 

20K.  M-0.6 

Symmetric 

pullup 

N/A 

6 

BB 

Roll  CAS  1 

Off  1 

1.  b-ib  Will 

oe  leAiporaril^  stabilised  At  It 

idde  an^fids  ac 

accack  co  scudy  wing  rock  onset  and  limit  cycle  behavior.  T-38  was  be 
stabilized  at  different  angles  of  attack  that  was  determined  in  the 
initial  phases  of  this  program. 


Table  02 
Data  Tolarancaa 


1  PARAMETER 

TOLERANCE  | 

1  Altitude 

i  2000- 

1  AOA 

±  1  degree 

1  Bank  Angle 

2  degrees 

1  Mach  Number 

^0.02 
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